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PREFACE 


This  report  presents  the  results  of  a  study  effort  to  develop  an 
improved  technique  for  the  flight  measurement  of  aircraft  antenna  radia¬ 
tion  patterns  at  the  Air  Force  Flight  Test  Center  (T^FTC) ,  Edwards  AFB, 
California.  The  study  effort  and  tests  were  authorized  by  the  Aircraft 
Antenna  Pattern  Synthesis  Technique  Study  Plan^  February  1974,  and  the 
Test  Plan  ^or  Evaluating  Techniques  of  Determining  Aircraft  Antenna 
feadiation  t>atterns,  august  1974.  tfhe  tests  were  conducted  in  October 
1974  and  December  1974  under  Job  Order  Number  SC6341.  The  results  of 
the  tests,  which  were  used  to  show  the  validity  of  the  data  processing 
techniques,  are  shown  in  appendix  D  of  this  report. 

The  format  of  this  report  was  developed  to  make  this  report  more 
usable  to  Project  Engineers  of  the  .Systems  Engineering  Branch  at  AFFTC. 
As  such,  information  is  presented  to  give  a  novice,  in  the  field  of 
antenna  radiation  pattern  measurement,  sufficient  background  and 
knowledge  to  perform  an  accurate  evaluation  of  aircraft  antenna 
radiating  systems. 

The  author  wishes  to  acknowledge  the  following  individuals  who 
were  instrumental  in  the  preparation  of  this  report:  Mr.  B. 

Lyle  Schofield,  Chief,  Flight  Test  Technology  Branch,  for  guidance 
and  editorial  comments:  Mr.  Alfred  H.  Boyd,  Systems  Engineer,  for 
providing  technical  expertise;  Sgt  James  N.  Robertson,  Engineering 
Aide,  for  assisting  in  the  development  of  the  computer  software; 
and  our  secretaries  Mrs.  Dorothy  M,  Shaffer  and  Miss  Mary  Jane  Gugllotte 
for  their  tireless  efforts  in  transforming  this  report  into  its  final 
form. 
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INTRODUCTION 


The  dynamic  measurement  of  aircraft  antenna  radiation  patterns 
has  been  a  time  consuming  task  which  usually  produced  results  of  marginal 
utility.  The  classical  cloverleaf  flight  pattern,  used  in  the  past, 
produced  one  antenna  pattern  with  12  to  24  data  points  and  required 
1  1/2  to  3  hours  of  flying  time.  The  severe  limitation  of  the  minimal 
data  saunples  and  the  high  cost  of  flying  time  has  rendered  this 
technique  obsolete  and  created  a  requirement  for  an  improved  method. 

The  techniques  reported  here  were  developed  to  provide  the  Systems 
Engineer  with  the  tools  necessary  to  conduct  an  accurate  and  efficient 
evaluation  of  aircraft  antenna  systems.  To  accomplish  this,  a  semi- 
automated  system  was  dev^-loped  which  permits  the  aircraft  to  fly  a 
variety  of  flight  patterns  and  produces  continuous  360  degree  antenna 
patterns. 

The  concepts  of  dynamic  measurement  of  aircraft  antenna  patterns 
presented  in  this  report  are  the  results  of  a  survey  of  the  field  of  i 

antenna  pattern  measurement.  The  technique  developed  at  the  Air  Force 
Flight  Test  Center  (AFFTC)  drew  upon  the  concepts  being  used  by  Rome 
Air  Development  Center  (RADC) ,  Naval  Air  Test  Center  (NATO,  and  others. 

The  result  has  been  the  development  of  an  efficient  dynamic  aircraft 
2intenna  pattern  measurement  system. 

This  system  was  developed  and  tested  in  the  Fall  of  1974.  The 
initial  testing  was  accomplished  at  the  Precision  Antenna  Measurement 
System  (PAMS)  at  Griffiss  AFB,  Nev;  York.  PAMS  is  the  most  modern 
dynamic  antenna  measurement  facility  in  the  Air  Force  and  data  collected 
there  was  used  as  a  standard  for  comparing  techniques.  Further  testing 
was  accomplished  at  Edwards  AFB  using  a  C-131  and  an  F-15  aircraft.  The 
results  of  these  tests  are  shown  in  appendix  D. 

The  body  of  this  report  presents  information  which  can  be  used  by 
a  Systems  Engineer  to  conduct  an  efficient  and  accurate  evaluation  of 
an  aircraft  antenna  system.  The  report  is  arranged  such  that  it  can  be 
used  as  a  handbook  for  planning,  testing,  analyzing,  and  reporting  tests 
by  Systems  Test  Engineers. 

Appendix  B  and  C  of  this  report  provide  documentation  of  two  sets 
of  computer  programs.  These  progreuns  constitute  a  complete  data  reduction 
package  to  assist  in  the  planning  of  flights,  processing  of  radar 
position  tapes  and  signal  strength  tapes,  and  plotting  of  measured 
effective  radiated  power  (ERP)  of  each  received  signal. 


AIRCRAFT  ANTENNA  PATTERN  MEASUREMENT  TECHNIQUES 


There  are  several  techniques  available  to  measure  aircraft 
antenna  radiation  patterns.  Each  technique  has  advantages  and 
disadvantages  and  each  vary  in  accuracy,  cost,  complexity,  and  usage. 
This  section  will  present  an  introduction  to  each  of  the  four  basic 
techniques:  (1)  theoretical,  (2)  modeling,  (3)  dynamic,  and  (4)  near¬ 

field.  No  attempt  is  being  made  to  present  all  the  details  of  each 
technique,  but  only  to  present  the  basic  concepts  of  each  and  how 
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they  are  used.  If  a  more  detailed  understanding  of  each  of  these 
techniques  Is  desired  additional  material  referenced  in  the  bibliography 
may  be  used. 

THEORETICAL 

The  theoretical  calculation  of  ^’.!ie  radiation  pattern  of  an 
antenna  mounted  on  a  complex  structure,  an  aircraft  for  example,  has 
only  become  feasible  with  the  development  of  the  modern  high  speed 
computers  (reference  1).  Since  the  computation  involves  the  solution  of 
numerous  differential  equations,  this  technique  is  dependent  upon  the 
speed  of  the  computer.  Normally,  the  use  of  this  technique  has  been 
restricted  to  initial  design  stages  where  the  goal  is  to  locate  the 
antenna  in  the  best  position  which  will  produce  an  acceptable  radiation 
pattern . 

The  main  advantage  of  this  technique  is  that  the  radiation  pattern 
of  an  antenna  can  readily  be  evaluated  at  various  locations  on  the  air¬ 
craft  once  the  aircraft  structure  is  modeled  into  the  computer.  This 
permits  the  rapid,  low  cost  evaluation  of  aircraft  antenna  radiation 
patterns  and  allows  the  designer  to  locate  the  position  on  the  aircraft 
which  has  the  best  probability  of  producing  the  most  optimal  pattern. 

The  prime  disadvantage  of  this  technique  is  its  sensitivity  to 
errors.  Aircraft  VHF  and  UHF  antenna  radiation  patterns  can  be  very 
sensitive  to  small  dimensional  errors  in  modeling  parts  of  the  aircraft 
structure.  The  area  within  the  near-field  region  of  the  antenna  (within 
a  few  wave  lengths)  is  especially  susceptible  to  these  discrepancies. 
Since  the  aircraft  structure  is  only  approximated  in  the  computer 
program,  the  calculated  patterns  can  have  numerous  discrepancies. 

A  number  of  computer  programs  have  been  written  which  calculate 
the  theoretical  radiation  pattern  of  antennas  mounted  on  complex 
structures.  Each  program  uses  approximation  techniques  to  model  the 
shape  of  the  aircraft.  The  most  common  technique  is  the  wire-grid 
approximation  which  approximates  the  aircraft  structure  by  a  series 
of  wires.  Figure  1  shows  a  typical  wire-grid  approximation  model. 

A  second  modeling  technique  is  the  geometric  approximation  which  app::oxi- 
mates  the  aircraft  structure  with  a  series  of  simple  geometric  shapes 
(reference  2). 

Both  of  these  methods  only  provide  approximate  shapes  of  the 
aircraft's  actual  structvire;  therefore,  the  accuracy  of  the  final 
result  will  depend  upon  the  significant  deviation  of  the  approximated 
aircraft  structure  from  the  actual  structure.  Hence,  the  accuracy 
of  the  prediction  will  depend  upon  the  user's  ability  to  include  all 
significant  portions  of  the  aircraft's  structure  in  the  model. 

Because  these  approximations  can  never  eliminate  all  errors  in  the 
shape  of  the  structure,  theoretical  patterns  should  be  considered  as 
producing  an  approximate  shape  of  the  actual  antenna  pattern. 

MODELING 

The  modeling  method  is  the  technique  most  often  used  to  measure 
ancenna  radiation  patterns.  This  technique  uses  a  model  of  the  aircraft 
mounted  on  a  pedestal  which  rotates  the  model  around  its  axis  while  the 
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Wire-grid  Model  of  Helicopter 


antenna  patterns  are  measured.  The  models  can  be  of  any  size  and  are 
normally  mounted  in  an  anechoic  chamber  or  an  outdoor  range.  RF 
absorbing  material  will  cover  all  surfaces  v;hich  may  cause  ur.v.’anted 
reflected  signals.  The  largest  of  these  facilities  is  located  at  RADC's 
Stockbridge  Range,  where  an  actual  B-52G  is  mounted  and  rotated  on  a 
pedestal  (figure  2). 

With  the  modeling  technique,  the  transmitter  is  normally  mounted 
within  the  aircraft  and  the  frequency  increased  by  1/N,  where  N  is 
the  fraction  of  the  model  size.  Therefore,  a  1/4  scale  model  will  require 
the  frequency  to  be  increased  by  a  factor  of  four  to  create  an 
equivalent  pattern.  To  obtain  an  accurate  model  the  conductivity  and 
permittivity  of  the  model's  surface  should  be  increased  bv  the  Seune 
amount.  In  an  effort  to  accomplish  this,  the  surface  of  the  model  is 
normally  coated  \;ith  copper  which  \;ill  increase  the  conductivity  as 
much  as  economically  feasible  (reference  3). 

The  cost  and  the  ease  of  measuring  antenna  patterns  are  the  main 
advantages  of  the  modelinq  technique.  Since  the  model  is  mounted  on 
a  pedestal,  all  depression  angles  can  easily  be  measured  at  a  relatively 
low  cost  per  pattern. 

The  most  accurate  antenna  patterns  are  produced  when  a  full  scale 
model  is  used  and  care  is  taken  to  avoid  the  effects  of  the  pedestal 
and  reflected  signals.  When  scaled  models  are  used,  the  accuracy  of 
the  patterns  v;ill  depend  upon  tlie  precision  of  the  model.  Small 
dimensional  errors  in  the  model  can  produce  major  errors  in  the 
measured  patterns.  Addit*  anally,  to  create  valid  patterns  careful 
consideration  must  be  g'  -en  to  effects  of  supporting  pedestal,  cables, 
and  other  equipment  w>  jh  are  not  part  of  the  actual  aircraft.  Due 
to  these  problems,  mo,iel  patterns  mav  contain  holes  v;hich  do  not  exist 
in  the  actual  antenna  system. 

DYNAMIC 

Dynamic  measurements  of  aircraft  antenna  radiation  patterns  are 
the  final  checks  of  tiie  antenna  performance.  VJith  flight  evaluation, 
the  aircraft  is  placed  in  its  operational  environment  and  is  completely 
isolated  from  external  structure,  like  supporting  pedestals  and  cables. 
Ideally  the  aircraft  should  le  placed  in  space  and  rotated  about  its 
axis  while  a  receiver  on  the  ground  records  the  signal  strength  at 
all  depression  and  azimutli  angles.  Since  this  is  impossible,  the 
aircraft  is  flown  through  a  variety  of  flight  paths  such  that  the  signal 
strength,  at  all  depression  and  azimuth  angles  of  interest,  is  measured. 
Normally  only  the  Inv/er  hemispherical  pattern  is  measured  with  this 
technique,  since  tlie  aircraft  v.’ould  be  reraiired  to  fly  inverted  to 
measure  the  upper  hemisphere  from  the  ground.  A  variety  of  flight 
paths  which  can  be  used  is  presented  in  detail  in  a  latter  section 
in  this  report. 

NEAR -FI ELD 

The  near-field  technique  of  measuring  antenna  radiation  patterns 
has  been  one  of  the  ne\’est  developments.  This  technioue,  which 
basically  consists  of  automatic  scanning  of  the  antenna  v'ithin  the 
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Note:  The  B-52  is  in  the  inverted  position  tor  lower  hemispherical 
antenna  pattern  measurement. 


Figure  2  Stockbridge  B-52G  Range 
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near-field  region  and  extrapolating  the  data  to  the  far-field,  was 
developed  at  the  National  Bureau  of  Standards  by  D.  M.  Kerns  (reference 
4).  This  method  has  been  used  at  microwave  frequency  range  producing 
excellent  results.  Figure  3  shows  a  typical  test  in  an  anechoic  cheunber. 
Concept ionally,  this  technique  can  be  applied  to  VHF/UHF  frequencies. 
Therefore,  it  may  be  possible  to  apply  this  technique  to  measure 
some  types  of  aircraft  antenna  radiation  patterns  in  the  future. 

One  advantage  of  this  technique  is  that  the  ground  reflections 
and  incidence  reflections  in  anechoic  cheunbers  are  eliminated  because 
of  the  close  proximity  of  the  measuring  probe  to  the  'antenna.  Where 
there  is  strong  atmospheric  attenuation,  such  as  in  the  millimeter- 
wave  range,  this  technique  is  very  useful  since  the  far-field  is  determined 
from  near-field  measurements. 

The  main  disadvantage  of  the  technique,  in  measuring  aircraft  VHF/ 

UHF  antenna  pattern,  is  the  location  of  the  probing  antenna.  Since 
the  aircraft  structure  will  alter  the  antenna  pattern,  the  probing 
antenna  must  be  located  beyond  the  outer  structure. 


MULTIPATH  RADIATION 


Multipath  signals  are  sources  of  some  of  the  most  significant 
errors  when  dynamically  measuring  aircraft  antenna  radiation  patterns 
in  the  VHF/UHF  range.  These  errors  result  from  the  reception  of  a 
direct  signal  and  reflected  signal  simultaneously,  with  the  reflected 
signal  changing  phase  due  to  the  difference  in  path  length  and  reflection. 
The  reception  of  the  reflected  signal  is  capable  of  producing  changes 
in  signal  strength  of  over  18  db.  Due  to  the  magnitude  of  this  effect, 
care  must  be  taken  to  avoid  influencing  the  antenna  pattern  with 
multipath  radiation.  Therefore,  this  section  is  included  to  give  an 
understanding  of  the  considerations  required  to  conduct  flight  evalu¬ 
ations  of  antenna  patterns  in  the  presence  of  multipath  radiation. 

In  the  study  of  multipath  radiation,  the  desired  result  is  the 
magnitude  of  the  vector  sum  of  the  direct  signal  and  all  the  reflected 
signals  received  at  the  antenna.  If  both  the  receiver  and  transmitter 
were  located  in  free  space,  only  the  direct  signal  would  be  received 
with  its  strength  reduced  by  the  amount  of  spherical  dispersion.  When 
the  earth  is  placed  in  this  environment,  the  received  signal  is  altered. 
The  direct  ray  will  become  curved  due  to  atmospheric  refraction  and 
multiple  signals  v;ill  be  received  due  to  surface  reflections.  The 
resultant  signal  may  show  a  loss  or  a  gain,  due  to  the  phase  differences 
between  the  direct  and  reflected  signals. 

Since  a  smooth  earth  and  true  standard  atmosphere  do  not  exist, 
the  simple  two  path  propagation  model  is  modified  by  several  factors. 

Reed  and  Russell  (reference  5)  list  a  number  of  factors  which  tend  to 
cause  variations  in  multipath  propagation.  These  include: 

1.  Irregular  terrain,  such  as  hills  and  valleys,  which  cause 

reflected  signals  1/2  wave  length  out  of  phase  with  the  direct 
ray. 
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2.  Absorbing  or  reflecting  objects  located  near  the  groxind  antenna, 

i.e.,  trees,  buildings,  etc. 

3.  Aircraft  structure  located  near  the  antenna  which  can  cause 
absorption  or  reflection  of  radiated  signal. 

4.  Non-standard  atmosphere. 

5.  Atmospheric  turbulence. 

6.  Changes  in  reflecting  surface  conductivity  and  permittivity  which 
cause  changes  in  reflection  coefficient. 

To  effectively  use  *the  method  of  calculating  effects  of  multipath 
propagation  presented  in  appendix  E,  care  must  be  taken  to  minimize 
these  factors.  The  following  are  a  number  of  the  precautions  which 
must  be  taken  for  the  accurate  prediction  of  the  effects  of  two-path 
propagation. 

1.  The  receiving  antenna  must  be  isolated  from  structures  which  may 
cause  reflections  of  signals  to  the  antenna  or  interfere  with 
the  reflected  signal  being  calculated.  Some  objects  which  may 
cause  this  are  locating  the  antenna  near  a  roof  or  locating  the 
antenna  near  buildings. 

2.  In  order  to  minimize  the  effects  of  non-scandard  atmosphere  and 
atmospheric  turbulence,  the  flights  should  be  planned  during 
periods  of  sta>)le  atmospheric  conditions.  This  will  mean  to 
avoid  flying  during  the  hot  summer  afternoons  when  rapidly  rising 
thermals  are  present. 

3.  Changes  in  conductivity  and  permittivity  of  the  reflecting  surface 
can  have  a  pronounced  effect  on  the  reflected  signal.  To  minimize 
these  effects  the  flights  should  be  planned  with  constant  surface 
conditions.  This  will  mainly  require  flights  to  be  made  only 
when  the  soil  is  dry. 

SENSITIVITY  TO  ERRORS 

When  the  technique  presented  in  appendix  E  is  used  to  predict  the 
effect  of  multipath  radiation  on  the  received  signal,  the  effects  of 
changes  in  the  input  parameters  should  be  understood.  Figures  4  through 
7  show  the  effects  of  c'r.an'fing  various  parameters  as  predicted  by  the 
multipath  option  of  the  ^‘.n'^enna  Data  Analysis  Progrzun  (ADAP)  described 
in  appendix  C.  The  param.riers  changed  in  these  figures  are: 

(a)  ground  antenna  height,  (b)  aircraft  height,  (c)  permittivity  of  the 
reflecting  surface,  and  (d)  conductivity  of  the  reflecting  surface. 

There  are  numerous  other  factors  which  can  affect  the  actual  reception 
of  multipath  radiation.  These  include:  surface  roughness,  turbulent 
atmosphere,  and  multiple  reflects.  Through  careful  flight  planning 
these  effects  can  be  minimized. 

The  importance  of  each  parameter  in  contributing  to  errors  in  the 
calculated  multipath  effect  is  a  function  of  the  ability  to  accurately 
measure  the  pareuneter  and  the  effect  the  parameter  has  on  the  calculation. 
Therefore,  if  changes  in  a  parameter  can  have  major  effects  on  the  pre¬ 
dicted  multipath  radiation,  but  the  parameter  is  easily  measured  with  a 
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high  degree  of  accuracy,  the  parameter  will  be  relatively  unimportant 
for  error  consideration.  Frequency  is  one  such  parameter  since  it  can  be 
accurately  measured,  is  very  stable,  and  has  a  major  effect  on  the  cal¬ 
culated  value. 

The  two  main  parameters,  which  produce  the  majority  of  the  errors  in 
the  prediction  of  the  multipath  effect,  are  the  height  of  the  antenna 
ad>ove  the  reflecting  surface  and  the  permittivity  of  the  reflecting 
surface.  The  reason  for  their  importance  is  that  they  are  difficult 
to  accurately  measure  and  they  have  major  effects  on  the  prediction  of 
multipath  radiation. 

The  effect  of  errors  in  the  height  of  the  ground  antenna  can  be 
seen  in  figure  4.  This  figure  shows  that  with  a  change  in  the  height 
of  the  antenna  above  the  reflecting  surface  of  eight  feet  will  produce 
a  significant  change  in  the  predicted  lobe  structure.  An  increase 
in  the  height  of  the  antenna  will  increase  the  number  of  lobes  and 
reduce  the  distance  betv/een  each  lobe.  These  errors  can  be  caused  by 
both  errors  in  measuring  the  antenna  height  and  changes  in  terrain. 

Since  the  terrain  is  not  normally  flat,  this  probably  is  the  major  cause 
of  errors  in  the  height  of  the  antenna. 

The  effect  of  errors  in  the  height  of  the  aircraft  above  the  reflect¬ 
ing  surface  can  normally  be  neglected.  Figure  5  shows  the  predicted 
effects  of  multipath  radiation  at  two  different  aircraft  altitudes. 

The  effect  is  to  increase  the  distance  between  lobes  as  the  aircraft 
distance  above  the  reflecting  surface  is  Increased.  Since  the  air¬ 
craft  altitude  can  be  easily  measured  to  within  +500  feet,  these  errors 
normally  have  little  effect  on  the  accuracy  of  tHe  predicted  multipath 
gain. 


The  effects  of  changing  the  permittivity  can  be  seen  in  figure 
6.  This  figure  shows  that  permittivity  affects  the  predicted  gain 
but  has  very  little  effect  on  the  location  of  the  lobes.  Additionally, 
permittivity  has  a  very  large  range  of  possible  values.  For  soil  the 
relative  permittivity  has  a  range  of  2  to  30. 

Conductivity  of  the  reflecting  surface  is  the  one  parameter  which 
affects  both  the  predicted  gain  and  the  location  of  the  lobes.  Figure 
7  shows  the  effects  of  changes  in  conductivity  on  the  prediction  of 
multipath  signals.  This  shows  that  only  for  high  conductivities 
does  the  change  in  conductivity  have  any  significant  effect  on  the  pre¬ 
dicted  gain.  Since  earth  normally  has  a  conductivity  between  0.1  and 
0.00001,  the  effects  of  changes  in  conductivity  can  be  minimized.  This 
can  be  accomplished  by  only  flying  when  the  surface  is  dry. 

The  use  of  Rogers  Dry  Lake  at  Edvrards  AFB  for  a  reflecting  surface 
appears  to  be  Ideally  suited  for  use  of  predicted  multipath  radiation 
effects.  The  surface  is  relatively  flat  and,  due  to  lack  of  moisture, 
has  constant  permittivity  and  conductivity  throughout  most  of  the  year. 
This  enables  the  direct  signal  strength  to  be  determined  to  within 
+2  db. 


In  case  the  flight  must  be  conducted  when  moisture  is  present  on 
Rogers  Dry  Lake,  the  conductivity  must  be  adjusted  to  produce  the 
correct  lobe  structure.  This  can  best  be  accomplished  by  flying 
outbound  and  Inbound  radials  to  measure  the  multipath  effect  and  use 
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Figure  4  Efiect  of  Antenna  Height  Errors  on  Multipath 
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Figure  5  Effect  of  Aircraft  Altitude  Errors  on  Multipath 


igure  6  Effect  of  Permittivity  Errors  on  Multipath 


the  MULT  option  of  ADAP  (appendix  C)  to  predict  the  multipath  structure 
at  different  conductivities.  The  conductivity  which  produces  the  best 
match  of  the  measured  nulls  and  lobes  should  be  used  in  the  reduction 
of  the  data  by  the  Antenna  Radiation  Pattern  Measurement  Program  (ARPMP 
appendix  B) . 

USE  OF  COMPUTER  PROGRAM 

The  technique  of  calculating  the  effects  of  multipath,  as  presented 
in  appendix  E,  is  available  in  the  computer  program  ADAP  (see  appendix 
C).  As  a  part  of  ADAP  the  computer  program  MULT  will  calculate  the  multi- 
path  effect  at  a  given  altitude,  ground  antenna  height,  frequency,  con¬ 
ductivity,  and  relative  permittivity.  The  program  creates  line  printer 
plots  plus  a  listing  of  all  data  points  calculated  (figure  8). 

This  program  can  be  useful  for  planning  flight  paths  for  the 
measurement  of  aircraft  antenna  patterns.  By  predicting  the  lobe 
structure,  the  flight  path  can  be  determined  to  keep  the  aircraft 
within  an  area  of  reinforced  signals. 

Figure  9  shows  the  effect  of  multipath  radiation  on  a  radial  flown 
at  7500  feet  AGL.  Figure  8  shows  the  predicted  change  in  the  signal 
due  to  multipath  and  figure  10  shows  the  radial  corrected  for  multipath 
radiation.  Vfhen  flight  evaluations  of  antennas  are  conducted,  the 
patterns  flown  must  avoid  the  null  areas  to  produce  meaningful  results. 
Even  though  the  Antenna  Radiation  Pattern  Measurement  Program  (ARPMP 
appendix  B)  has  the  option  of  calculating  the  multipath  effect,  accuracy 
is  lost  when  data  is  collected  within  the  areas  of  deep  nulls.  Therefore, 
by  using  program  MULT,  the  locations  of  deep  nulls  can  be  predicted  and 
the  flight  pattern  designed  to  keep  the  aircraft  out  of  the  null  areas. 


FLIGHT  TECHNIQUES 


The  key  to  efficient  dynamic  measurement  of  aircraft  antenna 
radiation  patterns  is  the  selection  of  an  appropriate  flight  pattern 
which  can  produce  the  data  desired  Ir.  *  r'l'^imum  of  flight  time.  To 
accomplish  thli^,  a  number  of  flight  techniques  are  available.  Depending 
on  the  capability  of  the  aircraft  and  the  type  of  antenna  patterns  to  be 
measured,  one  of  the  flight  techniques  listed  in  this  section  will  produce 
the  best  antenna  pattern  in  the  shortest  flight  time.  For  most  cases 
where  360  degrees  of  aspect  angle  are  to  be  measured  at  0  to  -5  degrees 
depression  angle,  the  skidding  turn  flight  pattern  should  be  used. 

CLOVERLEAF 

The  cloverleaf  flight  pattern  is  the  typical  pattern  used  in  the 
past  to  measure  aircraft  antenna  patterns.  The  majority  of  these  patterns 
are  12,  18,  or  24  point  cloverleafs  (figure  11).  The  pattern  should  be 
flown  at  a  range  and  altitude  required  to  measure  the  depression  angle  of 
Interest  and  at  the  same  time  keep  tlie  center  of  the  pattern  within  a 
major  lobe  of  the  ground  antenna  pattern.  This  will  Insure  the  strongest 
signal  is  received  and  will  minimize  the  errors  due  to  changes  in 
multipath  radiation.  The  normal  method  used  to  accomplish  thi£.  is  to 
pick  the  range  and  altitiide  first  and  then  adjust  the  height  of  the 
receiving  antenna  to  place  a  major  lobe  at  that  point. 
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FIGURE  9  riEfiSUREO  MULTIPRTH  EFFECT 


FIGURE  10  RflOIflL  CORRECTED  FOR  MULTIPATH 


Once  the  location  for  cloverleaf  is  determined,  the  aircraft  can 
be  either  radar  vectored  over  the  pxsint  at  the  proper  heading  or  the 
pilot  can  maneuver  the  aircraft  over  the  point  by  using  a  predominant 
landmark  as  the  center  point.  If  radar  tracking  is  not  used  the  aircraft 
should  pass  within  1/2  NM  of  the  center  to  produce  accurate  data.  An 
error  of  1  Nfl  can  produce  1  db  change  in  the  received  signal  as  a  result 
of  range  and  multipath  signal.  Therefore,  to  insure  the  data  is  not 
affected  by  changes  in  range  and  multipath,  the  signal  strength  should 
only  be  measured  v;hen  the  aircraft  is  within  1/2  NM  of  the  cloverleaf 
center. 

The  major  advantage  of  this  pattern  is  that  the  signal  strength  is 
measured  when  the  aircraft  is  over  the  same  point  each  time.  This  will 
tend  to  minimize  the  errors  caused  by  multipath  signals  and  atmospheric 
attenuation.  Another  advantage  of  this  pattern  is  that  the  time  required 
to  maneuver  the  aircraft  from  each  data  point  permits  the  data  to  be 
recorded  with  a  manual  system.  Consequently,  an  automated  system  is 
not  required  for  this  flight  technique. 

The  cloverleaf  has  two  major  disadvantages:  (1)  flight  time 
required,  and  (2)  quality  of  antenna  patterns  produced.  The  flight  time 
required  for  a  24  point  cloverleaf  can  vary  from  1  1/2  to  3  hours 
depending  upon  the  type  of  aircraft  flown.  Since  only  24  data  points 
are  measured,  the  vast  majority  of  the  flight  is  required  solely  for 
maneuvering  the  aircraft  and  not  collecting  data.  This  is  reflected 
in  the  quality  of  the  antenna  patterns  produced.  Figures  12  and  13 
show  two  antenna  patterns  measured  using  an  18  point  cloverleaf  and  a 
skidding  turn.  The  severe  limitation  of  this  flight  technique  can  be 
seen  by  comparing  the  tv;o  antenna  patterns.  Major  holes  can  exist 
in  the  actual  antenna  pattern  and  would  not  be  detected  using  the 
cloverleaf  technique.  Therefore,  this  flight  technique  is  of  only 
limited  use. 

PARALLEL  FLYBYS 

The  parallel  flyby  pattern  consists  of  a  series  of  offset  radials 
flown  past  the  ground  antenna  (figure  14).  The  length  of  each  leg  and 
the  distance  of  each  offset  is  determined  such  that  all  depression 
angles  and  aspect  angles  of  interest  are  measured. 

The  legs  of  a  parallel  flyby  pattern  measure  depression  angles  of 
the  decaying  sine  wave  form  (figure  15) .  As  the  altitude  of  the  aircraft 
is  lowered,  the  depression  angles  measured  will  approach  the  horizontal. 

When  flying  this  type  of  pattern  the  radar  controller  must  be  able 
to  visually  extrapolate  the  aircraft  course  from  a  small  flight  segment. 
This  will  be  required  to  vector  the  aircraft  to  the  proper  leg  and  keep 
it  from  drifting  off  course.  To  insure  the  aircraft  will  at  least 
parallel  the  plotted  course,  specific  headings  should  be  called  to  the 
aircraft  rather  than  just  the  changes  in  headings.  All  aircraft  heading 
changes  made  after  the  start  of  the  leg  must  be  recorded  with  the  IRIG 
"B"  time  for  use  in  the  master  data  reduction  program  (ARPMP) . 

The  most  critical  portion  of  the  flight  to  keep  the  aircraft  on 
course  is  the  legs  within  10  NM  from  the  ground  antenna.  Slight  errors 
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FIGURE  12  CLOVERLERF  RNTENNR  PATTERN 


Figure  14  Parallel  Flyby  Flight  Pattern 


in  the  course  on  these  legs  will  produce  major  errors  in  the  aspect 
and  depression  angles  measured. 

The  major  advantage  of  this  type  of  flight  pattern  is  that  the 
complete  lower  hemisphere  can  be  measured  with  a  high  number  of  data 
samples.  Therefore,  a  high  degree  of  accuracy  can  be  obtained  for 
effectively  all  the  depression  angles  in  the  lower  hemisphere  of  the 
aircraft  in  tv;o  to  four  hours  of  flight  time. 

The  major  disadvantage  of  this  flight  technique  is  that  for  valid 
data  to  be  collected,  the  multipath  signal  must  be  eliminated.  The 
best  method  to  accomplish  this  is  to  use  high  directional  receiving 
antennas.  This  is  satisfactory  for  frequencies  above  1000  MHz  in 
which  a  six  foot  diameter  dish  will  produce  a  sufficiently  narrow  beam. 
Ilov/ever,  for  frequencies  between  100  MHz  and  300  MHz  the  antenna  dish 
would  have  to  be  up  to  ten  times  the  size  of  the  1000  MHz  antenna. 

Since  an  antenna  this  large  is  impractical  to  use,  the  multipath  signal 
cannot  be  eliminated  by  the  receiving  antenna. 

An  alternate  approach  is  to  use  the  multipath  prediction  pro¬ 
gram  to  calculate  the  magnitude  of  the  effect.  This  method  of 
eliminating  multipath  signals  is  only  usable  for  reflection  angles 
of  less  than  12  degrees.  Approximations  used  in  deriving  the  equations 
makes  the  calculations  inaccurate  above  12  degrees.  Therefore,  the 
elimination  of  the  multipath  signal  by  prediction  can  only  be  used  for 
near  0  degrees  depression  angle  measurement. 

MODIFIED  PARALLEL  FLYBYS 

The  modified  parallel  flyby  technique  is  essentially  a  parallel 
flyby  except  the  flight  path  has  been  altered  to  keep  the  aircraft  within 
a  major  lobe  of  the  ground  antenna  pattern  (figure  16).  This  flight 
technique  is  designed  to  measure  a  specific  depression  angle  where  the 
parallel  flyby  technique  will  measure  effectively  the  complete  lower 
hemisphere.  The  modified  parallel  flyby  is  conducted  in  the  same  manner 
as  the  parallel  flyby  technique.  The  aircraft  should  be  radar  vectored 
to  the  beginning  of  each  leg  allowing  sufficient  time  to  Insure  the 
aircraft  is  on  course  before  the  start  of  the  leg.  To  insure  the  aircraft 
is  at  least  parallel  to  the  correct  course,  the  pilot  should  be  given 
specific  headings  rather  than  just  the  heading  change  when  being  radar 
vectored . 

The  advantage  of  this  type  of  flight  pattern  is  that  it  has  the  seune 
high  number  of  data  samples  as  the  parallel  flyby  technique;  yet,  the 
data  is  only  collected  in  the  reinforced  regions  of  the  ground  antenna 
pattern.  Therefore,  the  data  collected  is  not  biased  by  major  nulls. 

An  additional  advantage  of  this  technique  is  that  the  airspace  required 
to  fly  the  pattern  is  reduced. 

The  major  disadvantage  is  the  flight  time  required  to  complete  the 
pattern.  To  create  a  continuous  360  degree  antenna  pattern  up  to  32 
individual  legs  must  be  flown.  Another  disadvantage  is  the  difficulty 
to  insure  data  samples  are  obtained  at  all  aspect  angles.  Since  each 
leg  produces  a  small  Interval  in  aspect  angle,  a  small  error  in  the 
aircraft's  path  can  cause  specific  aspect  angles  to  be  missed. 
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Figure  16  Modified  Parallel  Flyby  Flight  Pattern 


SKIDDING  TURN 


The  flight  technique  which  is  capable  of  producing  a  360  degree, 
constant  depression  angle,  antenna  pattern  in  the  shortest  time  is  the 
skidding  turn.  This  technique  requires  the  aircraft  to  complete  a 
360  degree  turn  keeping  the  wings  level.  To  assist  in  making  the  turn, 
rudder  trim  and  asymmetric  pov/er  may  be  used.  The  location  of  the  center 
of  the  turn  should  be  picked  to  keep  the  aircraft  within  a  major  lobe 
of  the  ground  antenna  (figure  17).  As  with  the  other  techniques,  the 
areas  of  cancellation  should  be  avoided  to  minimize  the  multipath 
radiation  effects.  The  dleuneter  of  the  circle  will  depend  on  the  capa¬ 
bilities  of  the  aircraft.  On  previous  tests  a  C-131  flew  the  skidding 
turn  within  an  8  NM  circle  and  an  F-15  within  5  NM.  Since  the  objective 
of  this  flight  technique  is  to  rotate  the  aircraft  about  its  vertical 
axis,  the  actual  ground,  path  is  of  no  concern  except  to  insure  the 
aircraft  does  not  enter  an  area  of  signal  cancellation.  To  correct 
for  distortion  of  the  measured  antenna  pattern,  the  multipath  correction 
option  should  be  used  when  processing  the  data.  Figures  18  through  21 
show  a  comparison  of  skidding  turn  data  with  and  without  the  multipath 
correction  option. 

To  properly  conduct  a  skidding  turn  flight,  a  map  to  be  used  for 
radar  tracking  should  be  marked  with  concentric  circles  around  the  desired 
orbit  center.  By  using  a  one  inch  to  1  NM  scale  map  and  2.5  NM,  5  NM, 

10  NM,  15  NM  diameter  circles,  the  aircraft  can  be  vectored  to  a  starting 
point  which  will  minimize  the  effects  of  wind  drift.  The  major  concern 
is  to  keep  the  aircraft  from  drifting  into  a  region  of  cancellation. 

After  the  aircraft  is  vectored  to  the  starting  point,  the  aircraft  will 
be  setup  into  a  constant  rate  skidding  turn.  Every  15  degrees  during 
the  turn,  the  pilot  must  call  his  heading  as  read  off  the  directional 
gyro.  The  heading  and  IRIG  "B"  time  will  be  recorded  at  the  ground 
station  for  use  in  the  Antenna  Radiation  Pattern  Measurement  Program. 

This  progr2un  performs  a  linear  time  extrapolation  to  determine 
the  aircraft  heading  throughout  the  turn  for  correlation  with  radar  and 
signal  strength  data.  With  the  aircraft  in  a  constant  rate  turn,  the 
headings  can  be  determined  within  two  degrees  using  this  technique.  The 
major  source  of  error  is  in  reading  the  directional  gyro  and  recording 
the  proper  time. 

If  the  aircraft  cannot  complete  the  turn  within  a  10  NM  diameter 
circle  or  if  the  aircraft  drifts  too  close  to  the  ground  station,  the 
circle  can  be  flov.m  in  segments.  Figures  20  and  21  show  data  collected 
on  an  F-15  turn  flown  in  two  segments  with  and  v.’ithout  multipath 
correction. 

There  are  numerous  advantages  to  this  t  of  flight  pattern.  The 
main  advantage  is  the  short  time  required  to  ;„4'-\'nplete  a  360  degree 
antenna  pattern.  Depending  on  the  capabilities  of  the  aircraft,  this 
turn  can  take  as  little  as  four  minutes.  Since  the  time  required  to 
complete  a  pattern  is  so  short,  drift  will  be  minimized  in  the  test 
equipment . 

The  major  disadvantage  is  the  requirement  for  the  aircraft  to  be 
forced  into  the  turn  with  s.ideslip.  The  ability  of  the  aircraft  to 
accomplish  this  depends  on  the  type  of  aircraft.  In  some  aircraft,  this 
may  cause  engine  fuel  starvation  due  to  lateral  acceleration. 


37 


PQjeR  iAIN  PL3T  *I  ALIITUO?  3F 
7300.000  FT.  ASL. 


I 


I 

i 

I 

I 


O 

en  o 

o  • 


o 


CD 


O 


o 


o 


o 


o 


I  o 


o 


I  u> 
♦  <»• 


o 


I  o 

♦ 


o 


I 

♦  fo 


<0 

c;  • 

3 


>  X 
mO  < 
»-o  r 

O'- 

SC) 

QO  CM 

r-  •  ro 
a  s  • 
o  r 

z 

V 

►-0  r 

•-tO 

>  •  • 

»-  Z 


u 

UJ 


) 


Z 

< 

O 

o. 

Ui 

3 

a 

a. 


-f 


I  O  ti. 

I  •sc  »0  UJ 

I  OZ  XO  QL 

♦  ■  • 

I  lu  o<-»  n 
I  o  —  ►- 

I  Z  ~l 

I  <  o  < 

I  QC  Ql  Z 

I  r 

I  O  liJ 

I  •  fsj  ►- 
I  X  Z 

♦  CM  r  •? 

I  o 

I  v*!)  • 

I  O  •  ry 

I  zc^  u. 

I  ujr> 

I  r>-*  UJ 

10  o  o 
I  •  UI  z 
10  a.  < 

♦  ►- 

I  "0 


C3 


U.  O 
UJ 

►-  ^ 

T  < 
t*>  3 

*-•  r?”* 

ujo  m 
T*  •  CP 

r»-  ZfT* 
•jfo  ncv’ 
z  •-•«-• 
z  • 

tlJ  O  I 
»-  < 
z  cy 
<  u. 
u 
A 

C-Cj 
»  C.D 

z  o 

XO 
«J  U.  • 

«  •  oir» 

z 

(tir>  »- 

UJ 

►-  c 

z  UJ 


o 


31 


Figure  17  Location  of  Skidding  Turn 


FIGURE  21  SKIDDING  TURN  RNTENNfl  PRTTERN 


BANKED  TURNS 


Banked  turns  are  flown  using  the  seune  basic  technique  as  the 
skidding  turns,  except  the  aircraft  is  placed  in  a  constant  bank  turn. 

Once  the  turn  is  established  the  pilot  must  call  the  heading  at  the 
start  and  every  15  degrees  throughout  the  turn.  The  headings  and  IRIG  "B" 
times  should  be  recorded  at  the  ground  station  for  use  in  the  Antenna 
Radiation  Pattern  Measurement  Program.  The  seune  precautions  to  avoid  the 
areas  of  signal  cancellation  should  be  taken  with  this  technique  as 
with  the  skidding  turns. 

Banked  turn  orbits  produce  a  series  of  decaying  sine  wave  depression 
angles  as  shown  in  figure  22.  Through  the  combination  of  a  series  of 
different  banked  turns,  partial  upper  and  lower  hemispheric  patterns  can 
be  measured.  The  actual  depression  angles  covered  are  a  function  of 
bank  angle,  range,  altitude,  and  orbit  diameter. 

Since  some  model  antenna  radiation  patterns  are  measured  as 
circular  cuts,  the  banked  turn  can  be  used  to  flight  measure  the  same 
patterns. 

A  disadvantage  of  the  banked  turn  technique  is  that  it  does  not 
produce  conical  antenna  patterns  which  are  the  patterns  normally  of 
concern.  Additionally,  due  to  the  high  turn  rate  in  a  banked  turn,  the 
number  of  averaged  0.5  second  data  samples  per  pattern  are  reduced. 

Where  a  skidding  turn  may  produce  up  to  five  saunples  per  degree 
azimuth,  the  banked  turn  may  produce  less  than  three  data  points  per 
degree. 

POLYGON  TURN 

The  polygon  turn  is  an  effective  technique  to  measure  an  antenna 
pattern  in  a  short  period  of  time?  yet,  avoiding  the  skidding  turn.  This 
technique  can  easily  produce  72  data  points  per  360  degrees  which  creates 
a  reasonable  representation  of  the  actual  antenna  pattern.  Figures  23 
and  24  show  a  comparison  of  an  antenna  pattern  measured  with  a  polygon 
and  a  skidding  turn. 

To  fly  a  polygon  pattern,  the  aircraft  is  required  to  be  at  a  wings 
level  attitude  when  the  heading  is  called  on  each  leg.  Once  the  heading 
is  called  a  banked  turn  should  be  executed  as  quickly  as  possible  to 
change  the  heading  approximately  five  degrees.  This  process  should  be 
repeated  until  a  360  degree  circle  has  been  completed  with  each  called 
heading  and  IRIG  "B"  time  recorded. 

The  technique  required  to  radar  vector  the  aircraft  to  the  proper 
starting  point  is  the  seune  as  for  the  skidding  turns.  The  starting  point 
should  be  such  that  the  required  diameter  of  the  turn,  in  conjunction 
with  the  drift  due  to  wind,  does  not  cause  the  aircraft  to  fly  into  an 
area  of  cancellation.  The  dleuneter  of  a  polygon  flown  by  a  C-131  on  a 
previous  test  was  12  NM. 

This  flight  technique  is  useful  in  the  event  the  aircraft  cannot 
be  turned  in  a  wlngs-level-skid.  Since  most,  if  not  all,  airplanes  can 
fly  a  skidding  turn  this  technique  will  not  normally  be  used. 
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FIGURE  23  POLYGON  fiNTENNfl  PATTERN 


FIGURE  24  SKIDDING  TURN  fiNTENNR  PATTERN 


The  main  disadvantage  is  the  large  diameter  of  the  flight  path. 
Since  only  half  the  flight  time  is  spent  in  a  turp,  the  pattern  will 
tend  to  be  large.  This  may  require  the  polygon  to  be  flown  in  segments. 


TEST  PLANNING 


ADVANCED  PLANNING 

Advanced  planning  for  a  test  of  an  aircraft  antenna  radiation 
pattern  should  determine  the  facilities  to  be  used,  the  antennas  to  be 
measured,  the  flight  technique  to  use,  and  estimates  of  flight  time  and 
computer  time.  With  the  determination  of  these  factors  the  initial  test 
planning  can  be  finalized. 

The  first  major  item  to  be  determined  in  planning  a  test  is  which 
facility  to  use.  There  are  three  facilities  which  can  be  used  to  perform 
these  tests.  The  choice  of  the  facility  v;ill  be  made  on  the  data  require¬ 
ments.  For  the  measurement  of  antenna  patterns  between  100  Mllz  and 
400  MHz,  building  275  at  the  South  Base  should  be  used.  As  planned, 
a  facility  v.’ill  be  developed  at  this  building  which  will  be  capable 
of  measuring  six  antenna  patterns  between  100  MHz  and  400  MHz  simul¬ 
taneously.  For  pattern  measurements  in  the  1.435-1.540  GHz  or  2.2- 
2.3  GHz  range  the  receivers  and  tracking  antenna  located  at  the  TM  site, 
building  5780,  can  be  used.  This  equipment  should  be  used  in  the  seune 
manner  as  the  South  Base  facility.  For  tests  wliich  involve  frequency 
ranges  unavailable  at  Edwards  AFB  or  for  tests  which  involve  a  large 
number  of  antennas,  the  facility  located  at  Griff iss  AFB,  New  York 
should  be  considered. 

The  major  dynaunic  measurement  facility  of  aircraft  antennas  is  the 
Precision  Antenna  Measurement  System  (PAMS)  of  the  Rome  Air  Development 
Center  (RADC)  at  Griffiss  AFB.  This  facility  is  located  at  the  Verona 
Test  Range,  18  miles  south  of  Griffiss  AFB.  The  PAMS  facility  is  capable 
of  measuring  12  CW,  AM,  FM,  or  pulsed  signals  simultaneously  on  any 
frequency  from  0.1  GHz  to  18  GHz.  This  system  is  also  capable  of  measuring 
integrated  power  spectral  density  between  pairs  of  frequencies  and  radar 
cross  sectional  areas.  Through  the  use  of  tracking  receiving  antennas, 
tracking  radar,  and  onboard  aircraft  attitude  recorders,  the  PAMS 
facility  is  capable  of  being  used  with  any  type  of  flight  pattern.  If 
the  frequencies  to  be  measured  are  above  1  GHz  and  are  not  between 
1.435-1.54  GHz  or  2. 2-2. 3  GHZ,  PAJIS  should  be  used  and  appendix  A  should 
be  referenced  for  an  additional  discussion  on  PAMS. 

The  flight  pattern  to  be  flovm  will  depend  upon  the  type  of  data 
required  and  the  Capabilities  of  the  aircraft.  For  most  tests,  a  conical 
antenna  pattern  v/ill  be  desired  at  a  low  depression  angle  (-5  degrees). 

In  this  case  the  skidding  turns  will  produce  the  best  antenna  patterns 
in  the  shortest  flight  time.  For  other  types  of  data  requirements, 
different  flight  techniques  may  Le  used.  If  higli  depression  angles 
are  required,  banked  turns  or  parallel  flybys  sJiould  be  used.  An  esti¬ 
mate  of  the  depression  angle  covered  at  different  altitudes  and  ranges 
can  be  found  from  figures  25  and  26. 
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ELEVATION  ANGLE 


NOTE:  Earth  Radius  ■  3A36.99  fW,  I  *  6076.1  FT 
K  -  4/3 


(ij  0001*)  aoniinv 

Figure  25  Vertical  Coverage  Chart 


ELEVATION  AHOLE  (DEO) 


The  correct  values  for  point  A  are: 


Ground  Range 
Slant  Range 
Elevation  Angle 

0 

Altitude 

Depression  Angle 


8.65  NM 
10  NM 

30.0  ^ 
0.14  ° 

30,430  FT 

30.14  " 


Given  any  two  of  the  first  five  variables,  the  others  can  be  found 
using  the  chart. 


Given  the  Depression  Angle  and  any  one  of  the  first  four  variables, 
the  others  can  be  found  using  the  chart. 

For  ranges  greater  than  40  NM  use  the  other  Vertical  Coverage  Chart. 


Figure  25  (Continued) 
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ELEVATION  ANGLE  (DEG) 


(DEG)  CAdd  to  Elevation  Angle  to  get  Depression  Angle] 


El  EVATION  ANGLE  (DEG) 


The  correct  values  for  point  A  are: 


Ground  Range 
Slant  Range 
Elevation  Angle 

0 

A1 titude 

Depression  Angle 


80  NM 
80.2  NM 
3.0  ° 

1.34  ° 

29,747  FT 

4.34  ° 


Given  any  two  of  the  f’Yst  five  variables,  the  others  can  be  found 
using  the  chart. 


Given  the  Depression  Angle  and  any  one  of  the  first  four  variables, 
the  others  can  be  found  using  the  chart. 

For  Ground  Ranges  greater  than  40  NM,  the  Slcint  Range  is  within  0.8  M4 
of  that  range. 

For  ranges  less  than  40  NM  use  the  other  Vertical  Coverage  Chart. 


Figure  26  (Continued) 
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Tc  estimate  the  flight  time  requirement,  the  number  of  flight 
patterns  flovm  must  be  determined.  Since  the  test  may  require  antenna 
patterns  with  various  aircraft  configurations,  the  minimum  number  of 
flights  required  must  be  determined.  Flying  multiple  configurations  on 
each  flight  and  measuring  multiple  frequencies  wi  11  reduce  the  total 
flight  time  required  to  complete  the  test.  Therefore,  one  flight  can 
produce  antenna  patterns  with  the  landing  gear  up  and  dovm  and  flaps 
extended  and  retracted. 

If  the  multipath  radiation  is  to  be  eliminated  with  the  Antenna 
Radiation  Prediction  Measurement  Program,  data  collected  on  radial 
flights  will  optimize  the  calculation.  This  data  can  be  used  if  the 
antenna  pattern  is  expected  to  be  essentially  constant  from  0  to  -10 
degrees  depression  angle  at  the  nose  or  tail  of  the  aircraft. 

The  total  flight  time  for  each  flight  can  be  estimated  by  summing 
the  follov/inq  times:  (1)  time  to  fly  radial,  (2)  time  to  fly  each  orbit 
(or  other  patterns),  and  (3)  time  to  vector  aircraft  to  start  of 
patterns  (five  minutes  x  number  of  patterns) . 

DETAILED  PLANNING 

After  advanced  planning  has  been  completed  and  flight  time  estimates 
provided  to  the  Project  Engineer,  the  Systems  Engineer  should  begin  the 
detailed  planning  of  the  test.  This  section  will  present  the  test 
preparation,  data  collection,  and  data  reduction  recjuirements  for  the 
detailed  planning  phase. 

Test  Preparation: 

In  preparing  for  a  test,  the  Systems  Engineer  should  first  review 
all  existing  pertinent  data  on  the  antenna  system  being  evaluated. 

Previous  model  patterns  or  previous  dyneunically  measured  patterns 
should  be  reviewed  to  locate  specific  antenna  patterns  which  may  be 
of  significant  interest.  This  would  be  any  large  holes  v^hich  may  be 
present  or  any  unexpected  holes  in  the  pattern  v^hich  may  require 
special  emphasis.  The  physical  location  of  the  antenna  on  the  aircraft 
should  be  reviev/ed  to  determine  the  potential  effects  the  aircraft 
structure  and  external  stores  may  have  on  the  antenna  pattern.  Holes  in 
patterns  may  be  caused  by  shielding  of  the  antenna  or  reradiation  by  the 
aircraft  structure.  Pilot  reports  of  unexpected  loss  of  communications 
are  one  means  to  locate  potential  holes  in  antenna  patterns.  These  reports 
can  provide  insight  into  deficiencies  of  antenna  systems. 

Planning  the  tests  will  involve  estimating  the  amount  of  support 
to  be  required  by  Technical  Services  (AFFTC/DOET) .  This  support  will 
consist  of  CDC  6500  computer  time,  Calcomp  plotter  time,  VIIF/UHF 
Antenna  Radiation  Pattern  Measurement  Facility  (building  275), 

FPS-16  radar  tracking  (SPORT)  including  radar  recording,  radar  tapes 
created  by  AFFTC/DOETDA  (Data  Analysis  Section) ,  and  digitized  signal 
tape  created  by  DOETDD  (Digital  Data  Processing  Section) . 

Estimates  of  the  time  required  for  each  of  the  above  organizations 
can  be  determined  with  the  following  equations: 
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Computer  CPU  time  in  minutes  (T^) 

■  (Flight  time  in  hours) (number  of  frequencies  recorded) 
X  (2)  -f  (number  of  plots  to  be  produced)  x  (.2) 


Calcomp  plotter  time  in  minutes  (T^) 

Tp  *  (Number  of  plots  to  be  produced)  x  (2,5) 


Time  estimates  for  tho  radar  tapes  and  digitized  signal  tape  will 
be  equal  to  the  flight  test  time.  The  trac)(lng  radar  time  estimates 
will  be  equal  to  the  flight  time  plus  15  minutes  per  flight  to  allow 
time  for  the  radar  to  initiate  traclcing. 

Time  estimates  for  the  VHF/UHF  Antenna  Radiation  Pattern  Measure¬ 
ment  Facility  will  be  equal  to  the  flight  time  plus  three  hours  per  '* 
flight.  This  will  allow  two  hours  before  the  flight  for  system  warm¬ 
up  and  pre-calibration.  One  hour  is  allov;ed  after  the  flight  for  post- 
calibration. 

Frequencies  must  be  assigned  to  the  project  from  Frequency 
Management  (1925  COMM  Sq/DCXF) .  Each  frequency  must  be  assigned 
to  the  test  project,  since  any  interference  from  other  transmitters  will 
create  erroneous  signal  strengths.  One  comiTiUnication  frequency  must  also 
be  assigned  which  cannot  be  jointly  used  with  other  projects.  Since  the 
aircraft  will  be  radar  vectored  and  the  pilot  will  be  calling  headings, 
interference  from  other  projects  will  cause  loss  of  data  and  difficulty 
in  radar  vectoring.  If  the  aircraft  does  not  have  a  second  communication 
radio  which  can  be  used  for  communicating  with  RAPCON,  control  of  the 
aircraft  can  be  maintained  by  requesting  RAPCON  to  monitor  the  test 
communication  frequency. 

After  the  test  has  been  planned  in  detail,  the  System  Engineer  should 
prepare  a  Test  Information  Sheet  (TIS),  AFFTC  Form  0-128,  for  inclusion 
in  the  final  test  plan.  This  form  should  be  completed  according  to  the 
Systems  Engineering  Handbook,  page  I-II-6.  This  should  be  accomplished 
at  least  two  weeks  prior  to  the  test. 

A  test  card  (figure  27)  should  then  be  prepared.  This  card  should 
contain  information  which  will  be  required  by  all  test  participants  and 
be  distributed  at  the  pilot's  briefing  prior  to  the  flight. 

Data  Collection t 

The  data  collected  during  a  test  will  consist  of  manually  recorded 
data  by  the  Systems  Engineer  and  magnetic  tape  recording  of  signal 
strength  and  radar  positioning.  The  data  recorded  by  the  Systems 
Engineer  will  consist  of  transmitter  power,  before  and  after  each 
flight,  pre-  and  post-calibration  information,  flight  times,  heading 
information,  and  signal  tape  identification  Information. 
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ANTENNA  PATTERN  EVALUATION 

JON  998300 
OPS  336-77 

A/C  -  F-15  No.  20113 
A/C  Call  Sign  -  Eagle  11 
Radar  Call  Sign  -  SPORT 
South  Base  Call  Sign  -  BARRIER 
Communication  Freq  -  315.9  MHz 
Takeoff  -  0800  12  Dec  1974 
Test  Altitude  10,000  MSL 
TestFreqs  -  264.6MHz 
139.8  MHz 
TEST  SEQUENCE 

1 .  Contact  SPORT  to  initiate  radar  tracking 

2.  Contact  BARRIER  for  test  freq  ID 

3.  Fly  outbound  radial 

4.  Fly  right  skidding  turn 

5 .  Repeat  4 

6.  Repeat  4 

7.  Fly  inbound  radial 

8.  Terminate  test 


Figure  27  Sample  Test  Card 
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The  power  delivered  to  each  test  antenna  at  each  test  frequency 
should  be  measured  before  and  after  each  flight.  This  measurement 
should  be  made  on  the  aircraft  with  a  calibrated  Radio  Frequency  (RF) 
power  meter.  These  measurements  will  indicate  if  the  transmitters 
are  providing  constant  power  to  the  antenna.  Additionally,  the 
measured  power  can  be  compared  with  the  antenna  patterns  to  show 
gain  or  loss  from  the  isotropy. 

The  power  measurement  should  be  in  dbw  and  the  line  losses  in  db.* 
If  power  is  measured  in  dbm  or  v/atts,  it  can  be  converted  to  dbw  with 
the  following  equations: 

dbw  *  10  log  P  -  db,.. 

W  til} 

or 

dbw  *  dbm  -  30  -  db^^^j 

where 

dbw  B  power  to  antenna  in  db  above  a  watt 

P..  *  power  measured  in  watts 
w 

dbm  =  power  measured  in  db  above  a  milliwatt 

db.r  t  ^  line  loss  in  db  between  the  power  measurement  point  and  the 

'  '  antenna 

One  of  the  most  critical  tasks  of  measuring  antenna  patterns  is 
the  calibration  of  the  receivers.  Since  the  final  data  will  depend 
upon  the  accuracy  of  the  calibration  points,  care  must  be  taken  to 
insure  the  points  are  valid.  This  will  involve  insuring  the  calibration 
equipment  is  operating  correctly,  the  correct  signal  strengths  are  used, 
and  connecting  cables  are  in  good  condition. 

In  calibrating  the  receivers  the  object  is  to  place  a  known  signal, 
from  a  calibrated  RF  generator,  at  the  input  terminals  to  the  receiver. 
This  will  produce  a  signal  which  is  recorded  on  the  FM  tape  and  can  be 
correlated  with  the  known  signal  with  time. 

To  accomplish  the  calibration,  an  RF  pov;er  meter  is  connected  to 
the  RF  signal  generator  in  place  of  the  receiver  and  the  output  of  the 
signal  generator  is  calibrated.  The  power  meter  is  removed  and  the  RF 
generator  is  connected  to  the  input  terminals  of  the  receiver.  The 
output  attenuator  of  the  signal  generator  is  reduced  in  5  db  steps  for 
each  calibration  point  with  the  signal  to  the  receiver  and  time 
recorded  for  use  by  the  data  reduction  progrzutn. 


*  Note:  Power  is  normally  measured  in  dbw  except  for  very  low  pov/er 
levels  (typically  receiver  input  terminals)  when  dbm  is  normally 
used. 
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Normally  the  calibration  points  and  the  signal  input  to  the  receiver 
should  range  from  approximately  -50  dbm  to  -100  dbm.  The  actual  range  of 
the  calibration  can  be  changed  once  the  received  signal  strength  is  known 
from  test  flights. 

For  each  flight  pattern  flown  the  aircraft  magnetic  heading  and 
TRIG  "B"  time  must  be  recorded  by  the  Systems  Engineer  as  called  by  the 
pilot.  This  data  will  be  entered  in  the  Antenna  Radiation  Pattern 
Measurement  Program  for  final  processing  of  the  data.  The  headings  may 
also  be  recorded  on  the  radar  plot  hoard  for  future  reference. 

Data  Processing; 

Processing  of  data  collected  during  the  test  will  begin  with 
the  creation  of  a  digital  signal  tape  and  radar  tape.  These  tapes 
provide  the  strength  of  each  signal  recorded  and  the  location  of  the 
aircraft  relative  to  the  receiving  antenna  for  use  by  ADAP  and  ARPMP. 

Requests  for  the  creation  of  the  radar  tapes  are  made  through 
the  Data  Analysis  Section  (AFFTC/DOETDA)  by  submitting  the  Radar/ 
Postflight  Processing  Request  (figure  28).  The  reference  point  for  the 
processing  of  the  radar  data  should  be  the  location  of  the  receiving 
antenna . 

Requests  for  digitizing  the  FM  signal  tape  are  made  through  the 
Digital  Data  Processing  Section  (AFFTC/DOETDD)  by  submitting  the 
Digitizing  Request  (figure  29).  When  digitizing  the  tape,  care  must 
be  taken  to  insure  the  data  is  associated  v/ith  the  correct  source  number 
(see  appendix  B)  .  To  insure  the  correct  data  is  used  by  ARPMP,  the 
frequency  of  the  measured  signal  should  be  included  in  the  source 
identification  (DECOM)  card  used  in  digitizing  the  signal  tape  (see 
DOETDD  personnel).  The  contents  of  this  card  are  recorded  on  the 
digitized  tape  and  printed  by  ARPMP.  To  improve  the  resolution  of  the 
signal  strength  data,  the  signals  of  interest  should  be  expanded  to 
upper  and  lover  band  edge  of  the  data  channel.  The  use  of  a  6  Hz 
lowpass  filter  may  be  used  to  reduce  the  noise  content  of  the  data. 

Once  the  radar  and  signal  digital  tapes  are  created,  the  data  may 
be  reduced  and  plotted  using  the  data  reduction  programs  ARPMP  and  ADAP. 
The  coding  of  the  proper  computer  punch  cards  to  accomplish  this  are 
found  in  appendix  B  and  C.  Through  the  use  of  these  programs  the  final 
plots  can  be  produced. 

DATA  ANALYSIS 

Analysis  of  antenna  pattern  data  will  consist  primarily  of 
analyzing  the  location  and  significance  of  holes  found  in  the  antenna 
pattern.  Since  empirically  measured  data  is  never  exact,  some  errors 
are  to  be  expected  in  the  measured  antenna  pattern. 

Figures  30  and  31  show  a  comparison  of  tv;o  measurements  of  the 
same  antenna  pattern.  These  figures  show  that,  even  though  the 
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Flgup*  29  S«»npU  Raquaat  for  Signal  T«p«  Procaaiing 
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measured  patterns  are  not  identical,  the  basic  characteristics  of 
each  pattern  are  the  s£une.  The  repeatability  of  these  patterns  are 
as  close  as  can  be  reasonably  expected. 

The  siqnificance  of  holes  in  antenna  patterns  can  be  seen  through 
an  analysis  of  the  one  way  range  equation  for  electromagnetic  radiation 
propagation. 

‘r  -  Pt  *  “"r  ♦  S  *  S  -  *•  -  109  -  20  log  F,,„,  -  7.8 

=  pov/er  to  receiver  in  dbm 
=  output  po^.’er  of  transmitter  in  dbv; 

Go  =  gain  of  receiving  antenna  in  db 
=  gain  of  transmitting  antenna  in  db 
=  gain  due  to  multipath  radiation  in  dL 

L  =  sum  of  all  line  los.ses  in  transmitting  and  receiving 
antennas  in  db 

Rj^j^  =  range  to  aircraft  in  NJl 


^IlHz  “  frequency  of  signal  in  Mllz 

By  assuming  values  for  the  minimum  discernible  signal,  gains  of  the 
antenna,  and  line  loss;  the  effect  of  holes  in  the  antenna  pattern 
can  be  seen  on  the  usable  range  of  the  transmitter.  For  a  comparison 
the  follov.’ing  values  are  assumed:  Pj^  =  -95  dbm,  P^,  =  30\.’  =  14.8  dbw, 

Gj^  =  0  db,  =  0  dbm,  L  =  8  db,  =  264.0. 

Figure  32  shows  the  effect  of  a  hole  in  the  pattern  of  an  omni¬ 
directional  antenna  on  the  maximum  usable  range  of  the  transmitter. 

For  each  10  db  loss  the  maximum  usable  range  is  reduced  by  a  factor  of 
1/3.2.  Since  most  transmissions  occur  v/ithin  100  IIM,  a  6  db  to  10  db 
loss  will  be  insignificant.  Losses  greater  than  this  can  significantly 
degrade  the  capabilities  of  the  transmitter  and  when  combined  with 
multipath  cancellation  can  result  in  complete  loss  of  communications. 

The  final  analysis  of  antenna  patterns  of  omni-antennas  must  con¬ 
sider  the  location,  magnitude,  causes  of  deviations  from  the  isotropic 
level,  and  MIL  SPEC  requirements.  For  example,  a  20  db  loss  due  to 
landing  gear  extension  may  be  2K;ceptable  but  a  10  db  loss  due  to 
external  stores  may  not  be  acceptable.  In  evaluating  the  acceptable 
losses,  the  use  of  the  aircraft  in  the  particular  configuration  should 
be  considered. 

REPORTS 

The  Systems  Engineer,  responsible  for  conducting  antenna  pattern 
tests,  will  be  required  to  prepare  progress  reports,  a  final  report,  and 
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Figure  32  Antenna  Attenuation  Effects  on  Maximum  Usable  Range 


deficiency  reports  (DR)  on  the  test  project.  These  reports  should  be 
prepared  v;ith  care  to  present  the  data  in  a  clear,  concise  manner. 

The  progress  reports  are  required  to  keep  the  System  Program  Office 
(SPO)  and  management  informed  on  the  status  of  the  test.  This  report 
should  summarize  the  results  of  the  tests  for  the  reporting  period, 
indicate  any  problems  encountered,  and  briefly  describe  what  ta'-’-" 
remain. 

The  final  report  should  contain  all  the  information  kno\/n  about 
the  antenna  system.  The  main  body  of  the  report  should  shov/  hov;  the 
test  was  accomplished;  typical  antenna  pattern  on  each  frequency, 
emtenna,  and  aircraft  configuration;  and  conclusions  and  recommendations 
on  the  ability  of  the  radiating  system  to  perform  its  designed 
function.  The  conclusions  should  specify  the  presence  of  any  large 
null  areas  and  the  probable  effect  on  the  overall  performance  of  the 
emtenna  system.  Line  losses  and  voltage  standing  wave  ratio  (VSWR)  at 
each  test  frequency  should  also  be  presented. 

The  appendix  to  the  report  should  contain  all  the  data  v;hich  lead 
to  the  conclusions  in  the  report.  This  would  include  all  the  plotted 
antenna  patterns  which  v/ere  determined  to  be  valid. 

In  preparing  tlie  final  report,  the  Systems  Engineer  should  alv;ays 
consider  the  final  user  of  the  data.  As  such,  the  information  presented 
should  contain  the  data  which  may  be  required  to  effectively  identify 
the  characteristics  of  the  radiating  systems. 
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APPENDIX  A 


PRECISION  ANTENNA  MEAvSUREHENT  SYSTEM  (PAI^S) 


The  Precision  Antenna  Measurement  System  is  a  recently  developed 
facility  v;hich  represents  the  present  state-of-the-art  in  dynamic  measure¬ 
ment  of  airborne  RE  antenna  radiating  characteristics.  This  system  is 
presently  operated  by  the  Technical  Support  Division  of  the  Rome  Air 
Development  Center  (RADC)  at  Griffiss  AFD,  Nev  York.  The  PT^IS  facility 
is  located  at  the  Verona  Annex  of  RADC  at  Verona,  Ne\’  York. 

PAMS  is  an  automated  system  xdiich  can  measure  up  to  12  CVJ,  7VM,  F.m, 
or  pulsed  signals  simultaneously  and  is  capable  of  measuring  signal 
strength,  integrated  power  spectral  density  between  pairs  of  frequencies, 
and  radar  cross  sectional  areas.  Data  may  be  recorded  with  horizontal 
or  vertical  linear  polarizations  or  left  or  right  circular  polarizations 
simultaneously.  This  gives  the  system  the  capability  of  measuring  up 
to  24  antenna  patterns  simultaneously.  Through  the  use  of  a  HP  2116B 
computer  each  frequency  is  automatically  calibrated  in  one  dbm  steps 
producing  high  accuracy  in  the  calibrated  signal. 

The  PT^IS  concept  is  illustrated  in  figure  Al.  The  system  consists 
of:  (1)  an  airborne  monitoring  recorder  system,  which  records  on  seven 

track  magnetic  tape  the  aircraft  heading,  pitch,  roll,  and  time 
for  use  during  postflight  data  reduction,  (2)  six  airborne  signal  sources 
used  to  provide  stable  10  watt  signals  betv/een  .2  GHz  and  18  GHz,  (3) 
tracking  receiving  antennas,  (4)  FPS-IG  tracking  radar,  (5)  console  to 
control  the  data  sampling  and  monitor  the  received  signals,  and  (6) 
computer  system  to  reduce  the  recorded  data. 

Through  the  use  of  the  Airborne  Monitoring  System  (AMS)  and  the 
FPS-16  radar,  the  aircraft  can  fly  any  flight  pattern.  This  enables  the 
flight  pattern  to  be  designed  to  measure  specific  data  requirements. 

The  PAMS  facility  is  the  best  system  presently  available  to  perform 
antenna  pattern  testing  above  i  GHz.  If  patterns  are  required  on  antennas 
above  this  frequency  the  PWlS  range  should  be  used. 

To  perform  a  test  at  PAMS  the  following  organization  should  be 
contacted : 


RADC/TUTV  (M.  E.  Cook) 

Griffiss  AFB  NY 

Mr.  Ed  Cook  is  the  P7VMS  engineer  and  can  be  contacted  by  autovon  at 
587-4227  to  arrange  testing  on  the  PAMS  range. 

Additional  information  on  PAMS  can  be  found  in  the  following  tv;o 
reports  from  RADC:  (1)  Precision  Antenna  Measurement  System  (PAMS), 
RADC-TR-73-233,  August  19'^3,  AD  9l558lL,  (2)  PAMS  Airborne  Instrumentation, 
RADC-TR-74-144,  July  1974. 
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APPENDIX  B 


ANTENNA  RADIATION  PATTERN  MEASUREMENT  PROGRAM  (ARPMP) 

Deck  Number  4701 


INTRODUCTION 

The  Antenna  Radiation  Pattern  Measurement  Program  (ARPMP)  is  a 
software  package  developed  to  provide  the  Systems  Engineer  with  an 
automated  system  for  measuring  aircraft  antenna  radiation  patterns.  The 
progreun  is  used  for  post  flight  reduction  of  data  recorded  on  magnetic 
tape  during  an  antenna  pattern  measurement  flight  and  calculates  the 
effective  radiated  power  (ERP)  transmitted,  aspect  angle,  depression 
angle,  horizontal  range  of  the  aircraft  to  the  receiving  antenna,  and 
time.  The  data  read  off  the  magnetic  tapes  are  in  0.1  second  Intervals. 
For  continuous  patterns  (orbits,  radials,  and  parallel  flybys)  this 
data  is  averaged  for  0.5  second  and  for  discrete  patterns  (cloverleaf 
and  polygon)  the  data  is  averaged  for  2  seconds.  The  processed  data 
is  both  printed  and  written  to  permanent  storage  for  use  by  the  Antenna 
Data  Analysis  Progr^un  (ADAP) . 

SYSTEM  OVERVIEW 

ARPMP  consists  of  a  main  program  and  14  subprograms  which  are 
each  designed  to  perform  specific  tasks  (calibration,  read  signal  tape, 
read  radar  tape ,  etc . ) .  The  main  program  and  executive  subprogram 
READD  are  used  to  control  the  operation  of  the  program. 

The  input/output  requirements  of  the  program  are  shown  in  figure 
Bl.  The  signal  tape  is  the  digitized  PM  magnetic  tape  which  contains 
the  calibration  signals  and  received  signal  strength  of  each  signal 
source  being  recorded.  The  tape  also  contains  the  IRIG  "B"  time  for 
use  in  correlating  the  data  with  radar  and  aircraft  heading  data. 

The  radar  tape  is  the  tape  created  by  the  Radar  Computer  Program  of 
AFFTC/DOETDA .  The  tape  contains  the  X,  Y,  Z  location  of  the  aircraft 
relative  to  the  receiving  antenna  and  IRIG  "B”  time.  This  tape  contains 
other  information  on  the  aircraft's  flight  path  but  only  the  location  is 
used  by  ARPMP. 

The  card  input  to  the  program  defines  all  the  parameters  necessary 
for  the  processing  of  the  tapes.  This  Includes  calibration  data,  flight 
times,  flight  pattern  type,  and  options  to  use  in  reducing  the  data.  The 
options  available  Include:  (1)  processing  radial  flight  data  without 
radar  tapes,  (2)  using  both  pre-  and  post-calibration  points  on  the 
signal  tape,  and  (3)  correcting  signal  strength  for  multipath  radiation. 
The  input  cards  select  one  of  three  basic  flight  paths  flown  by  the 
aircraft.  The  flight  paths  are:  (1)  radial  or  parallel  flyby  (constant 
heading  legs),  (2)  cloverleafs  or  polygon  (discrete  points),  and  (3) 
orbits  (continuous  changing  heading). 

The  output  data  from  ARPMP  consists  of  an  output  file  and  a  computer 
listing  (figure  B2).  The  output  file  contains  all  the  processed  data 
and  is  stored  on  either  magnetic  tape  or  disc.  The  data  file  can  be 
used  by  ADAP  to  produce  polar  plots  or  range  versus  signal  plots.  The 
output  listing  contains  all  the  input  cards  read  and  a  copy  of  all 
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the  processed  data.  This  can  be  used  to  check  the  validity  of  the 
data  prior  to  plotting. 

PREPARATION  FOR  USE 

The  program  deck  of  ARPMP  will  be  permanently  stored  on  magnetic 
tape  number  5540  in  the  tape  library  of  building  3940.  This  tape  contains 
four  files  which  are:  (1)  ADAP  LGO  (compiled  file),  (2)  ARPMP  LGO 
(compiled  file),  (3)  ADAP  source  file,  and  (4)  ARPMP  source  file.  Prior 
to  the  first  use  of  these  progreuns,  the  compiled  file  must  be  copied  to 
disc  and  stored  as  a  permanent  file  as  shown  in  figure  D3.  Once  the 
program  is  stored  as  a  permanent  file,  it  may  be  attached  and  executed. 

US/IRS  GUIDE 

The  execution  of  ARPMP  requires  a  signal  tape  (TAPElO),  radar 
tape  (TAPEll),  output  tape  (TAPE12) ,  and  a  data  deck.  The  mass 
storage  files  may  be  either  magnetic  tape  files  or  disc  files.  For 
efficiency  of  operation,  it  is  suggested  these  files  be  located  on  disc. 
This  will  reduce  the  execution  time  and  free  the  program  from  requiring 
magnetic  tape  drives. 

The  input  magnetic  tape  files  may  be  copied  to  disc  with  a  COPYBF 
macro  and  the  MERGE  option  of  ADAP.  If  more  tlian  one  radar  tapes  are 
merged,  the  file  should  be  created  with  increasing  times.  An  example 
merge  deck  is  shown  in  appendix  C. 

The  deck  developed  for  running  ARPMP  is  composed  of  two  parts. 

The  first  part  is  the  job  control  cards  to  attach,  request,  catalog, 
and  execute  the  progreun.  Examples  of  these  cards  are  shown  in  figures 
B13  and  B14. 

The  second  portion  of  this  deck  is  the  pariuneter  identification 
cards  which  supply  all  the  required  card  input  data  for  the  program. 

There  are  nine  different  types  of  data  cards  used  by  ARPMP.  The 
selection  and  coding  of  the  cards  are  described  in  the  following 
section  and  shown  in  figures  B4  through  B12.  Figures  B13  and  B14 
are  example  data  decks  used  for  ARPMP. 

Input  Cards  for  ARPMP: 

ARPMP  card  1  (figure  B4)  is  used  to  provide  a  descriptive  heading 
on  the  ARPMP  listing.  Any  comment  may  be  written  in  the  80  column 
field  to  describe  the  data  being  processed.  Normally  items  like  date, 

OPS  number,  aircraft,  and  flight  type  will  be  recorded. 

ARPMP  card  2  (figure  D5)  is  used  to  identify  the  options  to  be 
used  and  to  identify  the  run.  By  coding  a  Y  for  YES  or  an  N  for  NO 
in  the  appropriate  columns,  radar  tape,  pre-  and  post-calibration,  or 
multipath  correction  options  can  be  selected.  The  radar  tape  option  N 
is  only  selected  for  radial  flights  without  radar  tracking.  This  will 
cause  the  progreun  to  calculate  range  based  on  time  and  range  data  supplied 
by  ARPMP  card  5.  If  tvro  sets  of  calibration  points  are  to  be  averaged 
to  calibrate  the  signal  tape,  a  Y  is  coded  in  column  2.  If  only  one 
set  is  to  be  used  an  N  is  coded.  Corrections  of  the  signal  strength 
for  multipath  radiation  can  be  selected  by  coding  a  Y  in  column  3.  If 
this  option  is  selected  the  multipath  effect  is  calculated  by  the 
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Note : 


Figure  B3 


Exact  Format  of  Cards  Depends  Upon  the  Version  of 
SCOPE  Being  Used  on  the  CDC  6500  Computer 


Deck  Setup  to  Copy  ARPMP  and  ADAP  to  Permanent  File 
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technique  described  in  appendix  E.  Selecting  this  option  will  require 
ARPMP  card  4  to  be  coded. 

The  number  of  flight  patterns  to  be  processed  is  coded  in 
columns  4-10,  For  each  flight  pattern  processed  ARPMP  cards  8-9  must 
be  coded.  The  number  of  different  signals  recorded  is  recorded  in 
columns  4-10.  The  Julian  date  (0-365)  of  the  flight  is  coded  in  columns 
21-30  and  if  a  radar  tape  is  not  used  the  number  of  range  entries  is 
coded  in  columns  31-40  and  the  altitude  of  the  flight  (AGL)  is  coded 
in  columns  41-50. 

ARPMP  card  3  (figure  B6)  is  coded  to  identify  each  signal  source 
recorded  on  the  signal  tape.  The  data  recorded  on  this  card  are  the 
frequency  of  the  signal  in  Mllz,  line  losses  of  the  receiving  system  in 
db,  and  the  gain  of  the  receiving  antenna  in  db.  Data  on  tv/o  sources 
are  recorded  on  each  card  v;ith  a  maximum  of  20  sources  permitted.  The 
order  of  recording  the  sources  is  the  same  as  the  order  of  the  CAL 
records  on  the  signal  tape.  This  will  normally  be  in  the  order  of  IRIG 
channels  used  to  record  the  data  on  FM  tape. 

APJ’MP  card  4  (figure  B7)  is  used  only  if  the  multipath  correction 
option  is  selected.  For  each  signal  the  receiving  antenna  height  in 
feet,  relative  permittivity,  and  conductivity  are  entered  in  the  same 
order  as  in  ARPMP  card  3.  The  conductivity  and  antenna  height  can 
either  be  measured  or  picked  such  that  program  MULT  of  ADAP  predicts 
the  location  of  each  lobe  to  coincide  with  the  measured  lobes.  The 
relative  permittivity  which  best  eliminates  the  multipath  effect  can 
be  found  by  program  MATCH  of  AD7VP, 

ARPMP  card  5  (figure  B8)  is  only  used  if  the  No  Radar  Tape  option 
is  selected.  This  card  records  the  range  at  the  end  of  the  radial  in 
NM  and  the  time  the  aircraft  was  over  the  receiving  antenna  and 
at  the  end  of  the  radial.  To  determine  the  range  of  the  aircraft,  time 
will  be  used  for  linear  extrapolation  of  range  between  the  start  and 
end  of  the  radial. 

ARPMP  card  6  (figure  B9)  will  record  the  number  of  calibration 
points  to  be  read  on  ARPMP  card  7  (figure  BIO) .  Cards  6  and  7  will  be 
coded  for  each  set  of  calibration  points  recorded.  The  calibration 
points  should  be  recorded  in  sequence,  starting  V7lth  the  signal  identi¬ 
fied  as  source  1.  If  pre-  and  post-calibration  points  are  used,  the 
set  of  calibration  points  with  the  lowest  times  should  be  coded  first.  A 
maximum  of  50  calibration  points  per  source  is  permitted,  with  a 
meucimum  of  200  total  points.  Normally  10  calibration  points  per  source 
will  be  sufficient. 

ARPMP  card  8  (figure  Bll)  is  used  to  provide  data  on  the  flight 
type.  The  first  entry  on  this  card  Identifies  the  basic  flight  technique. 
If  the  flight  has  legs  of  constant  heading  (radial  or  parallel  flyby),  a 
1  is  coded.  If  the  pattern  flown  consists  of  discrete  points  (cloverleaf 
or  polygon) ,  a  2  will  be  coded.  For  flight  pattern  with  constant  turn 
rates  and  continuous  data  sampling,  a  3  is  coded. 

In  columns  11-20,  the  number  of  headings  recorded  on  ARPMP  card  9 
is  recorded.  A  maximum  of  90  is  permitted  per  pattern  flown.  The 
aircraft's  attitude  (roll  and  pitch)  is  recorded  in  columns  21-40. 

Right  roll  is  recorded  as  positive  and  left  roll  as  negative  as  measured 
in  degrees  from  the  horizontal.  Pitch  up  is  recorded  as  positive  and 
pitch  down  as  negative. 


12 


Card  Identifying  thd  Data  Deck 


8A10 


Parameter  Identification  Card 


trimuiiww»i!ii' 


Parameter  Identification  Card 


goaDoaBaPBiXDDDaaEaBBc-DDDDPnaapaaQOD 


_ Parameter  Identification  Card 

Height  of  Source  #1  above  Antenna  Reflecting 
Surface  (ft) 


F10.3 


Relative  Permittivity  of  Reflecting  Surface 


F10.3 


Conductivity  of  Reflecting  Surface 
mho-meter/sq  meter 


F10.3 


Antenna  Height  of  Source  #2  above  Reflecting 
Surface  (ft) 


F10.3 


I  Relative  Permittivity  of  Reflecting  Surface 


F10.3 


Conductivity  of  Reflecting  Surface 


F10.3 


Repeat  Card  for  Additional  Sources  until  data  Is 
Recorded  for  all  Sources  (Maximum  of  20) 


Figure  B7  ARPMP  Card  H 
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Column;;  41-60  are  used  to  record  the  start  and  stop  times  to  be 
used  for  processing  of  flight's  data  for  flight  type  3.  For  flight 
types  1  and  2,  the  start  and  stop  time  entries  are  not  renuired  since 
the  times  are  used  for  each  heading  entry  on  ARPMP  card  9. 

The  last  card  is  ARPItP  card  9  (figure  B12).  This  card  contains 
all  the  headings  and  heading  time  intervals  to  le  used  in  data 
processing.  The  headings  recorded  are  magnetic  headings  as  read  off 
the  aircraft's  directional  gyro.  The  headings,  start  time,  and  stop 
time  are  recorded  four  to  a  card  V7ith  a  maximum  of  90  entries.  For 
flight  types  2  and  3  the  stop  times  are  left  blank  and  if  radar 
tapes  are  not  used  the  heading  may  be  left  blank. 

Program  Description; 

ARPMP  was  developed  in  a  modular  design  to  facilitate  future  pro- 
greun  modifications.  As  such,  the  prograun  uses  14  subprograms  of  which 
READD,  CAL,  SICNAL,  RADAR,  MDLTPTII  are  the  main  subroutines.  The  main 
program  is  responsible  for  initialization  of  data,  reading  control 
data  cards,  and  initiating  the  data  processing. 

Subroutine  CAL  creates  the  arrays  required  by  VCAL  to  calibrate 
the  data  on  the  signal  tape  in  dbm.  This  subroutine  reads  data  cards 
which  give  the  time  and  signal  recorded  on  the  signal  tape.  Data 
are  read  off  the  tape  for  tv'o  seconds  after  the  time  and  the  averaoe 
number  read  is  stored  in  th.e  array  CALSE.  Each  calibration  signal  for 
each  source  is  saved  in  this  array  v/hich  is  used  by  VCAL  to  linearly 
extrapolate  signal  tape  data  to  obtain  calibrated  data. 

Suliroutine  READD  is  the  master  su^  routine  for  processing  data.  This 
subroutine  is  called  once  for  each  set  of  data  to  be  processed  and  controls 
the  calling  of  other  subroutines  to  create  the  final  data.  The  effective 
radiated  power  is  calculated  by  the  one  v/ay  range  equation: 

-  Pr  *  L  *  20  log  *  20  log  ♦  37.8  -  -  G„ 

ERP,.  -  dlgi  -  30 
db\’ 

P^  =  pov^er  received  (dbm) 

L  =  line  loss  in  receiving  antenna  (db) 

Rj^  “  slant  range  to  aircraft  (NM) 

^MHz  *  frequency  (MHz) 

G-  >  gain  of  receiving  antenna  (db) 

R 

s  gain  due  to  multipath  radiation  (db) 

The  data  array  created  by  REAJ5D  is  v/ritten  to  mass  storage  and  the  line 
printer. 

Subroutine  SIGNAL  is  called  by  CAL  and  READD  to  provide  0.5  seconds 
of  averaged  data  off  the  signal  tape  for  each  source  of  data  recorded. 

This  subroutine  returns  uncalibrated  data  on  each  source. 
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Subroutine  RADAR  is  called  by  READD  to  provide  the  slant  range, 
horizontal  range,  aspect  angle,  and  depression  angle  of  the  aircraft  to 
the  receiving  antenna.  The  subroutine  v;ill  use  the  average  X,  Y,  Z 
location  of  the  aircraft  for  0.5  seconds  interval  and,  through  the  subroutine 
PAMS,  return  the  required  data  to  READD. 

Subroutine  MULTPTII  uses  the  technique  described  in  appendix  E  to 
determine  the  effect  of  multipath  radiation  on  the  received  signal.  If  the 
multipath  option  is  selected,  this  subroutine  will  be  called  to  return 
the  db  gain  of  the  multipath  signal  at  each  data  sample. 

Two  input  tapes  are  required  for  the  calculation  of  aircraft 
antenna  patterns  by  ARPMP.  These  tapes  are  created  by  AFFTC/DOETD 
from  an  FM  analog  tape  and  an  FPS-16  radar  tape. 

The  first  input  tape  is  the  digitized  signal  tape  created  from  the 
FM  analog  signal  tape.  This  tape  is  created  in  16  bit  words  and  is 
composed  of  three  types  of  records.  Function  IWRD  of  ARPMP  is  used 
to  extract  the  16  bit  words  from  the  60  bit  v/ords  in  the  CDC  6500 
computer . 

The  first  record  recorded  on  tape  is  the  heading  record.  This 
record  is  used  to  identify  the  tape  and  has  the  format  as  shown  in 
figure  B15. 

The  second  record  written  on  tape  is  the  time  Cal  record  (figure 
B16).  This  record  contains  the  tag  number  to  identify  the  millisecond 
time  in  the  data  record.  This  tag  number  will  be  either  a  177776g  or 
7g  depending  on  the  machine  being  used  to  create  the  tape. 

For  each  data  source  recorded  on  tape,  an  additional  Cal  record 
is  written  to  identify  the  tag  number  to  be  associated  with  that  data. 

All  Cal  records  are  read  by  ARPMP  v;hich  reads  the  tag  number  to  identify 
each  data  sample  in  the  data  record. 

After  the  Cal  records  are  written,  the  data  records  (figure  B17) 
are  written  which  contain  data  sampled  at  0.1  second  intervals  off  the 
FM  tape.  All  data  points  are  written  in  pairs  of  numbers  in  this  record. 

The  first  number  contains  the  data  value  and  the  second  number  contains 
the  data  tag  number  to  identify  what  data  is  reqprded.  At  every 
millisecond  change  a  millisecond  word  is  placed  in  the  data  record.  This, 
word  is  used  to  identify  the  time  in  the  data  record. 

The  second  tape  used  by  ARPMP  is  the  radar  tape.  This  tape  is 
created  by  the  Radar  Computer  Program  as  described  by  Ryge  Data  Processing 
System,  FTC-TIM-73-2.  The  tape  is  composed  of  four  different  record, 
types  but  only  one  record  type  is  used  by  ARPMP.  The  first  13  words  on 
record  type  3  is  used  to  provide  the  aircraft  position  relative  to  the 
receiving  antenna.  The  format  of  this  record 'is  shovm  in  figure  B18. 

The  output  tape  created  by  ARPMP  consists  of  records  of  482  words 
which  contains  data  in  six  word  groups.  For  each  data  sample  the  follow¬ 
ing  words  are  placed  in  the  record i  (1)  source  number,  (2)  signal 
strength  (dbw) ,  (3)  depression  angle  (degrees),  (4)  aspect  angle  (degrees), 
(5)  horizontal  range  (tm) ,  and  (6)  time  in  total  seconds.  Word  number 
481  in  the  record  contains  the  number  of  valid  words  in  the  record 
and  word  number  482  is  set  to  1  at  the  end  of  a  data  set. 

A  listing  of  the  ARPMP  source  deck  can  be  found  in  appendix  F. 
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Record  Code 


Word  No 


1 

2 

3 

4 

5 

6 

7 

8 
9 


N-4 


1 1- 3 


N-2 


N-1 


N 

N+1 

N+2 


16  Eit  Words 


SR  =  Source  Number  of  Digitizing  Equipment  (Not  Used  by  ARPMP) 

Format  Code  =  EBCDIC  Coded  Format  Identifying  Character  (Not  Used  by  ARPMP) 
Millisecond  Time  Word  is  Placed  at  Each  Millisecond  Change  in  Record 
Words  N,  N+1,  and  N+2  Represent  Time  of  Last  Data  Word  in  Record 


Figure  B17  Data  Record  Format 


IS 


WORD 

CONTENTS 

FORMAT 

1 

101 

Integer 

2 

Total  Time  (Seconds) 

Fp 

3 

Time  (Hours) 

Integer 

4 

Time  (Minutes) 

Integer 

5 

Time  (Seconds) 

Fp 

6 

Elapsed  Time  (Seconds) 

Fp 

7 

Status  Word  For  A  Valid  Record 

Equals  'lOH  /' 

AlO 

8 

X  (Feet  North  Of  Antenna) 

Fp 

9 

Y  (Feet  East  Of  Antenna) 

Fp 

10 

Z  (Feet  Above  Antenna) 

Fp 

11 

DX/DT 

Fp 

12 

DY/DT 

Fp 

13 

DZ/DT 

Fp 

Figure  B18  Radar  Record  Format 
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ATPENDIX  C 


AOTENNA  DATA  ANALYSIS  PROGRAM  (ADAP) 
Deck  Nianber  4700 


INTRODUCTION 

The  Antenna  Data  Analysis  Program  (ADAP)  is  designed  to  assist 
the  Systems  Engineer  in  planning,  processing,  and  reducing  data 
collected  during  antenna  radiation  pattern  tests.  To  accomplish  this, 
ADAP  consists  of  a  series  of  subroutines  which  can  be  used  fort  (1) 
predicting  the  effects  of  multipath  radiation  (MULT),  (2)  determining 
the  permittivity  which  best  eliminates  the  multipath  radiation  (MATCH) , 

(3)  merging  radar  tapes  into  one  file  (MERGE),  and  (4)  plotting  processed 
data  on  the  Calcomp  plotter  (PLOT).  Through  the  use  of  ADAP,  data 
processed  by  the  Antenna  Radiation  Pattern  Measurement  Program  (ARPMP) 
can  be  plotted  for  use  in  final  reports. 

SYSTEM  OVERVIEVJ 

7U)AP  consists  of  a  main  program  and  13  subprogreuns  vrhich  are 
called  to  process  the  required  data.  Through  the  use  of  the  main 
subroutines  MPLOTS,  RDMERGE,  MULT,  and  MATCH  each  of  the  four  available 
functions  of  ADAP  can  be  executed  (figure  Cl). 

Mass  storage  input/output  requirements  for  ADAP  will  depend 
upon  the  function  being  used.  The  mass  storage  requirements  for  each 
function  ares  (1)  MULT  -  requires  no  mass  storage,  (2)  MATCH  -  requires 
one  input  (TAPEll)  mass  storage  file,  (3)  PLOT  -  requires  one  input 
(TAPE12)  and  one  output  (TAPE13)  mass  storage  file,  and  (4)  MERGE  - 
requires  one  to  three  input  (TAPEll,  TAPE12,  TAPE13)  and  one  output 
(TAPElO)  mass  storage  file. 

The  card  input  to  ADAP  consists  of  11  different  cards.  The  cards 
to  be  used  for  each  run  will  depend  upon  the  options  of  ADAP  selected. 

A  decision  chart  for  the  selection  of  the  proper  cards  is  shovTn  in 
figure  C2  and  a  description  of  the  use  of  each  card  is  presented  in  the 
Users  Guide  section  of  this  appendix.  Example  deck  structures,  output 
listing,  and  plots  are  included  to  assist  in  the  use  of  this  data 
reduction  program. 

PREPARATION  FOR  USE 

The  program  deck  of  ADAP  v/ill  be  stored  with  ARPMP  on  magnetic 
tape  number  5540  in  the  tape  library  of  building  3940.  This  tape  will 
contain  the  follovrlng  four  files:  (1)  ADAP  LGO  (compiled  file),  (2) 

ARPMP  LGO  (compiled  file),  (3)  ADAP  source  file,  and  (4)  ARPMP  source 
file.  Prior  to  the  use  of  the  program  the  ADAP  LGO  file  must  be  copied 
to  disc  as  shown  in  figure  C3.  With  the  LGO  file  cataloged  on  disc 
it  may  be  attached  and  executed. 

USERS  GUIDE 

The  execution  of  each  function  of  ADAP  is  controlled  by  ADAP  card  1 
(figure  C4).  This  card  will  select  the  option  to  be  used  by  coding  the 
proper  name  starting  in  column  1.  When  using  the  functions  of  ADAP, 
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ADAP 


Figure  CZ  ADAP  Decision  Tree 


iwj  ij  i^itcrr 


j— m>HjiiUg  mTjfii" 


Note:  Exact  Format  of  Cards  Depends  Upon  the  Version  of 
SCOPE  Being  Used  on  the  CDC  6500  Computer 


Figure  C3  Deck  Setup  to  Copy  ARPMPand  ADAP  to  Permanent  File 


uijuNUjg^^ 
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Parameter  Identification  Card 


Pointer  to  One  of  Four  Major  Subroutines. Code 
Either  MERGE,  PLOT,  MATCH,  or  MULT  Starting  in 
Column  1 . 


Rest  of  the  Card  is  Blank 


g 


Figure  C4  ADAP  Card  Hi 

104 


care  must  be  taken  to  insure  the  proper  name  is  given  to  the  mass 
storage  files.  Figure  C5  lists  the  files  required  for  each  function 
of  ADAP. 

MERGE  Function; 

Execution  of  the  MERGE  function  of  ADAP,  which  combines  radar  tapes 
used  by  ARPMP,  reduces  the  data  storage  requirement  of  radar  data 
by  95  percent  and  reduces  the  computer  execution  time  of  ARPMP  by 
about  50  percent.  MERGE  can  be  used  to  either  store  the  radar  data  on 
disc  or  magnetic  tape.  Wlienever  possible,  it  is  suggested  that  MERGE 
be  used  to  store  the  radar  data  on  disc,  even  for  flights  having  onlv 
one  radar  tape,  due  to  the  savings  of  computer  time  and  the  elimination 
of  the  requirement  for  magnetic  tape  drives  by  ARPMP. 

The  only  data  card  required,  other  than  ADAP  card  1,  is  ADAP  card 
2  (figure  C6)  with  the  number  of  tapes  to  be  merged  coded  in  column  10. 
When  MERGE  is  used  the  tapes  should  be  merged  in  order  of  increasing 
times  with  the  tape  having  the  lowest  time  on  input  file  TAPEll.  The 
merged  data  v;ill  l)e  written  on  output  file  TAPEIO.  An  example  deck  for 
MERGE  is  shown  in  figure  C7. 

PLOT  Function; 

The  PLOT  function  of  ADAP  can  produce  either  polar  plots  of  the 
aircraft  antenna  pattern  at  selected  depression  angles  or  cartesian 
plots  of  range  (Nf<)  versus  signal  (effective  radiated  power  dbv')  of 
radial  flight  data.  The  input  data  is  read  from  file  name  TAPEi2  which 
was  created  by  ARPMP.  The  output  file  is  a  Calcomp  plot  tape  v/hicli  is 
written  to  file  name  TAPE13,  The  printer  listing  created  for  the 
range  plots  contains  all  the  data  points  used  for  the  plot.  The 
polar  plot  listing  contains  all  the  0.5  degrees  of  aspect  angles  which 
contain  data,  the  number  of  data  samples  in  the  aspect  angle  interval, 
average  signal,  and  standard  deviation  of  the  data  in  the  interval. 
Connecting  lines  will  be  dravm  through  all  data  points  on  the  polar 
plots  which  are  within  five  degrees  in  aspect  angle.  Example  decks, 
listings,  and  plots  are  shovm  in  figures  C14  through  Cl 7. 

When  calling  the  PLOT  function  the  first  card  coded  is  ADAP  card 
1  with  PLOT  and  the  second  card  coded  is  ADAP  card  3  (figure  C8). 

This  card  indicates  the  type  of  plots  (P  for  polar  and  S  for  range)  and 
the  time  interval  to  read  data  off  the  input  file. 

The  third  card  coded  is  ADAP  card  4  (figure  C9).  This  card 
identifies  the  data  to  be  processed  and  defines  the  scaling  to  be 
used  on  the  plots.  The  first  entry  on  the  card  is  the  source  number  of 
the  data  to  be  plotted.  The  source  numbers  are  set  by  ARPMP  and  can 
be  found  in  tlie  ARPMP  output  listing. 

The  next  items  required  on  this  card  are  the  minimum  and  maximum 
values  to  be  written  on  the  signal  axis  scale.  For  polar  plots,  only 
the  maximum  value  is  required.  For  range  plots  the  values  should  be 
picked  such  that  [(maximum  signal)  -  (minimum  signal)]  t  y  axis  length  = 
N 

1,  2,  or  4  X  10  ,  where  N  is  any  integer  number.  This  will  Insure  the 
Y  axis  is  scaled  to  an  easily  read  value.  The  range  Interval  is  used 
only  for  range  plots  and  should  be  set  in  the  stune  manner  as  the 
signal  interval. 
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FUNCTION 


FILE  NAME 


DESCRIPTION 


j - 

j  MULT 

1 

No  mass  storage  files  required. 

MATCH 

TAPEl  1 

Input  file  created  by  ARPMP  which 
contains  radial  flight  data. 

1  PLOT 

) 

( 

1 

} 

TAPE12 

Input  file  created  by  ARPMP  which 
contains  data  to  be  plotted. 

1 

! 

i 

TAPEl  3 

Output  file  plot  tape  to  be  used  by 
the  Calcomp  plotter. 

1 

j  MERGE 

TAPEll 

Input  file  radar  tape  to  be  merged 

I  or  condensed.  Contains  the  lowest 

I  times  for  merging. 


TAPE12 


Second  input  tape  for  merging  with 
TAPEll. 


TAPE13 


TAPEIO 


Third  input  tape  for  merging  with 
TAPEll  and  TAPE12. 

Output  file  of  merged  or  condensed 
radar  tapes. 


Figure  C5  Mass  Storage  Files  of  ADAP 
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Parameter  Identification  Card 


Number  of  Tapes  to  be  Merged 


Figure  C6  ADAP  Card 
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Exact  FoiTTaat  of  Cards  Depends  Upon  the  Version  of 
SCOPE  Being  Used  on  the  CDC  6500  Computer 

Figure  C7  Example  Deck  for  MERGE 


The  last  t\7o  items  on  this  card  are  the  vrild  point  check  parameters. 
The  minimum  and  maximum  signals  to  be  accepted  and  coded.  Any  signal 
which  falls  outside  this  interval  is?  ignored  by  the  prograun  This  is 
used  to  insure  data  is  not  plotted  v.’hich  would  be  off  the  plot  scale. 

The  fourth  card  required  is  a  heading  card  (ADAP  card  5  figure  CIO) 
which  provides  70  columns  of  heading  information  to  be  printed  on  the 
listing  and  witten  on  the  plots.  Information  to  identify  the  plots 
should  be  coded  on  this  card. 

The  fifth  card  to  be  coded  depends  upon  v.'hich  plot  routine  is  to 
be  used.  If  range  versus  signal  plots  are  to  be  produced,  ADAP  card 
7  (figure  Cl 2)  v/ill  be  coded  v’ith  the  length  of  the  X  axis  and  Y  axis. 

For  polar  plots  ADAP  card  6  (figure  Cll)  is  coded  v;hich  contains  tlie 
depression  angle  of  the  data  to  be  plotted,  scale  interval,  depression 
angle  interval,  type  of  plane  to  dravr,  and  polar  grid  drawing  code. 

The  polar  plot  v/ill  consist  of  all  data  vrhich  lies  v^ithin  the 
depression  angle  plus  or  minus  the  depression  angle  interval.  If  a 
999.0  is  coded,  for  the  depression  angle,  all  data  v.'ill  be  plotted 
regardless  of  the  actual  depression  angle. 

The  plane  type  should  be  set  to  a  "1"  to  draw  the  airplane  in 
the  center  of  the  plot.  The  paper  selection  should  he  set  to  a  "1"  if 
a  polar  grid  is  to  be  drav.m  on  plain  vrhite  plot  paper  or  set  to  a  "0" 
if  polar  grid  paper  is  used.  If  the  polar  grid  is  to  be  drav;n,  the 
plotter  request  card  should  specify  pen  2  as  v/et  black  and  pen  1 
as  black  ballpoint.  Instructions  should  be  given  to  center  pen  1  on 
the  perforations  and  to  use  200  plain  v;hite  paper.  If  polar  paper  is 
used,  instructions  should  h^e  given  to  center  pen  1  in  center  of  the 
polar  paper  (figure  C13),  Range  plots  should  specify  v’et  pen  1  black 
and  paper  201  or  202, 

ADAP  cards  3  through  7  are  repeated  for  each  plot  to  be  produced. 
Program  PLOT  will  continue  to  produce  plots  until  all  data  cards  are 
read .  Due  to  different  plotter  paper  requirements  polar  plots  and 
range  versus  signal  plots  cannot  be  created  on  the  same  run. 

HULT  Function; 

The  multipath  radiation  prediction  option  (nULT)  of  ADAP  can  be 
used  for  predicting  the  multipath  effects  to  be  expected  at  -t  given 
altitude  and  range.  This  program  uses  the  technique  presented  in 
appendix  E  to  calculate  the  multipath  effects.  Through  the  use  of 
predicted  multipath  effects,  the  range  and  altitude  can  be  determined 
such  that  the  aircraft's  flight  path  v^ill  avoid  areas  of  signal 
cancellation.  An  example  deck  and  output  listing  are  shov7n  in  figures 
C18  and  C19. 

The  first  card  coded  should  be  ADAP  card  1  v/ith  the  MULT  function 
selected.  The  second  card  coded  will  be  ADAP  card  8  (figure  C20)  v^ith 
each  altitude  to  be  processed,  A  maximum  of  eight  altitudes  per  run 
are  permitted.  For  each  altitude  listed,  a  printer  plot  will  be  produced 
for  each  ADAP  card  9. 
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The  card  that  controls  the  number  of  plots  to  be  produced  is  ADAP 
card  9  (figure  C21).  This  card  may  be  repeated  as  many  times  as  desired. 
The  entries  required  are  range  interval,  receiving  antenna  height, 
frequency,  polarization,  permittivity,  and  conductivity.  For  ease  in 
reading  the  printer  plots,  the  range  interval  should  be  picked  such  that 

N 

(maximum  range  -  minimum  range)  t  10  =  1,  2,  or  5  x  10  v/here  N  is  any 
integer  number.  Freauency  is  recorded  in  MHz  and  polarization  is 
equal  to  zero  for  horizontal  or  one  for  vertical.  The  relative  per¬ 
mittivity  and  conductivity  of  the  reflecting  surface  are  entered  as  the 
last  two  items  on  the  card. 

HATCH  Function; 

The  multipath  matching  option  (HATCH)  of  ADAP  can  be  used  to 
determine  the  permittivity  which  produces  the  best  elimination  of  tlic 
multipath  effect.  This  progreun  uses  the  output  file  created  by  ARPMP 
on  a  radial  flight.  The  radial  is  corrected  for  multipath  radiation 
by  MATCH  using  the  technique  presented  in  appendix  E.  Due  to  inaccuracies 
in  the  calculation  with  high  angle  reflections,  the  program  will  not 
use  data  which  has  a  reflection  angle  greater  than  12  degrees.  Through 
an  iterative  process  liATCH  will  determine  the  permittivity  v/hich  produces 
the  multipath  corrected  radial  with  the  least  deviation  from  the  data's 
mean.  When  tliis  value  is  found  the  measured  radial,  multipath 
corrected  radial,  and  predicted  multipath  effect  are  plotted  on  the  line 
printer.  An  example  deck  and  listing  are  shown  in  figures  C24  and  C25. 

To  run  the  HATCH  option  ADAP  card  1  is  coded  with  MATCH  and  ADAP 
card  10  (figure  C22)  is  coded  v;ith  an  80  colximn  comment  heading  to 
print  a  heading  on  the  output  listing.  ADAP  card  11  (figure  C23)  should 
then  be  coded  with  the  time  interval,  antenna  height,  aircraft  altitude, 
frequency,  and  data  tag  number.  This  information  v/ill  identify  the 
data  to  be  processed  and  parameters  required  to  calculate  the  multipath 
effect. 
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Figure  C8  ADAP  Card  #3 
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Parameter  Identification  Card 


HonnonannaronnanaoDg: 


Parameter  Identification  Card 


Parameter  Identicatlon  Card 


X  Axis  Length  In  Inches 


FIO.J 


Y  Axis  Length  In  Inches 


F10.3 


Rest  of  Card  Left  Blank 


N 


Figure  C12  ADAP  Card  #7 
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FIGURE  C17  EXflflPLE  RANGE  PLOT 
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Figur*  Cl?  (Contiiraad)’ 


Parameter  Identification  Card 


Parameter  Identification  Card 


Parameter  Identification  Card 


Parameter  Identification  Card 


Note:  Exact  Format  of  Cards  Depends  Upon  the  Version  of 
SCOPE  Being  Used  on  the  CDC  6500  Computer 


Figvjre  C24  Example  MATCH  Deck 
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Figure  CZ5  Example  MATCH  Ll.tlnf 


Of  Sl^at-L  CilTfl  luScJ-  £XU.  CiM>/S 


Figure  C25  (Continued) 
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APPENDIX  D 


ANTENNA  PATTERN  MEASUREMENT  TEST 


The  antenna  radiation  pattern  measurement  technique  dt,vex  ^ped  for 
this  report  was  the  result  of  tests  performed  in  the  Fall  of  1^74  at 
Griffiss  New  York  and  Edwards  AFB,  California.  These  tests  were 

designed  to  determine  the  ability  of  different  flight  techniques  to 
accurately  measure  antenna  pattern  and  to  determine  the  anterna  pattern 
measurement  system  to  be  established  at  the  AFFTC. 

This  test  was  conducted  in  tv;o  phases  using  the  Space  and  Missile 
Test  Center's  (SAMTEC)  Radar  Calibration  C-131  (#803).  The  first 
phase  was  at  Griffiss  AFB  using  the  Precision  Antenna  Measurement  System 
(PAMS)  and  the  second  phase  was  at  Edwards  AFB  using  equipment  at  the 
Barrier  Instrumentation  Facility,  building  275. 

GRIFFISS  AFB  PHASE 

The  Griffiss  AFB  phase  v;as  conducted  as  a  joint  project  betv;een 
SAMTEC  and  AFFTC.  The  aircraft  was  scheduled  by  SAMTEC  for  a  complete 
test  of  its  antennas  on  the  PAMS  range.  Therefore,  to  reduce  cost 
AFFTC  and  SAMTEC  combined  both  test  requirements  into  one  test  program 
which  was  conducted  in  October  1974 , 

During  the  tests  at  PAMS,  fifteen  antennas  were  evaluated  at 
eight  different  frequencies.  This  resulted  in  over  800  antenna  pattern 
plots  being  produced.  To  excite  the  emtennas,  the  aircraft  was  modi¬ 
fied  at  Edwards  AFB  with  racks  and  wiring  for  the  installation  of 
tremsmitters  and  the  PAMS  Airborne  Monitoring  System  (AMS)  for  the 
recording  of  the  aircraft  altitude  during  the  PAMS  tests.  The  final 
installation  of  equipment  was  made  at  Griffiss  AFB.  The  exciters  used 
for  t^e  antennas  being  evaluated  for  SAMTEC  were  the  airborne  RF  signal 
sources  provided  by  PAMS.  The  exciter  used  for  the  antennas  being 
evaluated  for  AFFTC  v/ere  standard  communication  transmitters  which  were 
rated  for  continuous  transmission.  The  136.8  MHz  CVfl  signal  was  provided 
by  a  Wilcox  807A  transmitter  and  the  225  MHz  CW  signal  by  an  AN/PRC-6G 
portable  transmitter. 

The  fliglit  patterns  flovm  at  Griffiss  consisted  of:  (1)  radials, 

(2)  parallel  flybys,  (3)  cloverleaf,  (4)  orbits,  and  (5)  polygon. 

Due  to  difficulties  with  the  fiHS,  data  on  the  orbit  and  polygon  flights 
were  lost.  The  data  on  139.8  MHz  and  225  MHz  were  collected  using  a 
circular  polarized  receiving  antenna  even  though  the  transmitted  signals 
were  vertically  polarized.  This  enabled  data  to  be  collected  on  nearly 
all  the  flights  at  Griffiss  AFB,  since  the  ground  antenna  was  being 
used  as  circular  polarization  on  all  flights  for  SAMTEC.  This  caused 
all  the  plots  of  the  139.8  MHz  and  225  MHz  data  to  be  3.01  db  lov;. 

A  total  of  six  parallel  flyby  patterns  were  flov/n,  one  radial 
flight,  one  orbit  flight  and  one  cloverleaf  flight  v/ith  data  collected 
on  139.8  MHz  and  225  MHz.  The  flight  patterns  flown  on  the  parallel 
flyby,  cloverleaf,  and  radial  patterns  are  shown  in  this  report  in 
figures  Dl  tlirough  D3. 
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LocAtioa  of  first  225  HUz  Ijooa 


Figure  D2  Cloverleaf  Pattern  Griffis 


The  flights  flown  on  the  PAIIS  range  for  the  S/UITEC  and  AFFTC  tests 
required  50  hours  of  flying  time  which  v;ere  flov/n,  for  the  most  part, 
from  0030  to  0830  hours.  The  flights  at  night  were  required  to  avoid 
interference  of  tlie  UHF  signal  from  local  TV  stations. 

EDWARLS  AFB  PHASE 

The  second  phase  of  the  tests  conducted  at  Edv;ards  AFB  was  designed 
to  determine  the  optimal  technique  to  be  used  at  the  AFFTC  and  the 
equipment  required  for  an  antenna  measurement  system  at  Edwards  AFB.  As 
such,  the  Scime  type  of  patterns  were  flov.'n  at  Edwards  as  were  flown  on 
the  PAf^S  range  except  modified  parallel  flybys  were  flown  instead  of 
parallel  flybys  (figures  D4  and  D5). 

The  tests  v;ere  flov/n  in  December  1974  using  the  C-131  #803  v;ith  the 
Wilcox  OOIA  VHF  transmitter  on  139.8  MHz  and  the  AM/PRC-66  UHF  trans¬ 
mitter  on  264.6  MHz.  The  UHF  signal  was  not  the  same  frequency  as  the 
Griffiss  AFB  tests  because  the  225  MHz  frequency  was  not  among  the 
AFFTC  assigned  frequencies.  Therefore,  slight  differences  are  to  be 
expected  between  the  UHF  antenna  patterns  of  the  PAME  and  AFFTC  tests. 

The  equipment  available  at  the  Barrier  Test  Facility  located  at 
Soutli  Base  v;as  used  to  record  the  signal  strength  data  during  the  tests. 

The  instrumentation  system  provided  the  basic  equipment  required  to 
record  the  data  and  the  IRIG  "B"  time  required  for  time  correlation 
with  radar  position  data.  The  equipment  indicated  in  the  lov/er  right 
section  of  figure  D6  was  available  at  the  Barrier  Facility  and  was 
used  for  the  tests.  The  facility  provided  a  14  track  Ampex  FR  1600 
instrumentation  tape  recorder  which  was  used  to  record  the  data  and 
IRIG  "B"  time  code  as  received  from  the  base  timing  system.  The  VHF 
time  code  receiver  and  demultiplexer  provided  time  data  whicli  was 
recorded  on  tracks  11  and  12  of  the  FR  1600  tape  recorder.  The  signal 
strength  data  v;as  recorded  on  track  2  of  the  tape  recorder  using  an 
IRIG  channel  10  and  an  IRIG  channel  13  voltage  controlled  oscillator 
(VCO) .  The  automatic  gain  control  (AGC)  voltage  of  the  Stoddard  Model 
NM-30A  VHF  receiver  was  filtered  with  a  lowpass  filter  to  remove 
modulation  components  from  the  AGC  voltage.  This  filter  was  not  necessary 
with  the  TR-104  telemetry  receiver  used  for  the  UHF  signal  receiver. 

For  the  purpose  of  measuring  received  signal  strength,  the  Barrier 
Instrumentation  System  was  supplemented  v/ith  additional  equipment.  An 
omnidirectional  VHF  and  UHF  antenna  system  v/as  installed  at  an  elevation 
of  27  feet  above  ground  level.  The  antennas  v/ere  connected  to  the 
receivers  by  approximately  45  feet  of  RF-213  coai'.ial  cable.  The  signal 
loss  was  measured  and  recorded  for  each  cable  at  the  test  frequencies 
of  139.8  MHz  and  264.6  MHz. 

The  UHF  signal  was  received  by  an  Astro  Communication  Laboratory 
Type  TR-104  Telemetry  Receiver.  The  ACC  voltage  ranged  from  0  volts  at 
-110  dbm  signal  to  approximately  -5  volts  at  -40  dbm  signal.  This 
voltage  was  used  as  the  input  to  the  model  VC-50  VCO  operating  at  IRIG 
channel  13.  The  VCO  was  adjusted  so  that  -110  dbm  signal  was  approximately 
90  percent  of  upper  band  edge  and  -40  dbm  signal  was  approximately  10 
percent  of  upper  Land  edge. 
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Figure  IM  Modified  Parallel  Flyby  Flight  Pattern  Edwards  AFB  Test 


A  Stoddard  Radio  Field  Intensity  Meter  was  used  to  receive  the  VHP 
signals.  The  AGC  voltage  from  this  receiver  was  used  as  an  input  to  a 
VC-50  VCO  operating  on  IRIG  channel  10.  It  was  similarly  adjusted  so 
that  -110  dbm  signal  was  near  upper  band  edge  and  a  -40  dbm  signal 
drove  the  VCO  toward  lower  band  edge.  The  two  VCO  outputs  were  mixed 
and  both  were  recorded  on  track  2  of  the  tape  recorder. 

The  signal  loss  from  the  antenna  to  receiver  input  was  measured 
for  each  coaxial  line  at  the  test  frequencies.  This  was  accomplished 
by  first  measuring  the  output  of  the  signal  generator  v/ith  the 
power  meter  directly  at  its  output  terminal.  The  feed  line  was  then 
inserted  between  the  signal  generator  and  the  power  meter  and  the 
difference  in  pov?er  was  noted.  The  signal  loss  for  the  coaxial  lines 
used  between  the  signal  generator  and  the  receivers  were  similarly 
measured  at  the  test  frequencies. 

Prior  to  the  start  and  after  the  completion  of  each  flight,  a 
series  of  calibration  points  were  recorded  on  the  FM  tape.  The  signal 
generator  provided  the  calibration  signals  to  the  receiver  terminals  in 
5  dbm  increments.  This  produced  signal  calibrations  from  -40  dbm  to 
-90  dbm  for  the  SHF  receiver  ^uld  -60  dbm  to  -110  dbm  for  the  UHF 
receiver.  The  IP.IG  "B"  time  and  signal  strengtli  were  recorded  for 
each  calibration  point  for  use  by  the  Antenna  Radiation  Pattern 
Measurement  Progreun  to  calibrate  the  recorded  signal  strength  in  dbm. 

DATA  AIJALYEIS 


Data  collected  during  both  phases  of  this  test  progreun  were 
used  for  analysis  of  the  multipath  effects  and  flight  teclmiques. 

Through  the  use  of  radial  flight  data,  the  ability  to  elinrinate  multi- 
path  radiation  was  evaluated  and  through  the  use  of  polar  plots  of  the 
antenna  patterns,  flight  techniques  were  reviev/ed  for  accuracy  and 
repeatability  of  data. 

Radial  Flight  Data; 

The  radial  flights  were  designed  to  map  the  multipath  radiation 
effects  at  both  test  facilities  and  to  demonstrate  the  ability  to 
eliminate  the  multipath  effect  through  computer  prediction.  The  Edwards 
AFD  radials  were  plotted  as  signal  (dbw)  versus  range  and  the  Griffiss 
AFB  radials  were  platted  as  signal  and  range  versus  time.  Since  the 
receiving  antenna  at  PAMS  was  receiving  the  vertically  polarized  trans¬ 
mitted  signal  as  circular  polarized,  both  left  and  right  circular  polar¬ 
izations  are  plotted.  Both  polarizations  should  produce  similar  data, 
but  transmission  path  effects  cause  differences  to  exist. 

The  Griffiss  AFB  radials  show  that  a  multipath  lobing  structure 
exists  for  the  139.8  MHz  2uid  225  MHz  data  (figures  D7  through  DIO).  The 
lobes  are  large  enough  to  produce  significant  effects  on  the  polar 
plots  and  are  especially  significant  between  -10  and  -90  degree  depression 
£uigles,  where  the  aircraft  was  within  10  NM  of  the  receiving  antenna. 

The  radials  also  show  that  multipath  radiation  is  being  received  by 
the  six  foot  dish  antenna  on  14  35  MHz  (figures  Dll  and  D12).  At  27  NM 
and  35  NM  large  cancellations  occur  but  data  at  less  than  25  NM  does  not 
show  any  multipath  effect.  This  is  due  to  the  beeun  width  of  the  six 
foot  antenna  being  sufficiently  narrow  to  eliminate  the  multipath  radiatior 
at  the  higher  antenna  elevation  angles. 


The  Edwards  radial  data  (figures  D13  through  D16)  show  the  multi- 
path  lobing  structure  which  existed  during  the  Edv;ards  AFB  tests.  The 
lobing  structure  was  close  to  that  which  was  expected.  Since  all 
orbits  and  modified  parallel  flybys  v;ere  flown  from  20  TO*  to  30  NM, 
only  the  264.6  Mllz  data  would  be  distorted  due  to  the  large  null  which 
exists  between  15  NI'i  and  24  NM. 

The  ability  to  eliminate  the  multipath  signals  by  computer  predictions 
is  shovm  in  fiaures  D17  through  D20.  The  areas  of  concern  are  from 
20  Nf!  to  30  Nfl  where  the  orbit  and  modified  flvby  flights  were  flov'n. 

In  this  range  all  the  radials  are  reduced  to  within  +1  dl  .  Large  errors 
exist  up  to  16  NM  in  the  264.6  MHz  radial,  corrected  for  multipath 
radiation  (figures  DIO  and  D20),  due  to  approximations  used  in 

multipath  predictions  (appendix  E) .  (Rl)^ 

Additional  radials  were  flovTi  at  Edwards  AFB  to  determine  if  the 
TM  site,  building  57fl0,  could  be  used  to  measure  aircraft  antenna 
radiation  patterns  in  the  VIIF/UHF  range.  The  data  collected  shov^ed  the 
effects  of  multipath  radiation,  v»ith  an  omnidirectional  receiving 
antenna,  are  too  severe  and  unpredictable  to  use  this  site  for  VHF/ 

UHF  aintenna  pattern  measurement.  This  is  probably  due  to  the  height 
of  the  receiving  antenna  above  the  reflecting  surface  producing 
a  large  number  of  lobes  and  the  changing  terrain  of  the  reflecting 
surface  modifying  the  lobing  structure. 

Polar  Plot  Data  (Griffiss  AFB  Test) ; 

The  polar  plot  data  from  the  Griffiss  AFB  test  (figures  D21  tlirough 
D27)  was  collected  v;ith  the  parallel  flyby  flight  pattern  (figure  Dl) 
and  a  cloverleaf  (figure  D2).  Since  multipath  radiation  vrar  being 
received  during  these  tests,  the  effect  the  reflected  signal  has  on  the 
recorded  data  must  be  considered  when  analyzing  the  antenna  patterns. 

The  data,  used  in  the  0  to  -7  degree  depression  angle  polar  plots 
of  the  Griffiss  AFB  test,  consists  of  all  data  collected  during  the 
parallel  flvby  flight  pattern  except  the  data  when  the  aircraft  was 
within  10  Nf’  of  the  PAMS  site. 

The  data  used  in  creating  the  0  to  -7  degree  depression  angle  pxslar 
plots,  from  the  parallel  flvbv  flight  oattern,  was  collected  from  10  NM 
to  35  NM  from  the  receiving  antenna.  Data  from  all  legs  of  the  flight 
pattern  are  the  averaged  for  each  0,5  degree  in  azimuth  angle.  Since 
most  azimuth  angles  are  covered  in  more  than  one  leg,  data  are  being 
averaged  from  different  legs  of  the  flight  path. 

The  resulting  effect  of  averaging  data  on  different  legs,  in  the 
presences  of  multipath  radiation,  will  depend  on  the  lobing  structure 
of  the  signal.  For  the  139,8  MHz  data,  one  lobe  and  one  null  was 
included  for  the  data  collected  from  20  to  160  and  200  to  340  degrees 


Numerals  preceded  by  an  R  within  parentheses  at  the  end  of  a  paragraph 
correspond  to  the  recommendation  numbers  tabulated  in  the  Conclusions 
and  Recommendations  section  of  this  report. 


azimuth  angle  around  the  aircraft.  This  has  resulted  In  reducing  the 
multipath  effect  In  these  azimuth  angles. 

For  the  azimuth  angles  around  the  nose  and  tall  of  the  aircraft 
the  data  being  used  come  from  the  0,  1,  and  3  NM  legs.  Since  these 
are  essentially  radial  legs,  the  lobing  effect  can  be  seen  in  the  polar 
plots  (figures  D21,  022,  D25,  and  026). 

As  the  frequency  is  Increased  more  lobes  will  be  present  in  the 
data.  This  has  the  effect  of  increasing  the  jagged  appearance  of  the 
data  over  the  nose  and  tail  of  the  aircraft  and  possibly  creating 
fictitious  holes  in  the  antenna  pattern  between  30  to  160  and  210  and 
330  degrees. 

Polar  Plot  Pata  (Edwards  AFB  Test) ; 

The  polar  plot  data  from  the  Edwards  AFB  test  was  collected  using 
orbits,  cloverleaf,  and  a  modified  parallel  flyby  flight  pattern 
(figure  04).  All  data  was  collected  from  20  NM  to  30  NM  from  the 
receiving  antenna  and  in  the  presence  of  multipath  radiation. 

The  effect  of  multipath  radiation  on  the  antenna  patterns  can 
easily  be  seen  in  figures  046  through  049,  Since  the  orbit  turns  were 
conducted  from  20  NM  to  30  NJl,  the  aircraft  entered  the  first  null 
of  the  264.6  MHz  signal.  This  resulted  in  the  apparent  depression  of 
the  antenna  pattern  on  one  side  of  the  aircraft.  For  the  banlced 
turn  the  depression  v’as  not  so  pronounced  due  to  the  small  turn  radius 
keeping  the  aircraft  from  entering  the  deep  portion  of  the  null. 

Through  the  use  of  the  multipath  correction  option  of  the  Antenna 
Radiation  Pattern  Measurement  Program  (ARPMP  appendix  B) ,  the  multipath 
effect  was  effectively  eliminated  (figures  058  through  064) . 

A  comparison  of  skidding  turns  and  polygon  turns  of  the  139.8  MHz 
antenna  patterns,  corrected  for  multipath  radiation  (figures  037  through 
045),  shov/s  that  the  general  shape  of  the  patterns  are  the  same.  Since 
there  are  numerous  factors  which  can  affect  the  received  signal,  an 
exact  comparison  of  the  antenna  patterns  cannot  be  expected.  The  primary 
concern  for  comparing  the  antenna  patterns  is  the  location  and  magnitude 
of  losses  in  effective  radiation  power  (ERP) .  Comparing  the  polar  plots 
for  these  items  indicates  the  skidding  turns,  polygon,  and  modified 
parallel  flyby  produce  essentially  equivalent  antenna  patterns.  Due  to 
savings  in  flight  time  tad  high  quality  of  data,  the  skidding  turn  flight 
technique  was  judged  the  best  technique,  (R2) 

A  comparison  of  data  between  the  two  phases  of  the  test  must  be 
conducted  with  the  139.8  MHz  data,  since  this  was  the  only  data  collected 
on  the  same  frequency.  Comparing  the  general  shape  of  the  antenna 
patterns  shows  the  gains  and  losses  are  basically  in  the  same  location. 

The  major  difference  between  the  Griffiss  AFB  data  and  the  Edv’ards 
AFB  data  is  the  magnitude  of  holes  in  the  antenna  patterns.  The 
Griffiss  data  has  holes  of  approximately  twice  the  magnitude  of  the 
Edwards  data.  The  reason  for  the  difference  has  not  been  resolved  at 
this  time. 
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The  equipment  used  during  the  Edv'ards  phase  of  this  test  (figure 
D6)  has  produced  repeatable  results.  With  the  addition  of  airborne 
transmitters,  which  are  capable  of  transmitting  multiple  signals  on 
antennas  simultaneously,  the  equipment  can  provide  an  inexpensive 
permanent  antenna  radiation  pattern  measurement  facility.  (R3) 

CONCLUSIONS  AND  RECOMMENDATIONS 

Orbit  turns  have  shown  to  be  an  effective  flight  technique  to 
produce  continuous  coverage  of  an  aircraft  antenna  radiation  pattern. 
Through  the  use  of  the  muitipath  prediction  option  of  the  Antenna 
Radiation  Pattern  Measurement  Program  (ARPMP) ,  the  effects  of  multipath 
radiation  can  be  effectively  eliminated. 

1.  To  effectively  eliminate  multipath  radiation  effects  in  the  VIIF/ 

UHF  range,  the  multipath  option  of  ARPMP  should  be  used  (page  148). 

2.  Orbit  turns  should  be  used  when  measuring  aircraft  antenna  radiation 
patterns  (page  149). 

The  antenna  pattern  measurement  facilitv  to  be  developed  at  AFFTC 
should  be  patterned  after  the  eouipment  used  in  this  test.  This  facility 
can  produce  accurate,  repeatable  results  v'ith  eouipment  readily  available 
on  base.  With  the  addition  of  airborne  transmitters  to  excite  two  antennas 
with  three  signals  each,  aircraft  antenna  radiation  patterns  can  be 
accurately,  economical Iv  measured, 

3.  An  antenna  radiation  patterns  measurement  facility  should  be 
installed  at  Edwards  AFB  patterned  after  the  system  used  in  this 
test  (  page  150) . 
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OURE  D7  RfiOIOL  FLIGHT  PLOTS  (GRIFFISS  RFB  TEST) 


FIGURE  09  RflOIRL  FLIGHT  PLOTS  (GRIFFISS  OFO  TEST) 
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FIGURE  D13  RADIAL  FLIGHT  PLOTS  CEDHARDS  AFB  TEST) 
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FIGURE  D21  VHF  ANTENNA  PATTERNS  CGRIFFISS  AFB  TEST) 


FIGURE  D22  VHP  ANTENNA  PATTERNS  (GRIFFISS  AFB  TEST) 
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FIOURE  026  UMF  flNTENNR  PfiTTERNS  (ORIFFISS  RFB  TEST) 
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ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AF8  TEST) 


FIGURE  034  VHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 


FIGURE  035  VHF  flNTENNfi  PfiTTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 


FIGURE  036  VHF  flNTENNH  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 


FIGURE  038  VHF  flNTENNR  PATTERNS  HITH  MULTIPATH  CORRECTIONS  ( EOHAROS  AFB  TEST) 
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FIGURE  D42  VHF  ANTENNA  PATTERNS  WITH  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  043  VHF  ANTENNA  PATTERNS  WITH  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  046  UHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  047  UHF  HNTENNfl  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  048  UHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  049  UHF  flNTENNfl  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIOURE  D52  UHF  RNIENNR  PRTTERNS  WITHOUT  MULTIPRTH  CORRECTIONS  ( EOWRROS  RFB  TEST) 


FIOURE  D53  UHF  ANTENNA  PATTERNS  HITHOUT  MULTIPATH  CORRECTIONS  (EDWAROS  AFB  TEST) 


FIGURE  054  UHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  ( EOHAROS  AFB  TEST) 
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FIGURE  055  UHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 


FIGURE  057  UHF  ANTENNA  PATTERNS  WITHOUT  MULTIPATH  CORRECTIONS  (EONAROS  AFB  TEST) 
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FIGURE  062  UHF  flNTENNH  PATTERNS  WITH  nULTIPATH  CORRECTIONS  lEDWARDS  AFB  TEST) 
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FIGURE  063  UHF  ANTENNA  PATTERNS  WITH  MULTIPATH  CORRECTIONS  ( EDWARDS  AFB  TEST) 
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FIGURE  064  UHF  ANTENNA  PATTERNS  WITH  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 


FIGURE  065  UHF  ANTENNA  PATTERNS  WITH  MULTIPATH  CORRECTIONS  (EDWARDS  AFB  TEST) 
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FIGURE  066  UHF  flNTENNR  PATTERNS  WITH  NUlTIPHTH  CORRECTIONS  (EOWRROS  RFB  TEST) 


APPENDIX  E 


MULTIPATH  RADIATION 


TWO-PATH  PROPAGATION  CALCULATION 

The  theory  of  tv.’o-path  propagation  is  thoroughly  developed  in  most 
textbooks  on  electromagnetic  wave  propagation.  This  appendix  was  written 
to  provide  a  condensed  and  simplified  presentation  on  multipath  radiation 
theory  and  to  shov?  the  multipath  prediction  technique  used  by  the 
computer  programs  documented  in  appendices  B  and  C.  Since  a  complete 
development  of  the  theory  is  not  presented  in  this  report,  the  sjTtibology 
used  v/ill  follov'  that  used  by  Reed  and  Russell  (reference  5).  This  v;ill 
alloy;  the  reference  to  be  easily  used  if  a  more  detailed  explanation  of  the 
theory  is  required.  This  section  v;ill  include  four  basic  concepts  v.’hich 
make  up  the  theory  of  multipath  radiation  predictions.  They  are  as 
follovrs:  (a)  earth  gain  factor  (g(0)),  (b)  divergence  factor  (D) ,  (c) 

reflection  coefficient  (R)  of  the  reflecting  surface,  and  (e)  the 
geometric  relationship  betv:een  the  receiving  antenna,  the  transmitting 
antenna,  and  the  reflecting  surface. 

EARTH  GAIN  FACTOR 


The  earth  gain  factor  g(0)  is  the  ratio  of  the  resultant  vector 
(sum  of  the  direct  and  reflected  signal)  to  the  direct  signal. 


g{0) 


(1) 


>=  resultant  signal* 


Eq  *  direct  signal 


The  relationship  of  the  tv'o  vectors  is  shov^n  in  figure  El. 


*TJotel  Actors  are  represented  by  Italic  characters  in  this  report. 
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Figure  El  Rewultant  of  Direct  and  Reflected  Vector 


Where  DREq  ie  the  reflected  signal,  D  is  the  divergence  factor,  R  is 

the  reflection  coefficient  (R  |^|),  6  is  the  phase  change  of  the 
reflected  signal  due  to  the  difference  in  path  length,  and  -<>  is  the 
phase  change  due  to  the  reflection. 

By  the  law  of  cosines 

eJj  -  Eq*  +  (DREq)*  -  2EqDREq  cob  (3)  (2) 

since 

B  »  180  -(0  -  ♦) 

and 

cos  (B)  ■  cos  (180  -  (0  -  ♦))  -  -  cos  (0  -  ♦)  (3) 

then 


therefore,  the  earth  gain  factor  is  equal  toi 


This  is  the  predicted  gain  due  to  multipath  signals,  which  in  db  is 
equal  toi 

g(0)db  ■  10  9(0)  (0) 
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DIVERGENCE  FACTOR 


The  divergence  factor  (D)  is  the  ratio  of  the  spread  of  a  bean 
reflected  off  a  plane  surface  to  the  spread  of  a  bean  reflected  off 
a  spherical  surface.  When  a  diverging  bean  is  reflected  off  a  plane 
surface,  the  direction'  of  the  bean  is  altered  but  the  rate  of 
divergence  is  left  unchanged  (figure  E2).  For  exanple,  a  direct  bean 
and  a  reflected  bean  which  travel  the  sane  distance,  as  neasured  along 
the  bean's  axis,  will  diverge  to  the  sane  size. 

This  is  not  true  when  the  reflecting  surface  is  a  sphere.  Since 
the  angle  of  incidence  is  eqxial  to  the  angle  of  reflection,  a  curved 
surface  will  cause  an  Increase  in  the  rate  of  divergence  of  a  reflected 
bean  (figure  E3}. 


a*0 


✓ 

Figure  £2  Plane  Surface  Divergence 


Figure  E3  Curved  Surface  Divergence 
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A  complete  development  of  the  divergence  equations  can  be  found 
in  Reed  and  Russell  (reference  5) .  The  equations  are  developed  based 
on  comparing  the  cross  sectional  areas  of  the  diverging  beeuns.  The 
ratio  of  the  areas  will  reduce  to  (figure  E4): 


D 


rd  tan 


(7) 


When  d  <<  r  then  r^^  “  h^/sin  X2  -  h^/sin  and  for  small  grazing 
angles  sin  t{i  tan  <{/  which  yields 


D  = 


IT. 


Ih'h^ 


(8) 


rd  tan  *  tj/ 

The  final  usable  equation  is  obtained  by  the  approximation 


^2 


“  52(S6d 


(9) 


and  by  applying  4/3  earth  radius  to  eliminate  the  atmospheric  refraction 
(r  -  5280  miles),  expressing  h  in  feet,  and  d  in  statute  miles.  This 
reduces  the  divergence  equation  to 


P  * 


1  + 


2^2 ' 
dh' 


~h 


(10) 


Therefore,  this  equation  can  be  used,  for  small  grazing  angles, 
to  predict  the  approximate  reduction  in  energy  of  a  signal  reflected 
off  a  spherical  surface. 


REFLECTION  COEFFICIENT 

The  derivation  of  the  reflection  coefficient  (R)  is  more  complex 
than  the  divergence  factor  and  a  thorough  discussion  will  not  be 
attempted  in  this  report.  Rather,  the  basic  equations  will  be 
presented  and  if  a  more  detailed  explanation  is  required,  Reed  and 
Russell  may  be  referenced. 


The  reflection  coefficient  (R)  is  the  ratio  of  the  reflected 
ray  (E  )  to  the  incident  ray  (E.)  and  is  a  complex  vector.  It  can  be 

shown  for  a  vertically  polarized  ray  this  ratio  is  equal  to 


^  n*  sin  i);  -  V  -  cos^ 
^  sin  +  V  m^^^cos^^ 


(11) 
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where  is  the  angle  of  incidence  and  n  is  the  index  of  refraction. 

For  a  horizontally  polarized  ray,  thin  ratio  is  eoual  to 

R„  -  ^  =  sin  »  -  V  -  cos‘»  ,  iRulel'*!,  (12) 

sin  +  \  -  cos^\|i 

The  reflected  rav  can  be  expressed  as  =  E  lRle^*^  where 

/i  *  A.  <, 

$  is  the  angular  phase  change  in  the  reflected  ray.  When  ({>  is  positive 
E^  will  lead  E^  and  when  it  is  negative  E^  will  lag  E^.  By  convention 

(J)  is  normally  expressed  as  the  angle  of  lag  even  though  it  may  be 
over  180  degrees. 

To  complete  the  evaluation  of  the  reflection  coefficient,  the 

index  of  refraction  must  be  calculated.  It  can  be  shov.n  that  n^, 
expressed  in  MKG  units,  is  eoual  to 


= 


DO- 

'r2  toe 


O 


(13) 


'rl 


Eq  =  permittivity  of  free  space  1/(36ti  x  10’)  farad  per  meter 


=  relative  permittivity  of  air  =  1 

^r2  “  relative  permittivity  of  reflecting  surface 

O2  =  conductivity  of  reflecting  surface  MHO  -  meter/square  meter 

0)  =  2nf  angular  frequency  in  radians/second 
With  e^j^  =  1  and  frequency  f  in  MHz,  the  equation  reduces  to 

=  e^2  -  3  ^8  ^  ^0’°2/^MHz 

Some  typical  values  for  e^2  range  are  (reference 


Material 

!2 

sea  v;ater 

81 

3  to 

5 

fresh  v/ater 

81 

10"^ 

to 

m 

1 

0 

wet  earth 

5  to  30 

10"^ 

to 

10"^ 

dry  earth 

2  to  5 

10"^ 

to 

lo"^ 
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GEOMETRIC  RELATIONSHIP 


To  use  the  above  equations  fivn  parameters  must  be  calculated 
from  the  geometric  relationships  involved  in  the  problem.  These 
parameters  are  dj^»  d2t  h^r  6,  and  <{>  as  shown  in  figure  E4,  with  6  being 

the  difference  in  physical  path  length  of  the  direct  ray  r.  and 
reflected  ray  r^. 

Since  the  earth  radius  is  much  greater  than  d. ,  the  angle  i{i  can  be 
approximated  by 


Ip  s:  tan 


5280d2 


where  d^  emd  d2  are  measured  in  statute  miles  and  h^^,  Ah^,  h^r 
Ah2  are  measured  in  feet  with  ■  hj^  -  Ah^^  and  h'  ■  h2  -  Ah2,  From 
figure  E4 

cos  a  ■  »■■  ■  V—  It."  ■ 


which  gives 

ka  +  Ahj^  ■  ka/cos  a 

additionally 


Using  the  series  approximation  for  cos  a  which  is 
cos  a  =  1  -  (^)  +  (^)  •  •  •  • 


and  tidcing  the  first  two  terms  will  give 


cos  a 


cos  (^) 


^1  ' 


therefore 


^1' 

ka  +  -  ka  +  2^ 


'1  21ur 
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Expressing  Lh^  in  feet;  in  statute  miles;  and  ka,  using  the  4/3 

earth  radius  factor  to  correct  for  atmospheric  refraction,  as  5280 
statute  miles  the  above  equation  will  reduce  to  the  final  form  of  Ah. 
as 

^1* 

Ahj^  -  (19) 

A  similar  technique  will  give 

Ah2  »  (20; 

To  develop  the  equation  for  the  physical  path  length  difference 
the  quadrilateral  AAj^Bj^B  in  figure  E4  is  approximated  by  the 

trapezoid  AAj^B^^B  in  figure  E5.  The  difference  in  the  length  of  the 

direct  ray  r,  and  the  reflected  ray  r2,  in  degrees  of  a  wave  length 
is  given  by 

0  =  ^(r2  -  r^)  (21) 


where 


X  *  wave  length  of  signal 

The  length  of  both  rays  can  be  calculated  from  triangles  ABfB  and 
A2B2B  which  gives 


^1 

/  d*  +  (h'  -  hp* 

(22) 

^2 

/  d*  +  (h'  +  hp* 

(23) 

These  two  equations  can  be  expanded  using  the  binomial  expansion 
and  since  d  >>  h^  +  h^  the  higher  order  terms  can  be  neglected  which 
gives 

(h;  -  hn* 

ri  =  d  +  % J  ■  -  (24) 


^2 


d  + 


(h2  +  hp^ 
- 33 - 


(25) 


therefore 


^2  -  ^1 


2h'h2 


(26) 
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Figure  E5  Patii  L  ength  Difference 
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Which  after  converting  units:  h£  and  h^  to  feet,  d  in  statute  miles, 
and  f  in  MHz  gives  the  final  working  equation 


1.385  X  10"^  hrh;f„„ 

6  -  - a - LJJ«i£  degrees 

“miles 


(27) 


The  fii  al  equation  required  to  complete  the  solution  of  the  multipath 
effect  Is  the  relationship  between  h^  and  hi  (figure  E4).  This  can  be 
upproximated  by  the  following: 

Since 


ka  » 


“  ^1^2 

d2  =  Ep4 

and  angle  ^2  degrees;  therefore,  by  similar  triangles 

d. 


^2 


31 


(28) 


SOLUTION  OF  EQUATIONS 

The  solution  to  these  equations,  as  performed  by  program  MULT 
(appendix  C),  is  started  with  the  calculation  of  the  reflection  distance 
which  is  necessary  to  solve  the  geometric  problem  at  the  given 

distance  and  height  of  the  aircraft.  The  technique  used  is  an  iterative 
procedure  of  incrementing  d^  and  solving  for  h2  until  h2  equals  the 

aircraft's  altitude.  The  following  equations  are  used  during  the 
determination  of  d^: 


^2  "  ^2  ^^2 


(19) 

(20) 


(26) 
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Vrtien  h2  best  approximates  the  altitude  of  the  aircraft  the 
final  gain  is  calculated  v/ith  the  folloving  equations: 

-1  ‘^1 

1.385  X  10-“  hp.' 


j  18  X  10" 


(15) 


(27) 

(14) 


if  tlie  calculation  is  for  horizontal  polarization  use 


n  _  sin  Ip  ~  \  11^  -  cos^  Ip 
sin  Ip  +  -co?^~ijr 

if  the  calculation  is  for  vertical  polarization  use 


R  = 
u 


n"  sin  Ip 


-  yj  n  ‘^  - 


cos  Ip 


sin  Ip  +  ^  -  cos^  ip 


-1 

=  """ 

e 


(12) 


(11) 


where  is  the  imaginary  portion  of  R  and  Rg{R}  is  the  real 

portion  of  R. 


R  =  |R| 


D  = 


1  + 


2did2^ 


1 

'3 


g  (6 ) dj^  =  10  log  (  ^  1  +  (DR)=‘  +  2DR(e  -  ({.)  ) 


(10) 

(6) 


These  equations  can 
producing  a  graph  of  t)ie 


be  olved  for  each  range  and  altitude  desired, 
predicted  multipath  propagation  effect. 
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APPENDIX  F 
ARPMP  SOURCE  LISTING 
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1 

5 

10 

15 

20 

25 

30 

35 

AO 


A5 

50 


1 

5 

10 

15 

20 

25 

30 

35 

AO 


C 

c 

c 

c 

c 

c 

no 


120 


130 


lAO 


150 

16C 

C 

170 

180 

IOC 

200 

21C 

22C 


aVER(.AY(AOAI><0>0) 

fROGRAH  AOAPI INPUT, OUTPUT.  TAP f5«INPUT»TAPE6»OUTRUT,TAR£10*TAPEll* 
1  TAPE12) 

CONNON  TAGIlOOIfSIGNL (10C)tDPANG(100),A$PANG(  100) , TSIGC 720) , 

1  ANGAR(720),RANGAR(720),i5aTA(8),IBCD(  7),XARRAY( 20),YARRAY(20)> 

2  CNT(720),ATAG,RNG( 100)>Tin(100),STin,SPTin,SIG(2JOO) , 

3  R ANGE ( 2000 ), TINE( 2000), NX$IG,NXRNG,LL,RN INS.  RNAXS,S0EV( 720), 

A  ASP(720),NINSIG,NINRNG 

THIS  IS  AN  ANTENNA  DATA  ANALYSIS  SOFTWARE  PACKAGE, 

OPTIONS  AVAILABLE  TO  THE  USER  ARE: 

A)  PLOT  DATA  ON  POLAR  PLOTS  OR  RANGE  PLOTS 

B)  NATCH  PERNITTIVITIES  OF  PREDICTED  DATA  AND 

ACTUAL  FLIGHT  DATA 

C)  PREDICT  signal  GAIN  DUE  TO  NULTIPATH  EFFECTS 
O)  HERGE  UP  TO  Three  RADAR  TAPES  TOGETHER 

READ  ADAP  CARO  NO.  1 

READ  (5,180)  ITYPE 

IF  (E0F(5) .NE.O)  GO  TO  160 

ICK-O 

IF  (ITYPE. NE.AHPLOT)  GO  TO  120 
ICK*1 
IITE  (6.2 


WRITE  (6,200) 
CALL  OVERLAY 
GO  TO  150 
IF  (ITYPE. NE. 
ICK-1 

WRITE  (6,210) 
CALL  OVERLAY 
GO  TO  ISO 
IF  (ITYPE. NE. 
ICK.l 

WRITE  (6,220) 
CALL  OVERLAY 
GO  To  ISO 
IF  (ITYPE. NE. 
ICK«1 

WRITE  (6,140) 
CALL  OVERLAY 
IF  (ICK.NE.O) 
WRITE  (6,170) 
STOP 


(6HHPL0TS,1,0) 
5HNATCH)  GO  TO  130 

( SHNATCH,2,0) 
4HHULT)  GO  TO  lAC 

(6HNULT, 3,0) 
SHNERGE)  GO  TO  15C 


(7HRDN  :RGE,A,0) 
GO  TO  160 
ITYPE 


FQRHAT  (IH  ,16H  PROGRAN  NANE  ,AS,21H  IS  NOT  A  LEGAL  NANE  ,//,39h 

ilegal  names  are  merge, plot, hatch.nult) 

Nrnat  (Aid) 


format  (31H 
FORMAT  (27H 


RADAR  MERGE  PROGRAM  IS  CALLED  ) 
PLOTTING  PROGRAM  IS  CALLED) 


FORMAT  (37H  MULTIPATH  MATCHING  PROGRAM  IS  CALLED) 
-  - -  “I  PATH  PREDICTION  PROGRAM  IS  -  ' 


FORMAT  (60H  MULT 
END 


CALLED  ) 


10 

18 

AO 

50 

60 

70 

60 

9C 

100 

110 

120 

130 

lAO 

15C 

160 

170 

180 

190 

200 

210 

220 

118 

250 

260 

270 

280 

290 

300 

310 

320 

330 

3A0 

350 

36C 

370 

380 

390 

AOO 

AlO 

A2C 

A3C 

AAC 

A5C 

A60 

A70 

A30 

A90 

500 

510 

k- 


c 

c 

c 

10 


20 


30 


AO 

c 

c 


OVERLAY(MPLDTS,  1,0) 

THIS  IS  THE  MASTER  PLOT  SUBROUTINE  WHICH  INITIALIZES  DATA  AND 
CALLS  either  THE  POLAR  PLOT  PROGRAM  PPLOTS  OR  THE  RANGE  VS 
SIGNAL  PLOT  program  SPLITS 

PROGRAM  MPLOTS 

COMMON  TAG(IOO) .SIGNL (100),DPANG(100)>ASPANG(  100 ) , T n G( 720 ) , 

1  ANGAR(720),RANGAR(720) , ID AT A( 8) ,  I BC D (  7 ) , X ARRA Y ( 20 ) , YA R R A Y ( 20 ) , 

2  CNT( 720), ATaG,RNG( 100 ) , T  IH ( 100)  ,  ST  IN,  SpT IM, S IG (  2000 ) , 

3  RANGE! 2000), TIME (2000),NXSIG,MXRNG,LL,RN INS, RMAXS.SDEVI 720), 

A  ASP(720)f HINSIG.HINRNG 

REAL  MXSIG.MXRNG 
CALL  PLOTS  (DUN, DUN, 13) 

IP  ITYPE  IS  EQUAL  TO  P;  THEN  POLAR  PLOTS  WILL  BE  PRODUCED 
IF  ITYPE  IS  EQUAL  TO  S;  THEN  RANGE  VS.  SIGNAL  PLOTS  WILL 
BE  PRODUCED 

LL-1 

REWIND  12 
DO  2C  M»1 

tagTm)«o. 

SIGNL (M) • 

OPANG(M)* 

ASPANG(n) 

RNG(M)*0. 

TIM(H).0. 

CONTINUE 
DO  30  M*1 
S0EV(M)-0 
ASP(M),(M 
TSIG(M)«0 
CNT (M)»0. 

ANGAR(M) • 

RANGARI M) 

CONTINUE 
DO  AO  M,1 
SIG(N),0. 

RANGE (M)- 
TIME(M)»0 
CONTINUE 


-  ICO 

0 

0.0 

0.0 

■0.0 

0 

0 

,720 

0 

-!)♦.! 

.0 

0 

0.0 

•0.0 

^  2000 
0 

86° 


1C 

20 

30 

AO 

50 

60 

70 

80 

90 


B  100 
B  lie 
0 
B 
B 
B 
B 

a 


ITYPE,  the  start  AND  STOP  TIK'i  TO  BE  PLOTTED  ARE  READ  IN 


H8 

lAO 
15C 
160 
170 
8  180 
B  190 
B  200 
B  210 
B  220 
B  230 
B  2aO 
B  259 
a  260 
B  270 
a  289 
B  290 
B  300 
B  310 
a  320 
a  330 
B  3A0 
B  359 
B  360 
a  370 
a  380 
B  399 
B  AOC 
B  AiO 
3  A20 
A39 


3  AAO 


22( 


so 

55 

60 

65 

70 

75 

1 

5 

10 

15 

20 

25 

30 

35 

60 

65 

50 

55 

60 


60AP  CARO  NO.  3 

READ  (5>60>  ITYPE>IHRl,flINl,ISECl>IHR2>HIN2>ISEC2 
IF  (E0F(5).NE.0.0!  60  TO  50 

THE  SOURCE  NUHB£R»THE  Hlh  AND  MAX  SI6HAL  AMO  RANGE:  AND 
THE  NIN  ANO  MAX  VALUE  TO  3E  USED  FOR  WILD  POINT  CHECK  IS 
READ  in:  AOAP  CARO  NO.  6 

READ  (5f70)  ATAGtNINSIG.nxSie.NINRNGt NXRN6tRNINS»RNAXS 


KCAU  I3f  rui  ■  I  nbt  nin> ibt no  n«Kni>t  I'ninit  Kn*A2 

WRITE  (6.80)  IHR1>NIN1«IS£C1>IHR2>NIN2.ISEC2> ATAG.niNSIG.NXSIG.niN 
IRNG.nXRSG.RHINS.RNAXS 
STin>IHRl«3b30.«niNl*60.«ISECl*l. 
SPTIN>IHR2*1600.PHIN2*60.»IS£C2*1. 

IF  tITYPE.EQ.lHP*  CALL  PPLOfs 
IF  (ITYPE.EO.lHS)  CALL  SPLOTS 
GO  TO  10 

CALL  PLOT  (0.f0..<)99) 

WRITE  (6.90) 

STOP 


FORNAT  (A1.3X.312.6X.3I2) 

FORMAT  (8F10.3) 

format  (IHl.AiHIF  ITVPE  «  P. 


format  (IHl. 
IlX.SaHIF  ITY 
2.1X,17HTm|  S 


68HIF  ITVPE  •  P.  THEM  P 
PE  -  S.  THEN  RANGE  VS. 
TART  TIME  IS.2X.3I2.2X. 


6  OATA  THAT  IS  LCSS 
TILL  NOT  BE  PLOTTED) 


POLAR  PLOTS  WILL  BE  PRODUCED  ./. 
.  SIGNAL  PLOTS  WILL  BE  PRODUCED./ 
X.16HTHE  STOP  TINE  IS.2X.3I2./.1X 

a‘».Kis»!bS'js 


FORMAT  (26H 
ENO 


■NO  OF  DATA  TO  BE  PLOTTED) 


B 

650 

B 

660 

B 

670 

B 

680 

S 

690 

B 

500 

B 

SIC 

B 

520 

B 

530 

B 

560 

B 

55b 

B  576 
B  580 
B  59C 
B  600 

a  610 
a  620 
B  630 
B  660 
B  650 
B  660 
B  670 
B  680 
B  690 
8  TOO 

I 

B  730 
B  760 
B  750 

I  m 

1  780 


SUBROUTINE  PPLOTS 

PO^R  PLOTS  WILL  B|  PRODUCED:  EFFECTIVE  RADIATED  POWER 

wUl  be  plotted  versus  The  azimuth  ancle  of  the  aircraft 

COMMON  TAG(10a).SIGNL(ipC).DPANG(lC0).ASPAN6(100).TSIG(720). 

1  AN6AR(720).RANGAR(72S).10ATA(B).I|c5(7) .XARRAY(20).YARRAY(20)i 

2  CNT(720).ATA6>RN6(100).TIH(10&).STIN.SPTin.SIG(2000)> 

3  RANG|(|^00).T1HE(2000).>1XSI6.HXRNG.LL.RNINS.RMAXS.SOEV(720). 

real'mxsig.mxrng 

the  depression  angle,  signal  interval  to  be  plotted  is  read  I) 

OELD  IS  THE  DEPRESSION  ANGLE  INTERVAL;  IPLANE  IDENTIFIES 


SION  ANGLE.  SIGNAL  INTERVAL  TO  BE  PLOTTED  IS  READ  IN. 
THE  DEPRESSION  ANGLE  INTERVAL;  IPLANE  IDENTIFIES 
RAFT  TO  BE  DRAWN  aT  CENTER  OF  PLOT.  IPAPER  IDENTIFIES 
OR  NOT  THE  GRID  ROUTINE  IS  TO  BE  CAUFD  TO  DRAW  THE 
ID;  IBCO  IS  A  70  COLUMN  COMMENT  THAT  IS  DRAWN  AT  THE 


OELD  IS  The  oepres 
THE  AIRCRAFT  TO  BE 

whether  or  not  the 

POLAR  GRID;  IBCO  I 
TOP  OF  THE  PLOT 


)  (5.220)  ANG.SINT.DELD. IPLANE. IPAPER 
IIPAPER.NE.O)  CALL  HOT  (0..S.5.>3) 

i<ang>oel5 


TOP  OF  THE  PLOT 

DIMENSION  SG(720).CT(720).AP(720).S0( 720) 

TLAST«C. 

WRITE  (6.210) 

READ  AOAP  CARO  NO.  5 

READ  (5.230)  IBCO 

READ  AOAP  CARO  NO.  6 

READ  (5.220)  ANG.SINT.DELD. IPLANE. IPAPER 
IF  (IPAPER.NE.O)  CALL  HOT  (0..S.5.>3) 

ANG1<ANG>0EL0 

ANG2«ANGp0EL0 

RBIAS«(6-(NXSIG/SINT))*SINT 

DATA  IS  READ  FROM  OUTPUT  FILE  OF  ARPMP 

READ  (12)  (TAG( I).SI6Nl(I).0PANG(  I)>ASPANC(I)fRN6(I).TIM(I).I«l.S0 
D.WDnO 
NWO<NONQ/6 

IF  (EOF( 12).NE.0.0)  GO  TO  190 
TLAST-TIMINWO) 

DATA  IS  SUMMED  AND  STORED  FOR  EACH  DATA  CELL 
no  3C  J>1.NWQ 

IF  (STIH.GT.TIN( J) )  60  TO  3C 
IF  (TIM(J).GT.SPTIN)  GO  TO  60 
IF  (AN6.EQ.999.  )  GO  TO  20 

IF  (DPANG(J).LT.AN61.0R,0PAN6(J).6T.AN62)  GO  TO  30 
IF  (ATAG.NE.TAGIJ))  GO  TO  30 

IF  (signl(j).lt.rmins.or.signl(j).gt.rmaxs)  go  to  30 

!i>ASPANC( J)*2»1.S0 
IF  (II. GT. 7201  II»1 
CNT(II)-CNTTII)*1 


CNT(II)-CNTni)*l 
TSi6(il)«TSIG(  in^SIGNKJ) 
SDEV(l!)-SDEV(il)«SIGNL(J)*SI( 


CONTINUE 
GO  TO  10 

DATA  IS  AVERAGED  ANO  THE  STANDARD  DEVIATION  IS  CALCULATED 
DO  50  J-1.720 


[II 

C  360 

^  4 

C  370 
C  380 
C  390 

uu 

C  630 

C  68y 
C  690 

m 

C  560 

I  in 


I  ^18 


»7 


1 


IF  (CNT(  J) .EQ.O)  GO  TO  50 

1 

C  640 

65 

TSIG( Ji*TSIG( J) 7CNT( J  >  — 

- 1 

:  650 

IF  (CNTI jT.IT.2)  S0EV(J)«0. 

1:  660 

IF  (CNT( J).LT.2)  GO  TO  50 

:  670 

$OEV(  J)>SORT((SDEV(  J1-TSIG(  J)*TSIG(  J)*CNTt  Jn/(CNT(  J)-l)  ) 

:  680 

50 

CONTINUE 

:  690 

70 

L1>0 

:  700 

C 

c 

DATA  CELLS  WITH  NO  DATA  ARE  THROWN  OUT 

1  m 

c 

:  730 

DO  60  J«l(720 

IF  (CNT( J) .EQ.0.0)  GO  TO  60 

L1-L1*1 

‘  740 

75 

:  76? 

SGIL1)«TSIG( J) 

:  770 

CT(L1)»CNT(J) 

AP(Ll)*ASP(ji 

:  780 
:  790 
:  800 

8C 

SD(L1)«SDEV(  J) 

60 

CONTINUE 

:  810 

IF  (Ll.LT.ZOO)  GO  TO  70 

.  820 

Ki.l 

:  830 

K2»|00 

:  840 

85 

GO  TO  80 

1  850 

70 

Xl-l 

:  860 

K2*L1 

:  870 

C 

:  880 

C 

oata  to  be  plotted  is  printed  on  line  printer 

I  "890 

<)0 

C 

:  900 

80 

DO  90  NN>1*4 

:  910 

WRITE  (6t240)  IBCO.Ll 

WRITE  (6«250)  ANG.DELO* SINT> I  PLANE.  I  API  J )  ,  SGI  J ) , SD 1 J ) f C T ( J ) > J 

‘  920 

XI, K 

:  930 

12) 

:  940 

95 

K1«K1*200 

:  950 

K2.K2+200 

:  960 

IF  (K2.GE.L1)  GO  TO  100 

IF  (K2.GT.720)  K2-720 

- 

:  970 

:  980 

100 

90 

C 

CONTINUE 

:  990 

:iooo 

C 

IF  IPAPER  NOT  EOUAL  TO  ZERO;  POLAR  GRID  IS  DRAWN  BY 

:ioio 

C 

SUBROUTINE  GRID 

:i020 

C 

:io3o 

105 

100 

IF  (  IPAPER. NE.O)  CALL  GRID 

:io4o 

:io5o 

ioaTa(1)«ioh  azimuth  I 

IOATA(2)-10HS  RELATIVE 

:i060 

lOATAOl  •IqH  to  A/C  he 

IOATA(4)>iOHAOING 

CALL  NEWPENIZ) 

:io9o 

110 

CALL  SYMBOL  I  >3 . 5«  4 . 8.  .10.  I  BCD*  0.  >  70 ) 

:iioo 

YPACE»4.0 

:iiio 

110 

CALL  NUMBER  ( . 1 5. YPAGE. . 10. MXSIG. 0. > 2 ) 

:ii2o 

CALL  SYMBOL  I  999 . . 999*. # 10. 4H  D8W.0..4) 

IF  (YPAGE. LE.l.)  GO  TO  IZO 

MXSIG»MXSIC-SINf 

YPAGE-YPAGE-l. 

11130 

:ii4o 

115 

GO  TO  110 

:ii7o 

120 

CALL  SYMBOL  ( -1 . 8.-4. 8. . 10. ID ATA.O . .40) 

hl80 

12C 

c 

CALL  SYMBOL  ( 0 . * 0 .. . 10. 3. 0 . . -1 ) 

11210 

c 

PLANE  IS  DRAWN  aT  CENTER  Of  PLOT  — “ 

c 

CALL  PLANE  (IPLANE) 

C»LL  NEWPENll) 

:i2zu 

iwn 

125 

K>1 

11250 

NZERO.Q 

11260 

130 

00  180  J«1.720 

IF  (CNT( J) .NE.O)  GO  TO  130 

NZERO-NZEROPl 

11270 

11280 

11290 

1x300 

113{0 

c 

IF  5  DEG  ASP  WITH  NO  OATA.  PLOT  DATA  SEGMENT 

11320 

IF  (NZERO.EO.IO)  GO  TO  140 

11330 

GO  TO  170 

11340 

135 

c 

11350 

c 

OATA  IS  PACKED  INTO  THE  PLOTTING  ARRAYS 

11360 

11370 

c 

RANGAR(K)<TSIG(J)PRBIAS 

130 

11380 

ANGAR(K).-1*((J-1) *.5-90)757. 4 

11390 

140 

K.K  +  l 

.1400 

NZERO-0 

11410 

CO  TO  170 

11420 

140 

IF  (K.GE.Z)  GO  TO  150 

11430 

C 

1144Q 

145 

c 

\  POINT  IS  PLOTTED  WITHOUT  CONNECTING  LINE 

11450 

c 

11460 

CALL  POLAR  (RANGAR.ANGAR.K-1.1. -1.3.0. .SINT) 

Jl470 

GO  To  160 

11480 

c 

11490 

150 

c 

POINTS  ARE  PLOTTED  WITH  CONNECTING  LINES 

11500 

c 

11510 

15C 

CAIL  POLAR  (RANGAR.ANGAR.K-1. 1.1. 3.0. .SINT) 

h520 

I6C 

K.l 

11530 

NZERO-0 

11540 

155 

GO  TO  180 

11550 

170 

IF  (J.EQ.720)  GO  TO  140 

11560 

18C 

CONTINUE 

WRITE  (6.260)  LL 

IF  {  {PAPER. NE.O)  CALL  PLOT  ( 14..-5 . 5. -3 ) 

IF  (IPAPER. EQ.O)  CALL  PLOT  (1Z..0..-3) 

11570 

160 

11600 

^^tUn^ 

.  16i0 
fl620 
11630 

190 

WRITE  (6.270)  TLAST 

211 


DO  200  J*l>720 

.NE.O)  GO  TO  «0 


FORNAT 

^ORHAT 


ITVRE  «  R  i  ROLAR  RIOTS  WILL  BE  RRODUCEO  ) 


3F10.3>2I1) 

BAIO) 

lH{,28X.7AlO>/>i 


FORMAT  <lH{,28X.7AlO>/>2eX,14t.22H  ROINTS  TO  BE  RLOTTED  ) 

FORMAT  (IH  •ZOXtSAHRLOTTED  DATA  FOR  DERRFSSION  ANGLE  >Fi0.3>17H  RL 
lUS  OR  MINUS  .F10.3>6H  DEGS.>/»2AH  THE  SIGNAL  INTERVAL  IS  fFI0.3f 
2I0H  IRLANE  -  >I1»/»AC33H  ASR  ANG  SIGNAL  STOEV  COUNT  )>/»200(l 
3X.I6Fe.2|/)l 

FORMAT  (IXtiSHROLAR  RLOT  NUMBER  tlX«I2»9H  RLOTTED  I 

format  (A9H  END  OF  FILE  — TaN  12  LAST  TINE  RROCESSED  WAS  tFlO.3 

4nd 


8HR0LAR  RLOT  NUMBER  tlX«I2»9H  RLOTTED  I 


FORMAT  (IXtlSHROLAR  RLOT  NUMBER 
format  (A9H  END  OF  FILE  — TaN 


»1X«I2»9H  RLOTTED  I 

12  LAST  TINE  RROCESSED  WAS  tFlO.3 


SUBROUTINE  RLANE  ilTVR) 

THIS  SUBROUTINE  DRAWS  THE  A/C  IN  THE  CENTER  OF  THE  ROLAR  RLOT. 

IF  A  different  view  OF  THE  A/C  IS  NEEQEDt  IT  HAY  BE  ADDED  IN 

THIS  ROUnNE  WITH  ITYR  BEING  USED  TO  tOENTlFY  WHICH  A/C  IS  TO 

BE  DRAWN. 

DIMENSION  XARRAV(20i>VARRAY(20) 

XARRAY(n--.9 

XARRAYI2l«-.i666 

XARRAY(3I«0.0 

XARRAY  A)-.1666 

XARRAY(Sj<.S 

XARRAY(6i>.i666 

XARRAY(7I-.5833 

XARRAyIs  <.0833 

XARRAY(9)<».1666 

XARRAY(10j<-.1666 

XARRAY<II)<-.0833 

XARRAY(12)<-.0e33 

XARRAY(l3 j<-.i666 

XARRAY(Ui<-.f 

XARRAY(19I>0.6 

XARRAY(16I<2, 

YARRAY(iT<-.1666 

YARRAY(2)<0. 

YARRAY(3><.5 

YARRAYjAi<0. 

YARRAY  91— .1666 


YARRAY  6  <-.1666 
YARRAYi7i<-.l333 

yarrayIb)*-.? 

YAIRAYi?l  —  .A2 
YARRAyHOI— .62 
YARRAYiIl)<-.A 
YARRAY  12j<-, 3333 
YARRAY|l3i'-.i666 
YARRAY(l6»<-.l666 
YARRAY(19I>0<0 
YARRAY{i6)<2. 

CALL  LINE  (XARRAYtYARRAYf I6fl»0>0) 

RETURN 

END 


SUBROUTINE  SRLOTS 

RANGE  VS.  SIGNAL  (EFFECTIVE  RADIATED  ROWER)  PLOTS  ARE  DRAWN 

COMMON  TAG(100)*SI6NL(100)tDRANG(lQO)>ASRANC(  100) r TSIG( 720) > 

1  ANGAR(720T>SANGAR(72S)»IDATAI2irliC0(7)f XARRAYI25)>YARRAY(20)« 
I  SNT(720),ATAG.RN6J  10O),TIMli0C).STIM,SRTI»,.$iGT2O067..,,,,^, 

3  R ANGE (2000)  fTINET2000)»NXSIGtHXRNG>LL.RMI.YS»RMAXS>SOEV( 720)  > 

6  ASR(720TtNtNStG>HINRNG 
DIMENSION  ST(2l 


DIMENSION 
REAL  MXRN 


MXRNG»MXSIG>MINSIGf HINRNG 


READ  AOAR  CARO  9 
READ(9«190)  IBCD 

READ  AOAR  CARO  NO.  7 


tEAD(9>180> 

ST(i)<HlNRN( 

STi2)<MiNSI< 


XLENtYLEN 


U(2t<M!NSIG 

(INT<(HXRN6-SI(1))/XLEN 


YINT<(HXSIG- 

ICK<0 


r(2))/YlEN 


CONTROL  CAROS  ARE  PRINTED 
WRITE  (6»200>  XLEN.YLEN»ST(1)>ST(2)>XINT 
DATA  IS  READ  FROM  OUTPUT  FILES  OF  ARRMR 
READ  (12)  (TAG( I)>SIGNL(I).DRANG( I)tASRANG(I)»RNG(I)»TIM(I).I>l»80 


229 


45 

50 

55 

40 

65 

70 

75 

80 

85 

90 

95 

100 

105 

110 

115 

120 

125 

130 

135 


20 


30 


40 

V 

60 


90 

100 


110 


120 


1>>WONO 

.NE.0.0)  CO  TO  160 
LAST>TIH(Kt 

DATA  IS  STORED  IN  THE  ARRAYS  SICfRANCE  AND  TINE 

|st](nHt.TIH(JI)  CO  TO  20 
IF  jTlHt JI.CT.SRTIN)  GO  TO  30 
IF  jATAC.NE.tACI J))  GO  TO  20 

IF  (SICNLIJI.IT.RNINS.OR.SICNLIJI.GT.RNAXS)  CO  TO  20 
N-N«l 

SIC(N)>SICNL(Ji 

RANGE(nI>RN6(JI 

TinE(N)«TIN(J> 

IF  (N.EQ. 19981  GO  TO  30 

CONTINyE 

GO  TO  10 

DATA  LESS  THAN  THE  HININUN  SIGNAL  AND  GREATER  THAN  THE  NAXIHUN 
SIGNAL  THAT  HAS  READ  IN  IS  THROHN  OUT 

F  N.6T;n^G0  to  60 

F  (SIG(J).GE.ST(2T.AND.SIG(J).IE.NXSIC>  CO  TO  50 
RITE  (6t210l  StGTj)fRANGE(JI»TIN(J) 


nn«n 

|Q  50  JflfNN 
«R 

JK-J 
KK»N-1 

DO  40  KJ«JK,KK 

RirGEllK  jf  ■^ANGeU  JAl  I 
COUTINUE 
N«N-i 
CONTINUE 

DATA  TO  IE  FLOTTEO  IS  FRINTED 

HRiTi  IliHSi  N?f8lC(IlfRANGEII)f  I>1»N) 
IF  (tCK.EO.l.l  GO  TO  150 

x^agI>xlen 

ftSl  [irdv.zv. 


8all 


21 


X-AXIS  IS  DRAWN 


IF  (XFAGE.LE.O.OI  CO  TO  80 

CALL  SYnIOL  (XFACEf C.».l«l3»0.>-1} 

XFAGE<YFAGE-.3 

GO  TO  VO 

XFAGE<-.15 

X-AXIS  LA8EIING  IS  DRAWN 
VAL'STIll 

CALL  NUNBER  (XFA6E»-.2>.1»VAL*0.>21 

VAL-VALaXINT 

XFAGE«XFAGE4.9 

K>XLEN 

DO  90  I>1*K 

IF  (XFACE.GT.XLEN)  GO  TO  100 

CALL  NUNBER  (XFACE. -.2* .l»VALfO.»2) 

VAL«yAL«XlNT 

XFAGEaXFAGEAl.O 

CONTiNUE 

YFAGE  —  .e 

I0ATA(1I>10HRANGE  NN 
XFACE>XLEN/2.-.S 

CALL  SYNIOL  (XF AGE. YFAGE* .13» 10ATA>0.f 101 
CALI  FLOT  (0..0..3) 

Y-AXIS  IS  DRAWN 

CALL  FLOT  (0.>rLEN.2) 

|F*fYFl^|!!LE.0.0)  GO  TO  120 

CALL  SYNIOL  (0. .VFAGE. .1. 13.90. .-1) 

YFA«E>YFAGE-.5 

GO  TO  110 

XFA6E>-.4 

YFACE-0.0 


Y-AXIS  LABELING  IS  DRAWN 

( XF AGE. TRACE*. l.VAlfO.. 21 


KaIil^nuSUr 

K*YLEN 
VAL«yAL*YINT 
yfage.yfaGEfi, 
XFA6|>XFAGE-.1 

?F° 


IF  lYFAGE.GT.YLEN)  CO  TO  140  . 

CALL  NUNBER  (XFAGE* VFAGE. .l.VAl.O. *21 

val-valayint 


130 

140 


YFACEAl.O 


Toftl  .. 

XFAGE-XLEN;2.-6. 


NflNUE 


360 

370 

380 

390 

400 

410 

420 

430 

tn 

460 

tn 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 


i 


iff] 


11335 


2M 


IF  (XPAGE.lt. 01  XPAGE-0. 

YPA6E«rLENA.AO 

CALL  SYNBOL  ( XP AGE. VPAGE» .1 S. IBCD>0 .» 70) 


CALL  SYNBOL  (XPAGE.VPAC 
IDATA(1)>10HSI6NAL  OBW 
XPAGE  — .8 


YPAGE«3.9 

CALL  SYNBOL  ( XP AGE. YP A6E» .15. IDATA# 90 .. 10) 

FIRST  POINT  AND  SCALING  DATA  IS  LOADED  INTO  THE  PLOTTING 
ARRAYS  FOR  CALL  Of  LINE  ROUTINE  TO  DRAW  THE  DATA  CURVE 


190 

L 

c 

c 

199 

c 

160 

c 

169 

160 

170 

oooo 

179 

210 

180 

220 

230 

29C 

189 

290 

RANGE(NP1)-SI(1) 

RANGE(N4^2)-XINT 

SI6(N*1)«ST(2) 

SI6(N«2)-YINT 


DATA  CURVE  IS  DRAWN 

CALL  LINE  (RANGE. SIG.N.l.O.C) 

KN«HDN0/6 

N«0 

CHECK  IS  NADE  TO  SEE  IF  NORE  TIN 
TO  BE  PLOTTED--IF  SO  CONTROL  IS 


,  _  TO  SEE  IF  NORE  TINE  IS  LEFT  IN  THE  TINE  INTERVAL 

to  iE  Plotted— IF  so  control  is  transferred  back  to  statenent 

10  and  NORE  DATA  IS  READ  FRON  THE  OUTPUT  FILE  OF  ARPNP 

IF  (TIN(KN) .LT.SPTIN)  60  TO  10 
WRITE  (6.290)  LL.TLAST 
LL*LL*1 

CALL  PLOT  (XLENP6..-1..-3) 

RETURN 

WRITE  (6.2S0)  TLAST 
STOP 

FORHAT  (9QH  ITTPE.  Si  RANGE  VS  SIGNAL  PLOTS  HILL  BE  PRODUCED  ) 
FORNAT  (BF10,3) 

FORHAT  (BAlOl 

FORHAT  (26H  THE  THE  X-AXIS  LENGTH  IS  .F10.3.21H  THE  Y-AXIS  LENGTH 
lIS«F10.3./.lX«22HTt<£  STARTING  RANGE  IS  .F10.3.2CH  THE  START  SIGNAL 
2  IS>f10,3./.1X.29H  THE  X  AXIS  INTERVAL  IS  .FIO.B) 

FORHAT  (IH  .90H  THE  DATA  POINT  THAT  IS  OUT  OF  RANGE  IS  .F10.3.19H 


FORNAT  (BF10,3) 

FORHAT  (BAlOl 

FORHAT  (26H  THE  THE  X-AXIS  LENGTH  IS  .F10.3.21H  THE  Y-AXIS  LENGTH 
lIS«F10.3./.lX«22HTt<£  STARTING  RANGE  IS  .F10.3.2CH  THE  START  SIGNAL 
2  IS>f10,3./.1X.29H  THE  X  AXIS  INTERVAL  IS  .FIO.B) 

FORHAT  (IH  .90H  THE  DATA  POINT  THAT  IS  OUT  OF  RANGE  IS  .F10.3.19H 
lAT  A  RANGE  OF  .F10.3.19H  AT  A  TINE  OF  >F10.3) 

FORHAT  (1H1.20X.8A10) 

FORHAT  (1H0.39X.39H  RANGE  VS.  SIGNAL  ANTENNA  PATTERN  PLOT  ./.39X.3 
19H  THE  NUNBER  OF  DATA  POINTS  TO  PLOT  ARE  . 110. /. 39X.21H  THE  DATA  A 
IrRAYS  are  f //.1X.6(2X.6HSIGNAL.9X.9HRANGE.2X)»//.170I1X.12F10.3./) 

^FORHAT  (13H  PLOT  NUNBER  .I10.33H  PLOTTED  LAST  TINE  ON  PLOT  WAS  . 

^pi^AY  (IH  .99HEN0  OF  FILE  ENCOUNTEtEO  TAPE  12  LAST  «NE  PKOCtSSEO 
1  WAS  .F10.3) 

END 


£1990 

E1960 

E1970 

Im 

E1900 

E1930 

E1990 

£1990 

Wo 

£1980 

H990 

E1600 

E1610 

E1620 

£1630 

tl690 

E1690 

£1660 

E1670 

£)6e0 

£1690 

a7oo 

m 

E1730 

£1790 

ii??8 

El  770 
E1780 

iiut 

asio 

im 


SUBROUTINE  GRID 

THIS  SUBROUTINE  DRAWS  A  POLAR  COORDINATE  GRID 

DIHENSION  RAD(2).ANGL(2).IBCD(a) 

A  CALL  IS  HADE  TO  FUNCTION  NEWPEN  TO  DEFINE  A  NEW  PEN  FOR 
USE  IN  DRAWING  A  POLAR  GRID  (THUS  A  DIFFERENT  COLOR  HAY 
USED  FOR  THE  GRID) 

CALL  NEWPEN  (2) 

RA0(1)>.9 

RA0(2)<i9.9 

ANGL(l)-80./97.9 

ANGL(|)<80./97.9 

0EG«10. 

RAOULS  ARE  DRAWN  EVERY  10  DEGREES 
DO  10  Ulr36 

CALL  POLAR  (R AD. ANGL. 2.1.0.0.-9..  1.  ) 

XPAGE«9.6*SiNr90  /97.9-ANGL (h ) 
yPAGE>9.6*C0S(90./97.9-ANGL(l )> 

IF  (DEG. GT. 180. I  XPAGE>XPAG£-.93 

CALL  NUNBER  (XPAGE. YPAGE..10.DEG.O..O) 

ANGL(iT'AN61(1)-10./97.9 


CALL  NUNBER  (XPAGE. YPAGE..10.DEG.O..O) 

ANGL(iT-AN6i(1)-10./97.9 

ANGL(2)-ANGL<1> 

OEG«OEGPlO 
CONTINUE 
XPAGE*. 9 
Ra>.9 
RF».9 

DO  20  1*1.9 

A  CIRCLE  IS  DRAWN  EVERY  HALF-INCH 

CALL  CIRCL  (XPAGE. 0..0.. 360. .RO.RF.O.) 

XpX6E*XPAGE*.9 

R0*R0*.9 

RF*RF*.S 

CONTINUE 

A  BORDER  IS  DRAWN  ABOUT  THE  PLOTS 
CALL  PLOT  (-6. .-9. .3) 


F  225 
F  230 
F  290 
F  290 
F  260 

^  Hi 

F  290 
F  300 
F  310 
F  320 
F  330 
F  390 
F  390 
F  360 
F  370 
F  380 
F  390 
F  900 
F  910 
F  920 


211 


50 


1 

5 

10 

15 

20 

25 

30 

35 

^0 


«>3 


5C 

55 

60 


65 


70 

75 

80 


85 

90 


CALL  PLOT  (-6. >5. >2) 

F 

660 

CALL  PLOT  (6.>5.f2) 

f 

670 

CALL  PLOT  (6.>-5.*2) 

f 

CALL  PLOT  (-6.,-5.»2) 

CALL  NEWPEN  Tl) 

F 

690 

F 

S95 

RETURN 

F 

515 

END 

f 

520' 

0VERLAY(NATCHf2(0) 

G 

10 

c 

G 

20 

c 

OVERLAY  NATCH  OETERNINES  THE  PERNITTIvnY  WHICH  PRODUCES  THE 

G 

^8 

c 

BEST  ELININATION  OF  NULTIPATH  SIGNALS 

G 

c 

G 

50 

PROGRAM  HATCH 

G 

60 

REAL  LASTVR 

G 

70 

COMPLEX  REF 

G 

80 

DIMENSION  SI6(101)»PMULT(101I>IHEAO(8)>CNT(1D1).ISTR(3I>ISTP(3)# 

1  OSR(80)»OSIG!80l«OEP(80l»ASP(80)»RNG(eO)>TtN(80)>OIF(101) 

WRITE  I6«160) 

6 

90 

G 

100 

G 

110 

10 

REWIND  11 

G 

il8 

PERH>1. 

G 

C0N0>.C01 

G 

160 

ENT«1. 

G 

150 

c 

G 

160 

c 

READ  IN  HEADING  INFORMATION 

G 

170 

c 

G 

18C 

c 

READ  AOAP  CARO  NO.  10 

G 

195 

c 

READ  (5.1901  IHEAO 

G 

G 

m 

IF  (EOF(S).NE.O)  GO  TO  170 

G 

220 

WRITE  (6.200)  IHEAO 

G 

230 

c 

G 

260 

c 

READ  IN  START  TIME ( HR. MIN. SEC) . STOP  TIHE(HR. MIN.  SEC )> HEIGHT 

OF  antenna  above  reflecting  surface  in  feet.  HEIGHT  OF 

G 

c 

G 

265 

c 

AIRCRAFT  ABOVE  ANTENNA  IN  FEET.  FREQUENCY  IN  MHZ.  AND 

G 

c 

SOURCE  NUMBER 

G 

285 

c 

G 

290 

c 

c 

READ  AOAP  CARO  NO.  11 

G 

G 

l?8 

c 

G 

320 

READ  (5.210)  ISTR.lSTP.HT.ACHT.FREO.tSR 

G 

330 

SR»1SR 

G 

365 

c 

G 

350 

c 

calculate  the  START  AND  STOP  TINES  IN  TOTAL  SECONDS 

G 

365 

c 

G 

370 

START-ISTR(1)*3600.*ISTR(2)*6C.6ISTR(3) 

STOP.ISTP(l)*3600.*ISTP(2)*60.6lSTP(3) 

G 

380 

G 

390 

WRITE  (6.220)  ISTR. ST  ART. ISTP. STOP.  HT.ACHT. FREQ.  ISR 

G 

600 

DO  20  I«l.l5l 

G 

cnt7i)*o 

G 

620 

SIG(I)»0 

G 

635 

20 

PMULT(I)«0. 

G 

660 

c 

G 

650 

c 

READ  DATA  FROM  OUTPUT  FILE  OF  ARPMP 

G 

660 

c 

G 

670 

30 

READ  (11)  (OSR(  I).0SIG(1).0EP(1).ASP(I).RN6(I).TIM(  I)»1>1.80).ANUM 

G 

680 

IF  (EQF(ll).NE.O)  GO  TO  160 

G 

690 

NUM«ANUM/6. 

G 

505 

c 

G 

510 

c 

CHECK  IF  DATA  IS  WITHIN  TIME  INTERVAL  AND  PROPER  SOURCE  NUMBER 

G 

320 

c 

G 

^30 

DO  60  I«1.NUH 

G 

565 

IF  (TIM(I).LT. START)  GO  TO  60 

G 

550 

IF  (riM(I).GE.STOP)  GO  TO  50 

G 

560 

IF  (SR.NE.OSR(I))  GO  TO  60 

G 

570 

c 

G 

580 

c 

CALCULATE  RANGE  INDEX  TO  STORE  THE  SUM  OF  THE  SIGNALS 

G 

390 

c 

G 

600 

lI»(RNG(l)-3,8)/.6»l 

G 

610 

IF  (II.LT.l.OR.II.GT.lOl)  GO  TO  <0 

SIG(II)«SIG(II)aOS(G(I) 

G 

620 

G 

635 

CNT(II)«CNT(Il)tl 

G 

660 

60 

CONTINUE 

G 

650 

GO  TO  30 

G 

660 

50 

DO  60  I'l.lCl 

G 

670 

c 

G 

680 

c 

CALCULATE  THE  AVERAGE  SIGNAL 

G 

690 

c 

G 

700 

IF  (CNT( I ) .EO.O)  GO  TO  60 

G 

710 

SIG(I)>SIG(I)/CNT(I) 

G 

720 

60 

CONTINUE 

G 

730 

IFIN*0 

G 

760 

LASTVR'IOOOOOOOO. 

G 

750 

c 

C 

760 

c 

PERMITTIVITY(PERH)  IS  FIRST  SET  TO  1  TO  CORRECT  RADIAL  FOR 
MULTIPATH  RADIATION  PERN  IS  INITIALLY  INCREMENTED  BY  KENT) 

FOR  EACH  CORRECTION  OF  RADIAL  UNTIL  STANDARD  DEVIATION 

G 

770 

c 

G 

c 

G 

c 

INCREASES.  PERM  IS  THEN  REDUCED  BY  1  AND  ENT  IS  SET  TO  .1  AND 
PROCESS  IS  REPEATED  TO  FIND  BEST  PERM. 

G 

800 

c 

G 

810 

c 

G 

820 

c 

START  OF  DO  LOOP  THAT  DETERMINES  PERMITTIVITY  OF  BEST  FIT 

G 

830 

c 

G 

865 

DO  120  I-1.20 

850 

DO  70  J«1.101 

G 

860 

875 

PMULT(J)«0 

G 

70 

OIF(J)»C 

G 

880 

R EF*CHPLX( PERN. -CONDPl 8000. 7FREQ) 

G 

890 

c 

G 

905 

212 


<>5 

100 

105 

110 

115 

120 

125 

130 

115 

140 

145 

150 

155 

160 

165 

170 

175 

180 

185 

190 


CALCULATE  MULTIPATH  OB  6AIH  FOR  EACH  RAN6E  CELL 


IF  (fiG( jiifo.c 1  CO  TO  ao 


RANGE«3.6»J*.4 

CALL  MULTPTH  ( RANGE. REF* HT* ACHT» tCKNP, OBGAIN* FREO) 

IF  ICKHP  NE  ZERO  DB  GAIN  CANNOT  BE  CALCULATEO 

IF  (ICKHP. NE.O)  60  TO  80 

PMULTI Jl-OBGAIN 

CONTINUE 

$UH«0 

SCNT-0 

SUBTRACT  HULTIPuTH  GAIN  FROM  SIGNAL 
00  90  J>ltlCl 

IF  (SIG(J).EO.O.OR.PNULT(Jt.EO.O.I  60  TO  90 
OIFij|>SIG(J)-PHULTU) 

SCNT»SCNT»1 


|UHaSUH«OIF(  Jl 
CONTINUE 

IF  LESS  THAN  10  RANGE  CELLS  OF  DATA  —  END  RUN  — 

IF  (SCNT.LT.IO)  GO  10  ISC 

AtfER>SUH/SCNT 

SUH2«C 

DO  100  J>1*101 
IF  (OIF( J I.EO.O)  GO  TO  100 
SUN2>SUN2«(01F(  J)-AVER)**2 
CONTINUE 

CALCULAIE  STANDARD  DEVIATION  OF  CORRECTED  DATA 
VAR.S0R'(SUH2/(SCNT*1) ) 

IF  LAST  CALCULATION*  60  TO  PRINT  FINAL  DATA 
IF  (IFIN.EO.il  GO  TO  130 

IF  VARIANCE  IS  INCREASING  REDUCE  PERMITTIVITY 


IF  (VAR.GT.LASTVR)  GO  TO  110 
LAStVR<VAR 


LAStVR<VAR 
PERM»PERM*ENT 
GO  TO  120 

PERMITTIVITY  INCREMENT  IS  ALREADY  .1  SET  END  INDICATOR 

IF  (ENT.E0..11  IFIN-1 
Pf RM.PERM-ENT 

GO  TO  160 

PRINT  FINAL  DATA 

WRITE  (6.23CI  PERN*CONO*REF*AVER*VAR 
WRITE  (6*240) 

PRINT  GRAPH  OF  SIGNAL 

CALL  PRPT  (SIGI 
WRITE  (6*2501 

PRINT  GRAPH  OF  MULTIPATH  EFFECT 

CALL  PRPT  (PHULT) 

WRITE  (6.26CI 

PRINT  GRAPH  OF  CORRECTED  SIGNAL 

CALL  PRPT  (OIF) 

GO  TO  10 

WRITE  (6*270)  TIHINUM) 

STOP 

WRITE  (6*280)  SCNT 
STOP 

WRITE  (6*290 
STOP 

FORMAT  (1H1./*1HTI 
FORMAT  (IX*6A1C) 

FORMAT  (1H1*30X*26HMULTIPATH  MATCHING  PROGR AN* ///* 20X*a AlO) 


FORMAT  (1H1*30X*26HNULTIPATH  Mi 
FORMAT  (2(4X.SI2I*3F1C.2*I10) 
FORMAT  (IH  */*14H  start  TIME  ■ 


1  TOTAL  SEC  ■  *F10.0*/*14H  STOP  TINE  «  *I4*3H  HR*I4*4H  NIN*I4. 
2H  SEC*19H  TOTAL  SEC  ■  *F10.0*/>26H  ANTENNA  HEIGHT  - 

3tF10.3f25H  AIRCRAFT  ALTITUDE  AGL  •  *F10.0*16H  FREQUENCY  -  *F10 
43*17H  SOURCE  NUMBER  -  .13) 

FORMAT  (///*2(1H  .100aH*)./)>/.lH  *40H  BEST  FIT  HAS  FOUND  WITH  P 
IRMITTIVITY  -  .F10.2.16H  CONDUCTIVITY  •  .F10.3»/*14H  CALCULATED 
216H  REFLECTIVITY  ■  .2F10.2. IH J. /* IH  .17H  AVERAGE  SIGNAL-  *Fl5.3*2 
3H  STANDARD  DEVIATION  ■  .F10.3*///*H1H  tlOOIlHRl* /) I 
format  (1H1.40X.43H  PLOT  OF  SIGNAL  DATA  USED  FOR  CURVE  FITTING./F 


«  A 

im 

61460 

mil 


>14. 3H  HR.I4.4H  HIN.I4.4H  SEC*18H 
OP  TINE  «  *i4*3H  HR.14.4H  nIn.I4.4 
.0./.26H  ANTENNA  HEIGHT  ■ 

•  .F10.0.16H  FREQUENCY  -  .FIO. 

IH  .40H  BEST  FIT  HAS  FOUND  WITH  PE 
VITV  •  .F10.3*/*14H  CALCULATED* 
*1H  .17H  AVERAGE  SIGNAL-  *Fl5.3*24 


F0RMAT^^1H1*40X*43H  PLOT  OF  SIGNAL  DATA  USED  FOR  CURVE  FITTING./F/ 
^fDRfiAT°(lHl*40X*38H  PLOT  OF  BEST  FITTING  MULTIPATH  EFFECT* ///* SX*2 


Gl$tC 

Gn^6?8 

G1620 

G1630 

t\m 

61660 

G1670 

GloaO 

M 

iti 

mil 

61790 

Giaoo 

61810 

61020 

Gia|3 

G1840 

C1830 

61060 

61090 

G190D 


FORMAT  (IHl.AAX.aSHFiaT  OF  MULTIFATH  CORRECTED  SICNAL  >/»IQXtS9H  I  6^10 
I  FOR  A  MrFECT  nf  THIS  CURVE  SHOULD  IE  A  CONSTANT  DIN  >»7/<SX»]H  civio 


is  ;i«!s;s4’'ssj{»‘!iSfUr.iu‘FSSiis  i!»icSi°«ii  c  iim 

xuDye  ...  _  _ _  - . .  ...  6i?feo 


IlH  »A«HIEST  FIT  STILL  NOT  FOUND  FOR  FERNITTIVITY  1  TO  20) 


SUBROUTINE  FRFT  (SIC) 

THIS  SUIROUTINE  CREATES  A  FRINTER  FLOT  OF  THE  ARRAY  SIC 

DIMENSION  SIG(101)»NTRX(6I»l01)»FRaS)»RN(ll)tRRN(I01) 
FR(i)— 5. 

^R(!?-{Rt{-!)«2.5 


LOOF  WHICH  SETS  THE  FOINTS  TO  IE  FLOTTED 


liiil!!:;'  '• " 


LOOF  TO  FRINT  ARRAY 
DO  40  I«l*12 

HR^TE  (6*10)  FR(14-I)«(HTRX(L»K)»K«1»101) 
DO  40  J-if4 

(6*100)  (HTRX(L*K)iK>l*101) 

CONTINUE 

FRINT  LAST  LINE  OF  ARRAY 
UR|||  FR(1)«(MTRX(1»K)*K>1*101) 

SET  VALUES  OF  RANGE  TO  IE  FRINTEO 


00*|i"Ji2*ll 

RN(|)>RN({-1)*4. 


RRN(|)*4 

DO  60  I*2fl01 

RRn7I)>RRN(!>1)«.4 

KK«S 

JJ-I 

WRITE  THE  ARRAYS  USED  FOR  FRINTER  FLOT 

HR  J?  \Ui\l  )  (SIG( J)*RRN(J)*J>JJ»KK) 
RK-KKRS 

i^'iKK^GT.IOl)  KK«101 
CONTINUE 

FORMAT  (9X*10(9(1H-)»1H*)»IH-), 

FORMAT  (2X«FS,1*2H  ♦flOIAlflN*) 

FORMAT  llXtlHitlOlAlelHl) 

fSrhat 

FORMAT  (IHliSlliH 


IH*) 

X*9HRAN6E  NH) 


H 
H 
H 
H 
H 
H 
H 
H 
H 

H  100 

H  tif 

H  PO 
H  160 
H  170 
H  180 
H  190 

H  ^ 

N  Iso 

n  k 

H  280 

d 

H  |l0 

H  Ho 
H  340 
H  3S0 
H  360 
H  370 

IIS 


H  $30 
H  $40 
H  S$0 


NAL  03  RANGE  NN ) ) 


SIIS 

H  $90 
H  600 
H  610 


H  6S6 
H  660 
H  670 
H  680 


SUIROUTINE  NULTFTH  (RNGEtREF* HI* ALTACt ICKNF»OICAIN>FREO) 

SUIROUTINE  TO  CALCULATE  NULTIFATH  EFFECTS  FOR  NATCH 
REF  AFFTC  TO-75-3 

COHFpX  NNfRRiRTfREF 

ICKMF>0 

NN-REF 

talt-altac 

FOL  «  1  FOR  vertical  AND  FOL  •  0  FOR  HORIZONTAL  FOLARIZATION 
VERTICAL  IS  ALWAYS  USED  FOR  THIS  FROGRAM 

rSncI*i  .149FRNGE 
CALCULATE  DISTl  AND  H2FRI 


i  IS 

I 

i  so 

{  IS 

i80 
90 

t?8 

I  its 

I  iso 


so 


CiLL  DIST  (R*NGE*HltTALT,OISTltICKnP>H2PRI) 

IF  (ICKNP.EO.l.OI  GO  TO  SO 
*RG2>Hl/(S2SO.O*OISTli 

FIND  ANGLE  OF  INCIO/REFLEC 

SI>(*TAN(*RG2i> 

IF  ANGLE  OF  REFLECTION  GT  .2  RADIAN  DO  NOT  CALCULATE  GAIN 

IF  (SI.GT..2)  GO  TO 
RT>CSORT(NN-CaS(SIM 

FIND  REFLEC  COEF  AND  ANGLE 

IF  (l.C-PQU  ictio.ao 
RR-(NN*SINiSn  -RT»/(NN*SIN( SI )*RT ) 

GO  TO  30 

RR'ISInTsII-RTI/ISINISD^RT) 

DIVERGENCE  IS  CALCULATED 

D>S0RT(l.aA((2.0*DISTl)*(RANGE-DISTlt PP2)/(RANGE*H2PRI) ) 
reflection  COEFFICIENT  IS  CALCULATED 
RnAG*SORT|R|AL(RRi**2vAInAG(RR)P*2) 


PHI>ATAN2( 

DIVl.O/D 


HAG(RR>>REAL<RRit 
CALCULATE  THETA 


H3»H1-OIST1**2/2.0 

THETA> (0.0001 3aB*H3*H2PRI*FRE0) /57.3/RANGE 
AN6LE3-THETA-PHI 

DB  GAIN  IS  CALCULATED 

GAIN«1.0»(OIVPRNAG)**2«2.0*DIV*RNAG*COS(ANGLE3) 

OBGAIN>IO*ALOGIO(6AIN> 

RETURN 

iCKNP-l 

END 


I  19U 
I  200 
210 
220 
230 

Itl 
260 
270 
260 
290 
300 
310 
1  320 
I  330 
I  3SC 
I  350 


{ 

I  SOU 


SIO 
S20 
I  S30 
I  sso 
I  S50 
I  S60 
I  470 
I  460 
I  490 


520 

530 

158 


I 
1 

1  560 
I  570 
I  560 
1  590 
I  600 
1  610- 


SUBROUTINE  DIST  (RAN6E>H1> TAL T»DI ST1» ICKHPrHlPR I ) 
SUBROUTINE  TO  CALCULATE  OlSTl  AND  H2PRI 
OIST1>C.O 

DH  IS  LINE  OF  SIGHT  RANGE 

OH.1.414*{SORT(H1)4SORT(TALT) I 
IF  (DH.LT. RANGE)  GO  TO  70 


IF  (OH.LT.ftANGE 

8lS^?»i!lt8?STI 


OISTl  IS  THE  DISTANCE  FRDN  ANTENNA  TO  REFLECTIVE  SURFACE 
(DISTl.GT. RANGE)  GO  TO  70 


H2PRI.HI*(RANGE-DI STD /OISTl 
H2SEC>( (RANGE-DlSTl)**2)/2.0 
H2fOT»H2PRl4H2$£C 
IF  (TALT.GE.H2T0T)  GO  TO  20 
CONTINUE 


0IST1«0. ODDIS 


hIprI-hipIranEetOISTiddist^ 


_  ^-HIPIRANGE-OISTI 
H2SEC>(lRANCE-0IStl)P 
H2T0T-H2PRI4H2SEC 
IF  (TALT.GE.H2T0T)  GO  TO  40 
CONTINUE 

_ERl*Hip(ilANGi-ElSTl  )/0IST 
H2SEC«((RANGE-0ISTl)PP2)/2. 
H2f0T«H2PRI*H2SEC 
IF  (TALT.GE.hHoT) 

CONTINUE 
return 
ICKNP«1.0 
END 


OIS 

00 

01 

H2 


GO  TO  60 


10 

20 

30 

58 

^78 

80 

90 

t$8 

ii8 

140 

150 

160 

170 

160 

190 

200 

Hi 

230 

240 

250 

260 

270 

260 

290 

m 

320 

118 

350 

m 

380 

390 

400- 


OVERLAYtnULTf  3>0) 

OVERLAY  NULT  PREDICTS  THE  NULTIPATH  EFFECTS  ON  THE  DATA 

A  PRINTER  PLOT  OF  THE  GAIN  VS.  RANGE  IS  PRINTED  ALONG  WITH 
THE  CORRESPONDING  DATA  FOR  EACH  CELL.  REF  AFFTC  ~~  ~ 

PROGRAM  NULT 


TO-75-3 


COMPLEX  NN>RR>RT 

DINENSIQN  AMTRXilOOIfNTRXtSOtlOO 
DIMENSION  XNfRXaiT»ALt(8)tRNCTI 


00)>ELAN6(100)tREFO(100)>OB(100) 


10 
20 
30 
.40 
50 

^78 
60 
90 
K  100 
K  110 


III 


15 

20 

25 

30 

35 

♦  0 

^5 

50 

55 

60 

65 

70 

75 

so 

85 

90 

95 

100 

105 

110 


PRINT  FORNAT  WILL  BE  8  LINES  PER  INCH 


WRITE  <6>120) 

THE  ALTITUDES  ARE  READ  IN 

READ  (5«140)  ALT 
DO  5  I>1«8 


ADAP  CARO  NO.  6 


READ  (5«] 
DO  5  I>i| 
If lALTlf 1 
NALT«I 
CONTINUE 


.EQ.Ot  CO  TO  10 


THE  BEGINNING  RANGE  AND  END  RANGE*  ANTENNA  HEIGHT  IN  FEET* 

THE  FREOUENCT  in  HHZ*  POLARIZATION*  PERNITTI WITT*  AND 
CONDUCTIVITY  ARE  READ  IN 

READ  AOAP  CARO  NO.  9 

READ  (5*150)  RANCE1*RANGE2*H1*FRE0*P0L*E*S 

IF  (EOFIS)  .NE.O.)  GO  TO  IID 

XX-0.0 

YV«0.0 

A00>(RAN6E2-RANGEl)/10.O 

LOAD  THE  RANGE  ARRAY  FOR  PLOTTING  OH  THE  X-AXIS  Of  THE’  PRINTER 
PLOT 

DO  20  KK>1*11 

XNTRX(KK|aRANGElPXX 

XX>XX«AOO 

FIND  INDEX  Of  REFRACTION  _ 

NNaCnPLX(E>-S*lBOOO./'  REO) 

PRODUCE  PLOT  AND  PriNT  DATA  FOR  EACH  ALTITUDE 

DO  100  I-1*NALT 
ONINalOOO. 

DNAX'O. 

RANGEaRANGEl 

RANGE>1.1A9*RANGE 

talt»altTi) 

DO  BO  JJ«1*100 

FIND  CORRECT  OIST  1 


RANCE*H1*TALT»01ST1»CHECK*H2PR1*0NIN*0NAX) 

ANhALT(I)-Hl)/(RANGE*52S0.))*57.6 


ALL  ODIST  (RANCE*H1>TALT» 

N6(JJ)«RANGE/1.199 

LANG(JJ)aATAN((ALt(Il-Hl) 

EF0(JJ)a0ISTl/i,U9 

F  (CHECK. EQ. 1.0)  GO  TO  60 

RG2«H1/(S280.0*OIST1) 


FIND  ANGLE  OF  INCIO/REFLEC 

SI-tATAN(iARGf)1 

Rt«CSQRT(NN-CaS($I)P*2) 

FIND  REFLEC  COEF  AND  ANGLE 

IF  (l.O-POL)  30*30*60 

REFLEC  COEF  FOR  VERTICAL  POL _ _ _ 

RR-(NN*SIN(SI)-RT)/(NN*SIN(SI)*RT) 

GO  TO  50 

REFLEC  COEF  FOR  HORIZONTAL  POL 
RR>(SIN(SI)-RT)/(SIN(SI>'»RT) 

0«S0RT(1.0*((2.0P0IStl)*(RANGE-0ISTl)**2)/(RANCEPH2PRI) ) 
RttAGaSQRT(REAL(RR)**2«AINAG(RR)PP2)  — 

PHI>ATAN2(A(NAG(RR)*REAL(RR)) 

IF  (SI.GT..2)  GO  TO  60 
DIV'1.0/0 

CALCULATE  THETA 

H3-Hl-DISTl**2/2.0 

THETA-(0.0001385*H3*H2PRI*FREC)/S7.3/RANGE 

ANGLE3aTHETA-PHI 

CALCULATE  CAIN  DUE  TO  NULTIPATH  EFFECTS  AND  LOAD  IN  THE  ARRAY 

anTrx  for  plotting 

GA INa 1. 0* ( DiyPRNAG) **2t2. 0*01 VPRNAG*COS( ANGLE  3) 
DSGAtNaf0*ALCG10(GAINT 
ITaOBGAIN 


DSGAINalOPALCG 

OBtJJTaOBGAIN 

ANTRX( JJ)at  aIN 

GO  TO  70 
ANTRX(JJ)aO 
0B( JJ)a6. 


ANTRX( JJ 
0B( JJ)a6 
CONTINUE 
RANGEaRAI 
CONTINUE 


NG|aRANGEP(RANCE2-RANGEl) Pi. 169/100 

NTIHUE  _  _ 

CALL  PLOTTING  ROUTINE  THAT  PRODUCES  THE  PRINTER  PLOT 


K  120 

j  in 

K  150 

s  in 

K  180 
K  190 
K  200 
K  ^10 
K  220 

5  m 

K  250 


^  in 

K  600 
K  610 
K  620 
K  630 
K  660 
K  ^50 
K  660 
K  670 
K  6S0 
K  690 
K  500 

5  m 


K  580 
<  590 
K  600 
K  610 
K  620 
K  630 
K  660 
K  650 
K  660 
K  670 


K 

K 

4 

K 

730 

K 

760 

K 

750 

K 

760 

1C  T70 

K 

780 

K 

790 

K 

800 

K 

810 

K 

820 

K 

830 

K 

860 

850 

K 

860 

K 

K 

tn 

K 

890 

K 

900 

K 

910 

K 

920 

K 

930 

K 

960 

K 

950 

K 

960 

K 

970 

K 

980 

K  990 
KIOOO 
XIOIO 

kI020 


KllOO 

KlllO 


2K 


US 

120 

125 

130 

135 

UO 

1 

5 

10 

15 

20 

25 

30 

35 

^0 

1 

5 

10 

15 


CALL  GNPLOT  ( AnTRX> NT R X> XNT RX , T AL T ) 

PRINT  160>  R*NGEl>R*NGE2>Hl,FRE0>PaL> E>$ 

0NIN*DMIN/1.1A9 

0NAX«DNAX/1.149 

PRINT  170>  NNfONINtONAX 

PRINT  180 

DQ  90  N>1>10 

Kl»l*10*fM-ll 

K2«K1*9 

PRINT  190«  (RNG(KI,K«K1,K2) , t DB( JJ) J>Kl>K2)f (ELANG( J)# J>K1,K2)> ( 
lREFD(L)«L*Kl«K2i 
CONTINUE 
CONTINUE 
GO  TO  10 
STOP 

FORMAT  (IHTI 
FORMAT  (12) 

FORMAT  (15(F10.3)) 


FORMAT  (7F10.3.I2) 
FORMAT  (lXW>lX>20e 


f-UKHAI  (  IXf  /  •  IXfZOHUl  )  I  ANLt  INItKVAL  NH  •  H  A  >  lOMAN  I  C  NN  A  Ht  1  bn  I  r  I  •  >  3 
iXf 13HFRE0UENCY  MHZ, 2X, 1 3HP0L ( 0«H,  1 > V ) . 2X , 12HPERN1 TT I V I T Y< 2 X« 12HC 0 
2NOUCTIVITY>/>1X>F6.0«  IX •  2rtT0>  1 X.  F6 . 0.  13 X>  F  5 . 1>  15X , F d .  2 1  5Xt  F 8 . 2>  12  X 
3«F&.2>AX,F7.5) 

PORMAT  (IXf/f 3X«33HCALCULATED  INDEX  OF  REFRACTION  S DUAR E 0 a 4X« 60H0  I 
ISTANCE  FR.  antenna  TO  REFLECTION  PT.  NN .  MINIMUM  MAXIMUHt /• 13Xa F 
2).4a1x. IH  «1XaF9.4»1Xa1HJa52X>F8.2>  3XAFe.2) 
rORMAT  ( IHIaIOFIO.S) 


FORMAT  (IH  aISHRANCE  (NM) 


1  >10F11.2a/aIX,18HELEV.  ANG  (OEGI 

2NM  )  aIOFU.Za/) 

END 


a1^-11>2a/a1Xa18HGAIN  (DB) 

a10F11.2a/a1X,13HREF.  DIST.  ( 


K112C 

K1130 

K1140 

<lliQ 

^i!?8 

<1180 

<1190 

<I|00 

K1210 

<1220 

<1250 

<1260 

<1270 

<1280 

<1290 

<1300 

<1310 

<1320 

mu 

<1350 

<1360 

<1370 

<1380 

<1390 

<1400 

<1430- 


SDBROUTINE  ODIST  (R ANGEaHIa TAD > 0 1 STl aCHEC<aH2PR1 a OMIN, DMAX ) 

SUaROUTINE  TO  CALCULATE  OISTl  AND  H2PRI 

CHEC<»C.0 

0IST1*C.0 

DH  IS  LINE  OF  SIGHT  RANGE 

OH»1.414*(SQRT(H1)*SORT{TALT) ) 

IF  (OH. LT. RANGE)  GO  TO  70 
00  10  I<1a200 

OISTl  IS  THE  distance  FROM  THE  ANTENNA  TO  THE  REFLECTING  SURFACE 

0I5T1-C.1*0IST1 

IF  (OISTl. GT. RANGE)  GO  TO  70 

H2 PR I.hl* ( R ANGe-0 1ST  1) /OISTl 

H2SEC*( (RANGE-OISTI) •*21/2.0 

H2T0T«H2PRI*H2S|C 

If  ( IALT.GE.H2T0T)  GO  TO  20 

CONTINUE 

OISTl-OISTl-.l 

00  30  I«l{100 

OISTl-O.Ol*DISTl 

H2 PH *H1* ( RANGE-DISTI ) /OISTl 

H2SEC-( ( RANGE-OISTI) ••2) /2.0 

H2TOT.H2PRUH2SEC 

IF  (TALT.GE.H2T0TI  GO  TO  40 

CONTINUE 

OISTl«OISTl-.Ol 

DO  50  I-1a100 

OISTi*0.001*DISTl 

h2 PR r ‘HI* (RANGE-DISTI) /OISTl 

H2SEC*( (RANGE-OISTI ) ♦♦2) /2.0 

H2T0T*h2PRI*H2SEC 

IF  ( TALT.GE.HITOT)  GO  TO  60 

CONTINUE 

CONTINUE 

OMIN.AMINKOMINaOISTI) 

OMAX«AMAX1(OMAXaOIST1 ) 

RETURN 
CHEC<  -l.O 
END 


L  20 
L  30 


L  80 
L  9C 
L  lOO 
I  110 
L  120 


L  150 
L  160 


L  21C 
L  220 
L  230 
L  240 
L  250 
L  260 
L  270 


L  300 
L  310 
L  320 
L  330 
L  340 
L  350 
L  360 
L  370 
L  380 
L  390 
L  400 
L  410 
L  420 
L  430 
L  440- 


SUBROUTINe  GNPLOT  ( AMTRXa MTR X a XMTR X a T AL T ) 

SUBROUTINE  TO  PRINT  A  PLOT  Of  THE  GAIN  DUE  TO  MULTIPATH  EFFECT 
VS.  THE  RANGE 

INTEGER  A(5) 

DIMENSION  AMTRXIIOOIaMTRXISOa 100>a  XMTRX(ll) 

Al»7.0 

A2*6.0 

A3>4.B 

A4«3.0 

A5«0. 

DO  10  I*1a5 
A(  I)»I 
CONTINUE 
LL*10*LL-.10 
IF  (LL.EQ.O)  LL-1 
DO  20  M«1a50 


M  10 
M  20 
M  30 
M  40 
n  50 
M  60 
M  70 
M  80 
M  90 
M  100 
M  110 

3  il8 

M  140 
M  150 
M  160 
M  170 
M  160 


217 


00  to  N«l*100 

LOAO  THE  PLOTTING  ARRAV 


IF  (I.lo.O.t  CO  TO  30 
HTRXn>JI>lH. 

^RlSPsOf  TALT 
I-l»50 


ANTRX(J)*S0.0 
•  5 


00  AO 
K-51-I 
IF  IK, 


IFo(K.N|j^S^.AN^.K.NEjAgjANO.K.NE.30.ANO.K.NE 

WQ.SO.)  FRINT  6Q«  A(  N  If  (HTRXIK*  J  t  >  J>lf  100>>  A1 
O.AO.j  FRINt  6Qf  AlNW jNtRxjKf jjf J>  # l55|f At 

o.ao.i  rrInt  60t  AiniflnTRxjKfjjfj-  fioQifAi 

O.to.l  FRIRT  eOf  Alnlf (NTRXIKf J|»J>  flOOifAA 
O.IO.I  FRINT  60f  AINIftNTRXIKfilf J>ifl00lf AS 
CONTINUE 
FRINT  86 

FRINT  90#  XNTRX  - - 

RETURN 

F0RNAT(lHlfS0X.31HF0UER  GAIN  FLOT  AT  ALTITUOE  OF.f /fAOX.FlO.Sf  lOH 
J  FT.^A6L.t//f3Xf lOHFOWER  GAINfliOXfAHOIGAlNt/f ITXf iHlf lOliHAf 9(Ih 

iiliAiT^IlH  f  lOXfIltSXf  lH«(100AlflH*f3XfF3.1l 
18!!!}  - - 


AN0.K.NE.30.AN0.K.NE.20.AN0.K.NE.10)  FRINT 


6df  A(Nlt(NtRXIKfilf J>lfIO 


OVERLAYiRONERGEfAtO) 

OVERLAY  THAT  NERGES  OF  TO  THREE  RAOAR  TAFES  TOCETtTfR - 

FROGRAH  RONERGE 
OINENSION  I2l2)fT2(2itTA(6l 

REAO  AOAF  CARO  NO.  2 

REAQ  (SfSO)  NUN 

write  (OfOOl  NUN  - 

NN«11 

00  AO  I'XfNUN 

REAO  (NN)  IlfTlf I2f TtflS.TA 


00  AO  I'XfNUN 

REAO  (NN)  nfTlf I2f T2fI3tT4 

ININ-fifli 

GO  TO  10 

NNaNNAl 

CONTINUE 

WRITE  (6.701  IHRfININ 


FORNAT  (110) 

FORMAT  UH  fI3t29H  RAOAR  TAFfS  WILL  OE  «RGC^-7 - 

FORMAT  (2AH  LAST  TINE  MERGED  NAS  .lA.SH  HRflA.AH  MIN) 


N  80 


IN 


APPENDIX  G 


ADAP  SOURCE  LISTING 


PROGRAM  ARPnP( INPUTf0UTPUT*TAPE5aINPUT»TAP£6*0UTPUT*TAPE10»TAPEll> 
1TAPE12) 

THIS  PROGRAM  CALCULATES  A/C  ANTENNA  RADIATION  PATTERNS  (IN  ERP  OBU 


?i5*prii“tiEr" 

IKX-i 

ROSTJ-O.O 

iiwiNO  10 

REWIND  11 
NUND«0 

WRITE  (6>290l 


READ  INPUT  CAROS 


READ  HEADING  RECORD  (ARPNP  CARO  t  II 


READ  (5.9201  IHEAO 
WRITE  (6.9201  IHEAO 


READ  ARPNP  CARO  •  2 


READ  (5.3001  IR AOAR. IPCAL. I  HULT.NTTPP.NUNSR. I  STOAT. HUNRNG. ALTAC 
WRITE  (6.3101  NTYPP.NUHSR.IRADAR.NUNRNG. ALTAC. IPCAL. ISTDAT.INULT 


READ  ARPNP  CARO  •  3 


IF  NULTIPATri  OPTION  USED  —READ  ARPNP  CARO  I  9 

READ  (S.320I  ( ANTHK 1 1 .PERN( I  I.CONOl I It I-l.NUNSK I 
DO  110  I'l.NUNSR 


CALCULATES  REFLECTIVITY 


REF(II«CHPLX(PERN(Ilt-CgND(II*l«000./FRE0(III 
WRITE  (6.3901  ( It  ANTHTT  I  I.PERH(lI.CONO(  1 1  »REF 
IF  ( IRADAR.NE.IHNI  GO  TO  190 


(II.I-l.NUNSRI 


IF  NO  PAOAR  TAPE  —  READ  RANGE  AND  CORRESPONDING  TINES.  ARPNP 
CARO  I  5 

READ  (5.3501  ( RDRNGI 1 . I I . I0ATA(6* ( l-n An . IDATA (6*( 1*11 t2 I. IDATA( 6 
l*n-llAll.t6AtA(6*(l-i|t9i.I0AtA(6*(I*llAM.IDATA(6*(  I-lU6I.I>l.N 
2UNRNGI 

CALCULATES  TOTAL  SEC  OF  RANGE  TINES 
DO  130  I>1.NUNRNG 

RORNG(2.(I«IOATA(6*(I-1I»1|P3600.*IDATA(6P(I-1I aZ|*60.aI0ATA(6*(I- 
111A3I 

R0RNG(3. II«I0ATA(6*( I-llA9l*3600.♦IDATA(6*II-llA5l*60.AIDATA(6*(I- 
lll  A6I 

WRITE  (6.3601  ( IOATA( 6A( I-l  I AlleI09TA(6*( I-ll >2 I. IDATA( 6* ( I-l I a3I > 
lI0ATA(6*II-ll«9I.IDAtA(6*(t-ll«9).l0ATA(6*(i-lTA6I.I>l.NUNRNGI 

CALCULATES  LOSSES  TO  BE  USED  IN  ERP  CALCULATION 

DO  150  |>l.NyNSR 

OELOB(II>SLO$S( II-RGAIN(I|a37.Ba20*AL0G10(FRE0( II I 
CONTINUE 

DO  160  i«l.NUHSR 
I0ATA6(!|>0 

READ  IN  THE  HEAOING^RECORO^ON^SIGNAL  JJJJJJJJJ**J*J*J*** 

BUFFER  IN  (10.11  ( IOATA( II . IOATA( 1)1 
IF  (UNlfdOil  180.250. IBC 
CONTINUE 
IK-II 

IltlWROIll 

tK-II 

tl-AN0(SH:FT(II.-13I.TBI 

IF  II  EQUALS  9  —  HEADING  REC0R3  FOUND 
IF  (II.E0.9I  60  TO  190 

IF  HEADING  RECORD  NOT  FOUND  BY  SECOND  TRY  —  END  RUN 

IF  (ICK.EO.ll  GO  TO  270 

ICK»1 

REWIND  10 
60  TO  170 

READ  IN  THE  TINE  CALIBRATION  RECORD  ••••♦♦♦t*»t*****9******** 


105 


110 


115 


120 


125 


130 


135 


1%0 


l<i5 


150 


155 


160 


165 


170 


175 


180 


185 


190 

200 


21C 


22C 

C 

c 

c 

230 


C 

c 

c 

24C 


250 
260 
2  70 
2  80 

285 

C 

2  90 
300 
310 


320 

330 

3*0 


3  50 
360 
370 
38C 
390 

tn 

*20 

*30 


SET  TIME  TAG 


BUFFER  IN  (10.lt  ( I  DAT*(  1 1 .  ID  AT  A(  12  I  I 
IF  (UNITdOlt  200.250.200 
CONTINUE 
Il-IWROh) 

IK«II 

n«AND(SHIFT<II>-13t.7BI 

IF  II  EQUALS  6--  TINE  CAL  RECORD  FOUND 

IF  ( II.NE  .61  GO  TO  260 

ITTAG>IUR0(*3i 

WRITE  (6.370) 

DO  220  I'l.NUnSR 

READ  DATA  CAL  RECORD  FOR  EACH  SOURCE  •**•*♦♦*♦♦•♦•♦♦♦♦♦♦♦**♦* 

BUFFER  IN  (10.11  (I0ATA(1).I0ATA(12)) 

IF  (UNITdOil  210.250.210 

CONTINUE 

IK>IWRDd  I 

KK"I 

II-IK 

1 1 -AND! SHIFT! 11.-131.78) 

IF  II  EQUALS  6  SET  DATA  TAG 

IF  (II.NE. 6)  GO  TO  280 
I0ATAGd)>IWR0(*3) 

CONVERT  80  EBCDIC  CHARACTERS  TO  CDC  CODE 

CALL  lOFUN  dCHR) 

WRITE  (6,380  I.ICHR 
CONTINUE 

CALL  SUBROUTINE  CAL  TO  CREATE  SIGNAL  CALIBRATION  ARRAY  ***♦♦•* 

CALL  CAL  dPCAL) 

REWIND  10 

NUND«0 

0AVS»0 

DO  2*0  II-l.NTYPP 


AlOlO 

\m 

AlOSC 
A 1060 
AlOTO 
AtOSO 


A1300 

AlSlO 


CALL  SUBROUTINE  READD  TO  START  PROCESSING  DATA  ♦♦♦*♦♦*♦♦♦♦* 

t****t*i*t*******t*  ********  »'»*»*»»*»»»**»****»*************  ******* 

;all  readd 


Continue 

END  FILE  12 
REWIND  10 
REWIND  11 
REWIND  12 
STOP 

WRITC(6.*30I 

STOP 

WRITE  (6.390) 

STOP 

WRITE  (6,*0CI 
STOP 

WRITE  (6, *10)  KK.IK 
CALL  settime 
00  285  J«1,NUMSR 
lOATAGd)  >1 
GO  TO  230 

FORMAT  ( IHl, 35X.*2HAIRCRAFT  ANTENNA  RADIATION  PATTERN  PROGRAM) 
FORMAT  (3Al.t7.3ItO.F10.3l 

FORMAT  dHO,**HTHE  NUMBER  OF  DIFFERENT  FLIGHT  PATTERNS  ARE  ,I*./.l 
IH  ,22HNUMBER  OF  SOURCES  ARE  . I *. lOX . 8 HR  AD AR  •  .Ai.22H  NUMBER  OF 
2RANGFS  ■  ,1*,11H  AC  ALT  •  .Fig,2./.23H  PRE  AND  POST  CALS  >  .Al 
3,13H  START  DAY  •  .I*./.31H  MULTIPATH  SIGNAL  ELIMINATED  «  .Al) 
FORMAT  (6F10.3I 

format  (1H  .1*HS0URCE  number  .I*.12H  FREQUENCY  ■. FS . 3. . 23HMHZ  SY 

ISTEM  LOSSES  •  .F5.2.20H  DB  RECEIVER  GAIN  •  .FS.Z.SH  DBI 
FORMAT  d5H  SOURCE  NUMBER  ,I2.18H  ANTENNA  HEIGHT  •  .F8.3.16H  PERMI 
ITTIVITY  •  .F8.3.16H  CONDUCTIVITY  «  .F8.3,/  ,31H  THE  CALCULATED  R 

2EFLECT1VITY  •  .2F10.3.2H  J) 

FORMAT  (*(F3.0.1X.3I2.*X.3I2) ) 

FORMAT  d9H  TIME  OVER  SITE  •  .312. 9H  TIME  AT  .F6.1.6H  NM  ■  .312) 

FORMAT  (IHO. 3HTAC.*X,6HSO0rC£ ) 

FORMAT  ilH  .13.*X.8A10) 

FORMAT  dX.25HTIHE  CAL  RECORD  NOT  FOUND) 

FORMAT  dX.2*HHEAD|NG  RECORD  NOT  FOUND) 

FORMAT  jlH  .6HS0URCE.I*.2DHCAL  RECORD  NOT  FOUND. 08) 

FORMAT  (SAlOt 

FORMAT!"  END  OF  FILE  ENCOUNTERED  ON  SIGNAL  TAPE  — TAPEIO—  NO  DATA 
IRECORO  FOUND") 

END 


Ai*5g 

Al*60 

Al*70 


Al6l0 

A16*0 
Ai650 
A1660 
A 1670 
A1680 
Ai690 
Ai700 
A1710 
A 1720 
A1730 
A17*0 
A1750 
Ai760 
A1770 
A17B0 
A1790 
AlSOO 
AIBIO 
AlBZg 
A1830 


1 

5 

10 


SUBROUTINE  READD 

•  THIS  SUIROUIINE  CREATES  THE  FINAL  OUTPUT  DATA  IT  ••••••*** 

*-  .CAIH***  RADAR  ANg  signal,  •?•!*!•!.••*•••• 

TYPE .NUMS*. 
ogitCAi 

P.Dl  iniSOB.XI.HUMUl 
....  .  .AITAC.OATS.OATR.IS 
FREOlZOtf IMULT 


COMMON  ITtPE.NUMSR.  IdaIa  dU).  ITTAG.igAT 

(ins 

RORN^|j^2^).AlTAC.6AT$.OATR.inOAl.Rastf.ANTHT(2dl,MFI20). 


COMP!  fX 

oTmef 


RFF 


ION  IHRll90).NlNll90).ISECl(90)#IHR2l9Ul*NtN2(90I.IS{C2(90). 


2«1 


15 


20 


2» 


30 


35 


«0 


55 


50 


55 


60 


65 


70 


75 


80 


85 


90 


95 


100 


105 


no 


80 


90 


1  V0*T(5821»SI6(20)*SNAX(2OlfSNINI'*0> 

•  . . . . . 

READ 


**•* 
R 


AD  IN  THE  HEADING  TYRE*  THE  NUNtER  OF  HfADI 
JHJ^Rg^^^AN6jE*^JH|^^lTJH^ANGtE~ARRIIR^CAR0 


NGS 

. . .  g  S  ******•••«•* 


UR 

If* 

•  •I 

READ  IN  THE  HEADING  AND  THE  C0RRESR0NDIN6  TINE  INTERVALS  ••••*** 
*•♦*««•******  ARPNR  CARO  8  9  ««*•**•***•*•*•****•«**«****•*** 

*ttt**********«A**R***«6***«**«***R********R******R*****«A*««**R** 

.  IHE0(II(IHRl(I)>NINl(II*irFCl<I>*IHR2II>(HIN2(I)*IS 
NUNHOT _ - 

~~'t  ((IHEO(I)f IHR1(I>>H1N1(II»1SEC1II)» 
■IfNUNHOI 

I  l(lHCOin*IHRl(I)tNINl(Iif  ISECKin* 


READ  15*2001  ((I 
1K2TI)^-!*NUNH 
!F(lHOTrR,EO,ll 
lfHRi(I)*NiN2ni* 


IFtfHOtYR.H) 

1I«1*NUNHD) 

.  TRlN-a. 


10 

c 

? 

20 


MI-1*NUNH0 

nain.i 


(n*60»iSECiin 

177 


ilN2(II>rHR2in*3600»HIN2in*604lSEC2(I> 

*  i!J"S;  ••*•*****«*#«*««*«•*•*  *  •  «  *  •* 

FIND  DATA  FOR  ALL  TINE  INTERVALS 

DO  20  NI-l.NUNSR 
SNAxThII— 9999. 

SNIN(HI)-9999. 

NN-1 

DO  30  J-1.582 
VDAfiJI-O. 

CONTIMUE 

IF  DISCRETE  DATA  POINTS  SET  NN  TO  NUN8ER  OF  POINTS  INUNHO) 

IF  (IHDTYP.NE.3>  go  to  50 

IF  CONTINUOUS  DATA  SANPLES  SET  NN  TO  NUN8ER  OF  .5  SEC  INTERVLS 


iT-|H|T*|6S8t{NiT*6St{ilT 

NN«ISP-STT*2.0 


50 


55 

? 

C 


50 

60 

65 


TIH-ST-.5 


SETS  NN  POINTS  OF  DATA  FOR  THE  INTERVAL  ♦♦*••*•*♦♦•*•♦•**•*♦•* 

*55*5«P*PP**«*«*5**5«555*****5***5********P**************P*****P* 

IF  (IH0TYP.EQ.2I  NN-NUHHQ 
IF  (fHOTYP.NE.n  GO  TO  55 

NH-<t{N2(NyNH5r-TTH . 

TIH«TIF 


Ni|NyHHO)-TlNl(in«2.0 


..  lHHn-.5 
??  itnO^Y^lEcl.E)  GO 
CONTINOUS  DATA 
-nftM.iN*,.? 


TO  90 


TO  70 

,AN0.TIN.LG*TlN2tin  GO  TO  70 

60  I-1*NUHH0 

TTlnttE.TIN2<II»  GO  TO  65 
CgNJIHU 


IF  (THQTYP.NE.lT  GO 

U  Irilult^^lNKD.A 

CONTINyE^  . . 


I^??T 


««««*•*«««««*••«««««««•««»«•««••««*«*«*»*««**•«**««  •  ** 

GET  HALF  OF  A  SECONO  OF  AVERAGED  DATA 

**«*g**«««***«*«*****««**A«*g««*ii***«***** ****•*••*• •*•*•*•****# 

CALL  RADAR  (  ASP*DEP*TIN*  TRIN*  RNG*  IROLL*  IPTCH.RNGXY*  IRCK I 
IF  (IRCK.EO.lV  60  TO  160 

Call  signal  (tih*sigi 

DO  80  Ln-lfNUNSR 

IF  (SIG(LNI.EO.O)  GO  TO  80 

iCKHP-0 

D8GAIN-0. 

CALL  NULTPTH  IF  OPTION  IS  SELECTED 

(RNGXY*LN*ICKNP*D8GAIN) 

IF  (ickhp!eo,ii  go  to’o?* 

CALCULATE  SIGNAL  IN  06U 

ESI6(LH)>VCALISIG(LH)>LH)«OEL08(LN|520*ALOG10(RNGI 
G(LH)«SIgTlH)-30.>0IGAIN 

miisiiiiiauihai'.iiEitRii 

NTINUE 
GO  JO  130 

DISCRETE  DATA  SANPLES  GET  AVERAGED  DATA  FRON  T1N-2.0  TO  TIN 

H52-IHE0(J) 

‘  K*1*NUNSR 


242 


115 

120 

125 

130 

135 

140 

155 

150 

155 

160 

165 

170 

175 

lao 

1S5 

190 

195 

200 

205 


TIN2(K)-0. 

CONTINUE 

TIN«TINl(JI 

CALL  RADAR  <  ASP«DER>Tin»  TRIM*  RN(*  IROLL*  IPTCH*  RNCXTt  IRCK I 
TIN-TIN-2.5 
00  lie  JJ-1.5 
TIN-TIN-. 5 

CALL  SIGNAL  (TIN.SIGi 

?0  110  K-l.NUHSR 
IN2Tk)-t!n2(KI-SIG(K) 

ONTINUE 


DATA  AVERAGED  FOR  2  SFCONDS 


00  120  K-l.NUNSR 

IF  (TIN2(K).E0,0.0)  GO  TO  120 

TIN2(K)-TIR2lKl/5.0 


ICKHP-0 

08GAIN-0. 


CALL  NULTPTH  if  OPTION  SELECTED 

IF  (INULT.EO.IHY)  call  NULTPTH  (RNGXYt K» ICKNPf OBGAIN) 
IFilCKNP.EO.lf  SIGlKi-O. 

IF  (ICKNP.EQ.l)  GO  TO  120 

CALCULATE  SIGNAL  IN  OBW 

SI6(K|-VCAL(TIN2(K),K|-OELOB(K|-20*ALaG10(RN6l 

SIG(K)-SIG(K)-30.-0BGAiN 

SNAXIKl-ANAXlISNAXfKT.SIGCK)) 

|NIN(KI-ANINl(SNIN(Ki»SIG(K>l 

CONTINUE 


CREATE  THE  OUTPUT  ARRAY  VDAT 


DO  150  K-l.NUNSR 
IF  (SIGIKT.EO.O.)  go  to  150 

PLACE  NON-ZERO  DATA  IN  OUTPUT  ARRAY 

V0AT(HNI-K 

VOAT|NN-il-SIG(K) 

V0AT(NN-2j«DEP 

VDAT(NN*3)-ASP 

V0AflNH*5»-RNGXY 

V0AT(NN-5I>TIN 

NN-NN-6 

IF  (J.Ea.NNI  VOATIAaZl-l 
IF  (NN.LT.AOOI  GO  TO  150 
V0AT(5ai»-585  , 

WRITE  (12)  VOAT 
WRITE  (6.2201 
WRITE  (6.230)  VOAT 
NN-1 

DO  150  L-1.582 
VDAT(L)-0. 

CONTINUE 

CONTINUE 

IF  (NN.EQ.l)  GO  TO  170 
VpAT(5ei)-HH-l 
VDAT(5e2i-l 
WRITE  (12)  VOAT 
WRITE  (6.220) 

WRITE  (6.230)  VDAT 

WRITE  (6.250)  ( NI.SHIN(NI )* SNAX(NI ).NI-1.NUNSR) 
CONTINUE 

B  c  Ti  I  Dm 

WRIT£(6,250)  IHDTYP 
STOP 


WRITE(6»260)I.IHR1{ D.NINK I ) . ISECl ( I )* IHRl ( I-l ) . NINl (I-l ). 
IISECI(I-I) 

STOP 

FORNAT  (5I10>2(5X.3I2)) 

FORNAT  (1H1.20HTHE  HEADING  TYPE  IS  .12. 


FORNAT  (1H1.20HTHE  HEADING  TYPE  IS  .12. 

A  "  (l-FLYBY.RAOIAL.  2-CLO  VEPLEAF.  KOYGON.  3- ORBITS  ) ". /.  IX. 

B  "The  nunber  of  hex  \ngs  arE". 

1  1X.I2./.1X.18HTHE  KOLL  ANGLE  IS  .IX. 13./.  "  THE  PITCH  ANGLE  IS 

2". IX. 12) 

FORNATdX."  START  TINE  IS  ".12."  HRS  ".12."  NINS  ".12."  SECS". 

1"  STOP  TINE  IS  ".12."  HRS  ".12."  NlNS  ".12."  SECS") 


2CS.10X.2HT0.10X.I2.1X.3HHRS.iX.l2.1X.5HNINS»lX.I2.1X.5HSiCS) 
FORHATIiNO.^THE  NAGNETIC  HEAPING  IS  "tlA.  /."  THE  CORRESPONDING  T 
IINE  IS  ".12."  HRS  ".12.*  NINS  ".12."  SECS") 

FORNAT  (1H1.51HTHE  DATA  ARRAY  TO  BE  USED  BY  ADAP  IS  - ./.2X.2(60 

IH  SOURCE  SIGNAL  pEP  ANGLE  ASP  ANGLE  RANGE  TINE  )./ 
2.2X.?(13X.5H(08W).5X.5H(DEC).6X.5H(DE6T.6X>5H(NN).1X.11H(T0TAL  SEC 

^Fi^HAl  (IH  .12F10.3) 

FORNAT  (IH  .3(50H  SOURCE  NUNBER  NIN  SIGNAL  NAX  SIGNAL )« /.5X. 7 ( 
13(5X.I3.7X.F10.3.2X.F10.3.5X)./)) 

FpPNAt("  illegal  heading  TYPE  "*110) 

F0(.NAT("  entry  ".13."  IN  THE  HEADING  ARRAY  HAS  A  START  TINE  ". 
1312."  WHICH  IS  LESS  THEN  THE  PRECEDING  START  TINE- ".312) 

END 


BlllO 

B1120 

Bll30 

iim 

B1160 

81170 

81160 

81190 

81200 

iim 

81230 

81250 

81250 

81260 

81270 

81280 

81290 

81300 

81310 

tm 

tun 

81360 

81370 

81380 

81390 

ntn 

IKI8 

81550 

81550 

81560 

81570 

81580 

81590 

81500 

81510 

81520 

81530 

81550 

81550 

81560 

81570 

81580 

81590 

81600 

81610 

31620 

81630 

81650 

81650 

81660 

81670 

81680 

81690 

B1700 

81710 

nwo 

81750 

81750 

81760 

81770 

81780 

81790 

Blaoo 

81810 

81820 

81630 

B185P 

81850 

81860 

81370 

81880 

81890 

31900 

81910 

B1920 

81930 

81950 

81950 

81960 

BI970 


82000 

82010 

82020 

82030 

82050 

§2050 

82060 

82070 

82080 

82090- 


2U 


I 

I 

I 


I 

5 

10 

IS 

20 

25 

30 

35 

40 


45 

50 

55 

60 


65 


70 


75 


so 


85 


<>0 


95 

100 


C 

C 

C 

r 

c 

c 

13 

40 


60 


13 


*«•*«***********•»*«*******«»******** 

THE  ARRAY  SI6  WITH  AVERAGED  »*♦•♦*•** 
THE  TINE  !nTERVAL--T1NE-.5  TO  TINE  • 

***«**{***«i****«*S  ************  •***•«* 

(S12i(ITTAG(I0ATAG(2OI»ISTU0»IH0> 

ii;5  tj;|8 

(ISTOATfROSTT* ANTHT(2ri»REFI20)« 
IG(20I 


90 

100 

110 

120 

C 

C 

C 

130 


140 

:  *0 
c 
c 
c 


SUBROUTINE  SIGNAL  (TINEvSIG) 

*****^*****i« •**•***•* •**•*•* 

THIS  SUBROUTINE  RETURNS 
DATA  FOR  EACH  SOURCE  FOR 

CONHON  ITYPE>NUNSR>  IDATA 

2  CALSC(200),CALHV<2QQi>CALSE 

3  NUNH0(IH0tYP«DTIN(5OOf2l»NU 

4  RDRNGIS.EOI/ALTACfDAYSfOAYR 

5  FRE0(20)f IHULT 
CONPLEX  REF 

DINENSION  SI6(20)tlSI6(20)>NS 
ICNT.O 

IF  .NUNO.EQ.O)  TLAST-1000000. 

IF  (NUND.EO.OI  GO  TO  30 

IF  TINE  IS  LE  END  TINE  OF  OTIN  ARRAY— 00  NOT  READ  OATA 

IF  (TINE.LE.DTININUN0t2) )  GO  TO  SO 

***************************************************************** 

BUFFER  IN  A  OATA  RECORD  INTO  IDATA 

TlAST-0TIN(NUN0t2l 

BUFFER  IN  (10,1)  ( IDATA!  II . 10 ATA( 512) ) 

IF  (UNIT(IO)T  40,220,40 
CONTINUE 

IF  NQRE  THAN  100  BACKSPACES— END  RUN 
IF  (ICNT.GT.lOOl  GO  TO  230 

IF  STATUS  WORD  NE.O  —  READ  NEW  RECORD 
IF  (IWRO(3).NE.O)  GO  TO  30 

****************************************  ***************** 

CALL  SUBROUTINE  SETTINE  TO  SET  UP  THE  TINE  ARRAY  DTIN  *••*•♦♦♦*♦ 

CALL  SETTINE 
GO  TO  10 

iF  TINE  GE  START  TINE  IN  ARRAY  OTIN  —  PROCESS  DATA 

IF  (TINE.GE.DTIH(l,2)-.5)  GO  TO  60 

IF  TIME  GT  END  TINE  OF  LAST  RECORD  AND  LT  START  TIME 
OF  THIS  RECORD -  THE  TINE  WAS  NOT  FOUND 

IF  ITINE.GT.TLAST)  GO  TO  70 

IF  TIME  IT  END  TIME  OF  LAST  RECORD  —  BACKSPACE  TAPE  10 

DO  6C  I«l,2 
BACKSPACE  lO 
CONTINUE 
ICNT-ICNTAl 
GO  TO  20 

WRITE  (6,240)  TINE,TLAST,0TINI1,2I 

IF»C 

IL»C 

NEXT«0 

DO  100  I-1,NUMD 

******************************************* fk********************* 

FIND  INDEX  INTO  OTIN  FOR  TINE  INTERVAL  ************************ 

IF  (TINE. GE. DTIN!  1,2) .OR. IF. NE.O)  GO  TO  90 
IF«OTIN( 1,11 

IF  (TiHE*.49,GT.0TIN( 1,2) I  GO  TO  100 
lL«OTiN(I-l,i) 

IF  (I.EO.l)  IL*0TIN(1,1I 

GO  TO  110 

CONTINUE 

NEXT.l 

DO  120  I>1,NUNSR 

isiG(n-o 

NSIGI l)-0 

****************************************************************** 

IF,IL  —  FIRST  AND  LAST  INDEX  IN  INTERVAL  **•♦♦***♦♦♦♦♦**♦•♦•*♦• 

J«IF 

K«IL 

IF  (K.EO.OI  K«IWRD(2) 

IF  (K.r.T.1920)  K»1920 
IF  (J.LT.5)  J»5 
00  150  I»J,K,2 
DO  140  l-i,NUNSR 

IF  |lWR0(I).NE.I0ATAG(L)l  GO  TO  140 
ISIC(LI>ISiG(Ll»IHRD(I-l) 

NSIGILUNSIGILItl 

CONTir.UE 

CONTINUE 

IF  NEXT  NE.  0  -  ENO  TINE  NOT  FOUND  —  READ  NEW  RECORD 

IF  (NEXT.EO.O)  GO  TO  200 

NEXT-0 

IF»0TIM(1,1) 

****************************************************************** 

««**^«**!****S*I***p$****l4*S*IS**** *******•*!***»*• ************* 

TLAST-OTIH(NUHO,2I 


C  330 
C  360 


560 

590 

m 

620 

630 

tn 

660 

670 

660 

690 

m 

740 

750 

760 

770 

760 

790 

600 

610 

820 

630 

l^sE 

S^3 

680 

690 


244 


BUFFER  IN  (lOfli  (IOATA(ll,IDATA(S12)l 

IF  lUNITIlO))  170»220,170 

CONTINUE 

IF  (IWROOt  .NE.O)  60  TO  160 


CALL  SUBROUTINE  SETTIHF  TO  SET  UR  THE  TINE  ARRAf  OTIN 


CALL  SETTINE 
00  ISO  i«i«Nuno 

IF  TtIMC».A9,6T.0TIH(I,2>I  60  TO  ISO 
IL*DTIN(I-l>i) 

IF  (I.EQ.ll  IL>OTin<l»li 
GO  TO  190 
CONTINUE 
NEXT-1 
GO  TO  no 

SET  AVERAGED  SIC  ARRAY 

DO  210  I-liNUNSR 
IF  (NSIGC n .EO.Oi  GO  TO  21D 
ISIG< It-FLOATI ISIGIIi l/FLOATIMSIGI  I>)^.5 


ISIGilt-FLOATd 

StG(I)-ISfG(I> 

RETURN 


URITE(6>2S0ITLAST 

STOR 

WRITE  (6,260)  TINE,DIln(l,2), lOAY 
STOP 

FORMAT  (6H  TINE  #F10.3«27H  NOT  FOUND  NO  DATA  BETWEEN  »F10.3,9H  AND 
1  ,F10.3) 

FQRNAT("  END  OF  FILE  ENCOUNTERED  -  TAREIO  -  LAST  TINE  FOUND  WAS”, 
1  F1C.3) 

FORMAT  (6H  T IME, FIO. 3, 26H  NOT  FOUND  FIRST  TINE  IS  ,F10.3,7H  DAY 

1-  ,110) 

END 


SUBROUT INE  RADAR  ( A  SR, DER, S TMSEC,  FT  MSEC , RNG,  I  ROLL,  I RTCH, RN6XY,  I RCK 

^ *,*,**«,*,,**,**•*»*****•***«****************•***** •*»*•*«****»* 

1HIS  SUBROUTINE  RETURNS  THE  A/C  DEPRESSION  ANGLE, ASPECT  *** 
ANGLE, slant  range,  AND  HORllONTAL  RANGE  AT  TIME  STHSEC  **** 
*,•*«*,*,,,*,,***«***•***, ,****»*,*****,,*•»*, *•****,*» 
COMMON  ITYPE,NUMSR,  ID ATA ( S12 ) , I TTAG,  I  DAT  AC  (  20) , I  ST WD, I WD, 

Z  CALSG(200),CALMV(200),CALSEI20,2),IHEDI90),  TINl(93).TIM2(90), 

3  NUMMO,  lHOTrP,OT IN ( BOO, 2 ) , NUMD,DELDB ( 20),HD2, 1  XX,  I  R ADAR , NUNRN6, 

A  RORNG(3,20),ALTAC,DAYS,OAYR,ISTDAT,ROSTT,ANTHT(20>,REF(20), 

5  FRE0<20),IMULT 
COMPLEX  REF 

DIMENSION  SUMV(6),0ArA(6) 

IRCK-0 


IF  NO  RADAR  TAPE 


GO  TO  200 


IF  (IRADAR.EO.IHN)  GO  TO  200 
DO  10  1-1,6 


DO  10  1-1,6 
SUMV(I)-0.0 

AN-0. 

ROLL-IROLL 

IF  (ETMSEC.NE.O.)  GO  TO  30 

READ  THE  RAOAR  TAPE  (TAPEll) 

READ  (11)  IWD.TTMS.IHR, ININ, TMSEC,ETMSEC,ISTHO,(DATA(I), 1-1,6) 
IF  (EOF(ll).NE.O.)  GO  TO  190 

USE  ONLY  DATA  RECORDS  WHERE  IWD  EQUALS  101 

IF  (IWD. NE. 101)  GO  TO  20 
IF  (OAYR.EO.1.0)  GO  TO  AO 
IF  (ROSTT.EQ.O.O)  ROSTT-TTMS 
IF  (POSIT. GT. TINS)  OAYR-1.0 
TT MS -TTMS-D AYR *26*3600 

IF  STATUS  WORD  NOT  VAI  ID  --  READ  NEW  RECORD 


IF  ( ISTWO.EO.lOH  /)  GO  TO  50 

GO  TO  20 

IF  (TTMS.LT.STMSEC)  GO  TO  20 

IF  END  OF  .5  SEC  DATA  REACHED 

IF  (TTMS.GE.(SIMSEC*.5))  GO  TO  BO 
DO  70  1-1,6 


PROCESS  DATA 


V?l}-^LIM 


UMV(I)*OATA(I) 


suMV(  n- 

AN«AN*1  ' 

GO  TO  20 


IF  NO  DATA  FOUND  - 

IF  (AN.EO.O.)  go  to  60 
DO  BB  1-1,6 
DATA! I )«SUMV(I ) /AN 
RNG-SQRT(OATA(  1 )**2*0A 


GE)  NEW  RECORD 


RNG-SQRT(OATA(  1 )**2*0ATA(2)**2*0ATA(3)**2)/6076. 
RNGXY-SORTf DATA(1)**Z*DATA(2)**2) /6076. 

IF  (IHOTYP.GT.l)  GO  TO  170 
HEAD-HD2 


69 

70 

79 

80 

89 

<>0 

V9 

100 

109 

no 

119 

120 

129 

130 

135 

140 

149 

190 

1 

5 


90 

100 


C 

I 

no 

c 

? 

120 

130 


140 

150 


I 

C 

160 


17C 

180 

190 

200 


210 

IIS 

C 

m 

260 


FOP  DISCRETE  POINTS  (CLOVER  LE4FS>P0LYG0N  FLIGHTSI  GO  TO  160 
IF  (IH0TYP.E0.2)  GO  TO  160 

FOR  ORBIT  FLIGHTS  GO  TO  STATEHENT  110 
IF  (IH0TTP.EO.3)  GO  TO  110 


THE  HEADING  IS  READ  IN  IN  TINE  INTERVAL 


FOR  FLYBYS  AND  RAOULS  START  HERE 

DO  100  UliNUlHD 

FINO  HEADING  FOR  TINE  SThSEC 

IF  (STHSEC. GE.TIH1(I>.AN0.STNSEC.LE.TIN2(IM  GO  TO  90 
GO  TO  100 
HEAO>IHEO( 1) 

GO  TO  160 
CONTINUE 

IF  TIME  INTERVAL  NOT  FOUND  GO  TO  STATEMENT  ISO 
GO  TO  180 

ORBIT  FLIGHTS 
DO  120  Ul.NUHHO 

II  >  UPPER  INTERVAL  OF  TIME 

II«I 

IF  (TIMldl.GT.STMSECI  GO  TO  130 

CONTINUE 

GO  TO  180 

H1-IHE0(II-1I 

HUlHEOlf  II 

0|LH-ABS(H1-H2I 

IF  (OELH.LT .200)  GO  TO  150 

IF  (HI. GT. 180. I  GO  TO  140 

H1«H14360. 

GO  TO  150 

DETERMINE  HEADING  BY  LINEAR  INTERPOLATION  OF  TIME 
H2«H2f360 

H£A0«(H2-H1I*(STMSEC-TIM1(II-1H/(TIMI(II)-TIM1I II-l) lAHl 
IF  (HEAD. GT. 360. I  HEA0«HEA0-360. 

IF  (HEAO.LT.O.i  HEA0«HEA0«360. 

SETUP  DATA  TO  CALL  SUBROUTINE  PANS 

TASP>ATAN2(0ATA(2I*0ATA(1I 1*57.3 
TOEP^ATANI  DATA!  31 /(BNCXY*60  76.11*57.3 
IF  (TASP.LI.O.I  TASP«TA$P*360. 

RNGM«RNG*1852. 

PTCH-IPTCH 

ALTAC-OaTAOI 

CALL  PAMS  FOR  ASPECT  ANGLE  AND  DEPRESSION  ANGLE 

CALL  PAMS  (ASP.OEPrRNGMtTASP.TDEP.HEAD.PTCH.ROLLI 
RFTUR 


IHOTYP 


(URN 

WRITE  (6.2401 
STOP 

WRITE  (6.250)  STMSEC. IHOTYP 
IRCK-1 
RETURN 

WRITE  (6.2601  TTMS 
STOP 
ASP«0. 

OEP-0, 

ETMSEC'O. 

^0^|20  I-1,NUMRNG 

IF  (RDRNG(2.I).GT.R0RN6(3.in  GO  TO  210 

IF  (STHSEC. GE.R0RNG(2. I). AND. STNSEC.LE.RDRNG(3.I)I  GO  TO  230 
GO  TO  220 
IF  (STHSEC. LE.F 
CONTINUE 

RNGXV«R0RNG(1. IN)*(  STHSEC-RDRNG(2. INI  I / (R0RNG(3« I NI-RDRNG( 2. INI  I 
RNG«SQPT(RNGXY**2*ALTAC/6076.**2) 

FORMAT  (1X.24H  ILLEGAL  HEADING  TYPE  >121 

FORMAT  (1X.4HTINE.F8.1.24H  NOT  IN  HEADING  INTERVAL. ITHFOR  HEADING 

F0RHAf^fix.66HEN0  OF  FILE  ENCOUNTERED  ON  RADAR  TAPE— LAST  TINE 
1  FOUND  WAS  .F10.3I 
END 


.R0RNG(2>II. AND. STHSEC. GE.RDRNGO. in  GO  TO  230 
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740 
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770 
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SUBROUTIME  CAL  (IPCAtl  ,  ^  ^  ^  ^  ^  E 


THIS  SUBROUTINE  SETS  UP  THE  SIGNAL  CAUtRATION  ARRAYS  ******** 
CALMV  AND  CALSG  **••*«{»**{*{•{{*****************••*«** 

COMMON  ITYPE.NUHSR.  10ATA(512)>  ITTAG. I0ATAG( 201 • ISTHO. IWO. 


2«i 


pttaMLMM 
r  c9IWiW!a 


10 


15 


20 


25 


30 


35 


50 


55 


50 


55 


60 


65 


70 


75 


80 


85 


90 


10 

C 

c 

c 

c 

c 


20 


30 

50 

50 

60 


70 

80 

no 

m 

130 


2  CALSC(200>>C*inV(200)*C*LSEI20»2l>IHEO(90l>Tini(90l,TIH2(90)>  £  70 

3  NUMHDtIH0TYP*0TinCS00f2i>NUN5tQELDB12gieHD2>IXX»iR5DAR»NUMRNS>  E  80 

5  RDRNG(3*2Qt<*LTACfDAYS>OAyR«iST6*T«RQSTT*ANTHT(20),REE(201>  E  90 

5  rREQ(20)<INULT  E  100 

CORRLEX  REF  E  110 

QlnENSION  IHR(5C)>ININ(50)>ISEC(5CI>NSRID(20>>SI6(20I.TRSG(2,100)  E  120 

INTEGER  CALSE  E  130 

IXX.l  I  J5J 

IF  PRE  AND  POST  CALS  —  SET  IXX  ■  2  E  160 

E  170 

IF  (IPCAL.EQ.IHY)  IXX-2  E  160 

JO  10  I-t>100  E  190 

DO  10  K-lt2  £  200 

TPS67K*lt«0.0  E  210 

00  50  lX«l>ixX  E  220 

J»1  £  230 

E  250 

READ  OBN  AKO  TINES  FUR  EACH  SOURCE  L  250 

E  260 

DO  20  I«1«NUNSR  E  270 

•t****************************************************************  t  280 

READ  IN  THE  NUHBER  OF  SETS  OF  CALIBRATION  SIGNALS  **••♦•*••*•*  £  290 

AND  THE  CORRESPONDING  TINES  *♦♦*♦♦♦**»♦*••*♦•♦♦♦♦♦♦*♦*♦*♦♦♦»♦  fc  Jto 
t*****************************************************************  E  310 

£  320 

REAO  NUMBER  OF  CAL  POINTS  FOR  EACH  SOURCE  -  ARPHP  CARO  16  E  330 

E  350 

READ  (S«eOI  CALSEIfeE)  £  350 

WRITE  (6>90l  CALSEfI»2>  E  360 

£  370 

SET  STARTING  AND  ENDING  ENTRY  OF  CAL  POINT  IN  ARRAY  CALSE  E  380 

E  390 

CALSEIItll-J  E  500 

JJ-J»CALSEtI*2)-l  E  510 

!tOT>JJ  £  520 

**********6****»**6****************6*A**6****************6*6«***«*  E  530 
REAO  IN  THE  SIGNAL  IN  DBN  AND  THE  £  AAO 

CORRESPONDING  TINES  **•*♦****•••*•••♦•••••*♦•♦♦*»••****♦•♦*****  £  550 

*******«****«*****«**«*««*********6********«*«*6*«**«********«****  £  <,60 

E  570 

REAO  SIGNAL  OBN  AND  TINES  ARPHP  CARD  17  E  580 

E  590 

READ  (5tltO)  (CCALNVIK)«IHR(K|>ININ(KI>ISEClKt)fK*J«JJ)  £  500 

WRITE  (6«100i  (TCALNV(K),IHRIK)>ININ(K>>ISEC(Kn>K>J, JJ|  £  510 

J«JJ*1  £  520 

DO  50  I*UNUNSR  E  530 

L«CALSE(i«l)  E  550 

K«CAl!E(I«2I  E  550 

LL«L*K-1  t  560 

E  570 

GET  SIGNAL  AT  CALIBRATION  TINE  TSEC  £  580 

t  590 

DO  50  8«1>LL  E  60C 

TSEC>IHR( J)*3600.AININ(J)*60.5I$EC(  Jl  E  610 

tT«0  £  620 

E  630 

REAO  2  SEC  OF  DATA  FRON  SIGNAL  TAPE  (TAPEIOI  E  650 

£  650 

DO  30  23*l/6  £  660 

CALL  SIGNAL  (TSECfSIG)  £  670 

TT-TT*SIG<I»  e  680 

TSEC«TS£C*.5  E  690 

TPSG(IX>J]-TT/5.0  E  700 

WRITE  (6»i30i  (IHRI J)>IHIN( J)»ISECfJ)>CALHV(J)>TPSG(IXf JIf J>1»LL)  b  710 
CONTINUE  b  720 

DO  60  I-lrlTOT  £  730 

CALSG(H-aPS6Il,I»»TPSCI2,I)  )/IXX  £  750 

DO  70  I-1,NUN$R  b  750 

L*CALSE(ifl)  E  760 

LL>L-l«CALSE(It2)  E  770 

WRITE  (6>120i  t»((CALNV(J)>CALSC(J>)>J<L,LL)  £  780 

CONTINUE  E  790 

b  80U 

FORMAT  (1101  E  810 

FORMAT  (1H0>35HTHE  NUMBER  OF  CALIBRATION  SETS  ARE  >1X.I2)  c  520 

FORMAT  (lH0t33HfH£  CALIBRATION  SIGNAL  IN  OBN  IS  t FIO. 5> /> 1X> 26HTHE  E  830 
1  CORRESPONOING  TINE  IS  »  I2»  fx,  3HHRS>1X,  I2»  IX,  AHNINSdX,  I2*1X.  5HSEC  E  850 
2S>  E  650 

FORMAT  (5(F10.6,5X,3I2))  fc  860 

FORMAT  (IHOflBH  SOURCE  NUMBER  ,I2,/f5(20H  SIGNAL  03N  SIGNAL  l#/<2  £  870 

10(10(1X,F9.1I,/)I  E  880 

FORMAT  (1H0,9(22H  TINE  Z  DBN  SIGNAL  I, /> 20 (5( IX, 312, IX. F6. 1# IX.  E  890 
lF7.1l,/,Ti  i  900 

END  E  910- 


1 

9 

10 


FUNCTION  VCAL(X,I)  f 

•  ******1**««««*6111*11*«*11»***»******1***1«,*1****,1*1*****A*1*«  f 

FUNCTION  VCAL  PERFORMS  SIGNAL  CALIBRATION,  I-SOURCE  NUMBER-**  f 

•**«**s«««**^i**i****i*i*$*Us*i»i******nn*********************  ^ 

CONNON  ITYPE,NUHSR,IDATA1912),ITTAG,IOATAG(20I,ISTWO. IWO,  F 

2  CALSC(200>,CALNV(2Qgi,CALSE(20,2),tHEO|90l,TINl(9g),TIM2(90),  F 

3  NyHHD,lH0TVP,0TIN79OO,2i,NUM0,0|L0BI2Oi,HD2,IXX,IRA0AR,NUNRNG,  F 

5  R0RNG(3.2gi,ALTAC,OAYS,OAYR,IST5AT,ROSTT,ANTHT<2O),REF(2O),  F 

9  FREQ<26),IHULT  F  100 

COMPLEX  REF  F  110 

INTEGER  CALSE  F  120 


10 

20 

30 

50 

50 

60 

70 

60 

90 


24T 


15 


20 


25 


30 


35 


1 


5 


10 


15 


20 


25 


30 


35 


<i0 


<.5 


50 


55 


60 


1 


5 


C 

c 

c 


10 

20 

C 

C 

c 

30 


<.0 

50 

60 


C 

c 


c 

c 

c 

10 

20 

30 

<.0 

50 

'0 

70 

SO 

00 

100 

110 

120 

130 

ISO 

15C 


C 

C 

C 

c 

c 


NN«CALSfc( !t2) 

MS-CALSEl 1,1) 

00  60  ^*1,NN 
J 

N'HM*  J  -1 

16  SIGNAL  CAL  A«DAr  IS  ASCENDING  —  GO  10  10 

IF  (CAlSGIfM)  .GT.CALSG(Hf'ANN-l)  )  GU  TO  10 
IF  (*.CT.CALSG(N) 1  GO  TO  50 
GO  TO  20 

IF  ( X .LI.CALSGI N) )  GO  10  50 
IF  (J.tl.l)  GO  TO  SO 

PEKF03N  LINEAB  E  K  T  H  A  t*OL  A  T  I  ON  TU  i  I  ND  SIGNAL  CALISRATION 

VCAL-CAL6V(N-n*(CALNV(N)-CALF>V(N-l)»*IX-CALSGINTl))/(CALSG(N)-CAL 
lSGIN-1 ) 1 

return 

N«N»1 
GO  TO  30 

IF  (J.ea.NNl  GU  1C  30 

CONI INUE 

END 


FUNCTION  IWBOdl 

THIS^FUNCTI0N^EXTRACTS^THE^I^TH^16^BU^W0RD^F jaj^THE^60^3I I^A* 

COMNON  ITVPF.NUMSR.  IDA7A (512 ) . ITTAG.  IOAI AGI 20) . ISIWD, 1 WD, 

2  CALSG( 200) ,CALMV( 2001 ,CAl  SE ( 20»  2 1 » 1  HE  D ( NO ) . T 1  Hi ( GO ) , 1 1 H2 1  U 0 ) , 

3  NUMHO.lHOIrP,nTIH(5CO,2).KLfiO,UELDBr20),HD2,U<,lRADAR,NUMRNG, 

S  RORNGI 3,20) . At  TAC.OAYSfOATR, lSTCAT,R3STT,ANTHTI2i) , REF (20) , 

5  FReu(20), InULT 

COMPLEX  REF 
112*77773 
Ml*177777a 

MS  *3773 

calculates  RHICH  60  BII  WORD  lUE  16  Bit  KURD  RILL  BE  IN 

•1>FL  JATI  I  )/15.0000l 
II«1-M«15 
JJ*SAM1l 
I  Wl  »  I  OATS  (  3  J  I 
IW2"  lOA  TaI  J  J*,  ) 
tW3» lOAt A( J J  *2) 

IWS>IOATA(JJO) 

IF  (  I  I.LT.I.UR. II.&I.15)  II»1 

GO  to  THE  CORRECT  STAIEMEM  TO  GET  THE  RL)RO 

GO  TO  (  10,  20, 30, SO,  50, 60,  70, BC,9C,  100, 1  10,  120, 130, ISO, 150),  II 
IRP0*AN0(SH1FT( IW1,16),M1) 

RETURN 

lURO* ANO( SHIFT ( 1WI,32;,M1) 

RE  turn 

IRRD«AN)ISHIFT(1W1,SB),M1) 

RF  Turn 

IWRD.ORtSHim  ANO(  IWl  ,M2),S  ),  ANO(  SHIFT!  IR2,S)  ,M3)  ) 

RETURN 

IWPO-ANO(SHIH(  IR2,20)  ,11) 

RF  TURN 

IRPD-AN0(SHIET( 1W2,36),m1) 

Return 

IRRC*AN0(SHlfT(lR2,52),Ml) 

RETURN 

I  RRC-Ok  (  iMI  I  I  (  ANl)(  I  R2,  MS  ) ,  B  ),  AND(  SMifT(lR3,8),MS)) 

return 

I  VRC-  ANn(  S  II  F  T(  I  R3,  2S  )  ,M1  ) 

PFTUPN 

tRRD*ANJ(5'UfT(  IR3,S0),M1) 

RETURN 

1WP0'AN0(SH1FT(IW3,56),M1) 

P  F  TUR  N 

I WRD-OR ( SUIF  T ( AND(  I  W J,  13)  ,  I  2) ,ANOI SHIFT  I  I  MS, 12) ,M2) ) 

RE  TUR  N 

IWRP*AilD(SNlFr(lWS,28),Ml) 

RE  TURN 

rWRO*  AMIJ(  SHI  FT  (  IWS,  SS  )  ,  Ml) 

R  E  turn 

IMR0-AND(  IMS, Ml) 

RETURN 

END 


SUBROUTINE  StTTlME 


This  subroutine  sets  array  diim  with  the  entry  location  and  ♦* 
total  time  in  seconds  at  each  HiLLIStCONU  ENTRY  IN  THE 
SlGNAl  RECORD  (FRC1  TAPFiC)  *♦♦•*****»**»»*•**•***♦♦♦♦*♦♦**** 


COMMON  ITYPE.NUMSR,  I C A T A < 512  I  ,  1 T T AG, I D AT  AG  I  2 0 ) , I  ST  WO, I WD , 

2  CAL  SG(200) ,C ALMV(2  00  )  ,  C  Al  SE  1 20,  2  )  .  1  HEO ( 90 ) , T I  Ml ( 90 ) , T I  M2 ( 90 )  , 

3  NUMHO, IHOIYP, DTIMISOC,  2)  ,  NUMO, 0 1 L DB I  20 ) , H02 , I  XX, I R ADA R,  NUMR NG , 


13C 

F 

IS'j 

F 

150 

f 

16C 

f 

170 

F 

iBO 

t 

190 

F 

200 

F 

210 

F 

220 

F 

230 

F 

2S0 

f 

250 

F 

260 

F 

2  70 

f 

280 

h 

290 

F 

300 

F 

3)0 

F 

32C 

F 

330 

t 

3S0 

F 

350 

f 

360- 

G 

10 

G 

20 

G 

30 

G 

SC 

u 

50 

G 

60 

G 

70 

G 

80 

G 

90 

G 

ICC 

G 

1  1  u 

G 

i20 

G 

1  30 

•J 

ISO 

150 

G 

IF  0 

0 

17.' 

0 

IMr 

y 

190 

N/ 

t’ 

J  G 
2iu 
2.-J 

G 

G 

23C 

0 

2S0 

G 

250 

G 

260 

u 

2  70 

230 

G 

2  9J 

G 

300 

G 

31C 

G 

320 

G 

330 

u 

3S0 

G 

350 

G 

360 

370 

G 

38C 

N* 

3  9., 

LA 

soo 

LA 

Sio 

G 

S2C 

G 

■,3c 

SS  J 

s50 

s  (.0 

G 

S70 

s  30 

SOO 

5  00 

■J 

510 

o 

52o 

530 

,, 

5  SO 

La 

55.: 

560 

(, 

57C 

\j 

580 

O 

590 

G 

6C0 

G 

61C- 

H 

10 

■i 

2C 

H 

30 

H 

so 

H 

50 

h 

oC 

H 

7C 

H 

80 

r« 

90 

248 


i 


lo 

<>  RDPNGI  3|20  ) ,  ALTAC>DAY$.  DA  YR  >  1  STDAlf  R3STT,*NTHT  1  201*  REFI  20  1  >  H 

ICO 

5  FREQ( 201 . INULT 

H 

110 

C 

CQMPIEX  REF 

H 

H 

ill 

c 

GET  END  TIME  OF  SIGNAL  RECORD 

H 

150 

15 

c 

H 

150 

CALL  TTINE  ( IHR > NIN . 1  SEC » TS EC  i 

H 

160 

IFT*0 

H 

170 

ASEC-0. 

H 

180 

lLAST»9<?9q 

H 

190 

20 

NN«IWPD(2I 

H 

200 

IF  (NN.GT. 19201  NN«1920 

H 

210 

c 

H 

220 

c 

c 

COUNT  THE  NUHSER  OF  SECOND  CHANGES 

IN  RECORD  H 

H 

m 

25 

DO  20  1*5. NN. 2 

H 

250 

IT-MILS£C(IW»0II-1»  » 

H 

260 

IF  (IWROt n.NE.ITTAG)  GO  TO  20 

H 

2  70 

IF  (  IT. 6T. HAST!  GO  TO  10 

H 

280 

3SEC*ASEC*1 

H 

290 

30 

10 

ILAST-IT 

H 

300 

20 

continue 

H 

310 

ILAST»9999 

H 

320 

ST*AINT( TSEC )-ASEC 

H 

330 

M*0 

H 

350 

35 

11*1 

H 

350 

C 

H 

360 

C 

SET  ENTRY  IN  OriN(n<l)  AND  TlnE  OF 

ENTRY  IN  0TIf1(N#2)  H 

370 

C 

H 

380 

DO  50  I*5.NN»2 

H 

390 

IF  (IWROdI  .NE.ITTAGI  go  TO  4C 

H 

500 

H.M*1 

H 

510 

IT«NILS£C(IWRDI  I-II  ) 

H 

520 

0TIN(M,1)*I 

H 

530 

0TIM(M,2I*ST*FLQATI  ITl/lOOC. 

H 

550 

<t5 

IF  (IT.GT.ILAST)  GO  TO  30 

H 

550 

OTIH(M,2I-DIIM01,21*1. 

H 

560 

ST.$T»1 . 

H 

570 

30 

ILAST-IT 

h 

580 

<i0 

CONTINUE 

H 

590 

50 

c 

H 

500 

c 

NUHO  •  NUH8EP  OF  VALID  ENTRIES  IN 

ARRAY  DTIN  H 

510 

c 

H 

520 

NUMD»N 

H 

530 

END 

H 

550- 

1 

SUflBOUTiNE  TTIME  (  I  T OTHR , I T OT hi N, I TOTSE C .  T TM) 

I  10 

C 

tt*********************************** ***************************** 

I  20 

C 

this  subroutine  EXTRACTS  THE  END  TINE  OF  THE  SIGNAL  ***•♦•»♦*♦ 

1  30 

C 

RECORD  FROM  THE  LAST  58  BITS  *»♦♦**♦»*******♦••*♦******»»*♦ 

I  50 

5 

r 

tt **************************************************************** 

COMMON  ITYPF.NUMSR.  I  0  A  T  A  (  5 12  ) »  ITT  AG»  lOAT  AG  I  201  »  I  ST«D»  I  W0» 

2  CALSG(200)«CALMV(2CO)>CALSE(20»2)elHEO(90)>Tlnil90) >TIM2(90>f 

3  NUMHD.IH0TYP.0TIM(500*2) .NUhO,DELDB(20)»HD2» IXX»  IRADAR.NUMRNG. 

5  R0RNG(3,2O),ALTAC.OAYS.OAYR,  IST0AT»R5STT.ANTHTI2OI»RtF(2O>. 

5  FRf0(20t. IMULT 

I  50 

I  60 

1  70 

I  80 

I  90 

10 

I  100 

complex  REF 

I  110 

MW01*78 

I  120 

l«mR0(2) 

I  130 

MW02«179 

I  150 

15 

C 

MWD3*3B 

1  ISO 

I  160 

c 

IF  DAY  EOUAL  START  OAT  PLUS  1  —  SET  DAY  •  1  TO  ADD 

I  170 

c 

25X3600  SECONDS  TO  TOTAL  TIME 

I  180 

c 

I  190 

20 

IF  (DAYS. EQ. 11  GO  TO  10 

1  200 

IUOY.ANDI  SHIFTUWROT  I».-6>  ,  MR 02) 

I  210 

ITDY.ANDI  SHIFK  IWRDI  I  )f-10)  ,Mt,D2) 

I  220 

I HOY* AND! SHIFT! I WRD 1  I )  ,-l 5  )  , M WD 3) 

I  230 

lOAY-IUDY*! TOY*10»IHOY»100 

I  250 

25 

IF  ((OAY.EQ.ISTOAT*!)  DAYS»1.C 

I  250 

10 

ISFC»AN0(  IWR0<1»1»,MMD2) 

I  260 

ITSlC»AND(SHIFTt IWRO! lAl) *-5) ,MUD1) 

I  270 

IiniSEC-lSEC+10»ITS£C 

1  26w 

IM IN* AND!  SHIFT!  IWRD ! 1 *1 ) * -7 ) . MWD2 ) 

I  290 

30 

1TMIN*AN0(SHIFT! IWRO! I»l ),-ll ),MWDn 

I  300 

ITatMlN*IMIN»10»ITMIN 

I  310 

IHR>ANO( IWRO! I ) ,HRD2» 

I  320 

ITHR-ANOISHIFT!  IWRD! n»-5),MWD3) 

1  330 

ITOTHR*IHR*10»ITHR 

I  350 

35 

TT M.ITCTHR* 3600*1  TO TMIH*60»ITCTSEC AFLOAT (  MILSEC! IWRD!  1*21 ) ) /1030.  ♦ 

1  350 

10AYS*25. 0*3600. C 

I  360 

return 

I  370 

END 

I  380- 

1 

SUBROUTINE  lOFUNTICR) 

J 

10 

DIMENSION  ICR!e»,  ICHAR!117I>1CH!1C) 

J 

20 

INTEGER 

SSCRIP 

J 

30 

COMMENT 

J 

50 

5 

COMMENT 

J 

50 

COMMENT 

"SSCRIP" 

J 

60 

COMMENT 

EXTERNAL  FUNCTION  TO  SECURE  TABLE 

SUBSCR IPT 

,1 

70 

comment 

JOMMENT 

TO  obtain  DISPLAY  CODE  EQUIVALENT 

OF  AN  EBCDIC  CODE 

J 

J 

80 

90 

10 

COMMENT 

THE  FOLLOWING  TABLE  DESCRIBES  THE  RELATIONSHIP  BETWEEN 

J 

100 

COMMENT 

ALPHANUMERIC  ICHARACTERSt  THEIR  EBCDIC 

CODE!  IN  BOTH  OCTA 

J 

110 

15 

?0 

^5 

30 

35 

50 

55 

50 

55 

60 

65 

70 

75 

80 

1 


COMMENT 

AND  DECIMAL 

FORM! 

•  INTERNAL  DISPLAY 

CODE  lOCTALI 

»  AND 

J 

120 

COMMENT 

THEIR 

ORDER 

NUMBER  tfITHIN  •CICHAR"  ARRAY... 

SORT 

IS  BY 

J 

T  30 

COMMENT 

ORDER 

NUMBER. 

J 

m 

COMMENT 

J 

COMMENT 

J 

160 

COMMENT 

ORDER 

KHAR 

EBCOI! 

DEC 

DSPC  ORDER  ICHAR 

EBCDIC 

DEC 

DS 

J 

170 

COMMENT 

J 

180 

COMMENT 

COMMENT 

xl 

• 

m 

n 

N 

0 

325 

326 

i 

1 

16 

17 

J 

J 

190 

200 

COMMENT 

i? 

) 

76 

52 

P 

32  7 

2] 

5 

20 

J 

210 

COMMENT 

♦ 

80 

55 

85 

0 

330 

6 

21 

J 

220 

COMMENT 

28 

t 

133 

91 

53 

85 

R 

7 

J 

COMMENT 

29 

% 

135 

92 

57 

95 

S 

352 

2< 

'6 

h 

J 

COMMENT 

33 

150 

96 

57) 

95 

T 

353 

2 

7 

25 

J 

250 

COMMENT 

35 

/ 

151 

97 

50 

96 

U 

355 

2. 

!6 

25 

J 

^2^78 

COMMENT 

55 

\n 

97 

V 

355 

229 

26 

J 

COMMENT 

55 

! 

155 

51 

98 

w 

356 

230 

27 

J 

280 

COMMENT 

6C 

■ 

173 

123 

55 

99 

X 

35  7 

21 

1 

TO 

J 

290 

COMMENT 

61 

A 

301 

193 

01 

00 

Y 

350 

232 

31 

J 

300 

COMMENT 

62 

B 

302 

195 

02 

101 

08 

1 

351 

2: 

3 

3? 

J 

310 

COMMENT 

63 

C 

303 

195 

03 

0 

360 

250 

33 

J 

320 

COMMENT 

65 

0 

305 

196 

05 

09 

1 

361 

251 

35 

J 

330 

comment 

65 

E 

305 

197 

05 

10 

2 

362 

252 

35 

J 

350 

COMMENT 

66 

F 

198 

C6 

r}2^ 

3 

363 

253 

36 

J 

3  50 

COMMENT 

67 

G 

199 

07 

h 

365 

255 

37 

J 

i^7g 

COMMENT 

COMMENT 

66 

H 

200 

13 

5 

365 

255 

50 

J 

69 

I 

311 

201 

11 

15 

6 

366 

256 

51 

J 

3H0 

COMMENT 

77 

J 

’21 

209 

n 

7 

367 

257 

52 

J 

390 

COMMENT 

78 

K 

322 

210 

13 

8 

370 

258 

53 

J 

500 

COMMENT 

79 

L 

323 

211 

1* 

17 

9 

371 

259 

55 

J 

510 

COMMENT 

80 

M 

325 

212 

15 

J 

^20 

COMMENT 

J 

530 

COMMENT 

J 

550 

COMMENT 

ALL  CtCHAR 

ELEMENTS  NOT  USED  ARE 

SET 

TO  \ 

J 

550 

C 

J 

560 

C 

J 

570 

05 


,lH(.l5*lH\>lH-> IHA, 1HB>IhC»1hS»1KE>1HF> 1H6> 
iHH«rHt>7MiU>lHJt  lHK«lHL>lHn.lHN»lH5>lHP>l.H0>  1145 1  8 » 1  H\>  |hS> 


1H7>1H8«1H9  / 


COMHENT 
COHMENT 
COMMENT 
COMMENT 
COMMENT 
COMMENT 
COMMENT 
^OMMENT 


IN  AOgiTinN,  THE  EOLIOWING  ARE  INCLUDED  TO  PROVIDE  024  KP 


in  «uuiiiun>  me 
CQMPATIBUITY.... 


ORDER  ICHAR  EBCDIC  DEC  OSPC 


ORDtR  ICHAR  EBCDIC  DEC  OS 


15 

15 


115 

116 


77 

78 


51 

55 


30 

59 


» 


135 

176 


93 

126 


52 

55 


10 

20 

100 


KK»2 

DO  20  j>ue 
DO  10  I>1>10>2 
KK«KK»i 
ITP-IWRO(KKI 

H2»AND«  ITP»  37781 

3*AN0(SHIFT(ITP.-8),377B) 

iCH<  n»  ICHAR  (sscR  IP  I  nan 

ICH(ni)«ICHAR(SSCRIP<IT2n 

ENCODE! 10. 100. IFMII <ICH(  II.  I>l>10i 

ICRIJniFMT 

CONTINUE 

FORMAT! lOAl) 

END 


FUNCTION  MILSECin 

THIS  FUNCTION  SETS  THE  MILLISECONDS  FROM  THE  12  BIT  DECIMAL  CODING 
K«17B 

MILSEC-ANO!  I . K I flO* AND!  SHIF T!  I.-5 1 . Kl R1 DO* AND  IS  HIF T 1 1 . -8) . K I 
END 


580 

590 

500 

520 

530 

550 

55G 

560 

570 

580 

590 

600 

610 

620 

630 

650 

650 

660 

67£ 

660 

690 

m 

720 

730 

750 

750 

760 

770 

780 

790 

800- 


10 

20 

30 

tl 

n. 


10 

1 

5 


INTEGER  FUNCTION  SSCRIPIIII 

THIS  FUNCTION  IS  USED  TO  FIND  THE  INDEX  IN  ICHAR  OF  lOFUN 
I-II 

IF  !I.GT.259.aR.I.LT.65l  n65 
SSCRIP-I-53 

IF  ISSCRIP.E0.63I  SSCRIP-59 

IF  ISSCKIP.GT.1Z9)  SSCRIP-SSCRIP-69 

END 


SUBROUTINE  PAMS  ! VAO.  VEO. TR . T A, T£ , VH> VP > VR  I 

*«»***»********«***#****#•**»**  **********************  ARRRRRRRRRRR* 

THIS  SUBROyTIME  CALCULATES  THE  DEPRESSION  ANCLE  AND  ASPECT  »♦♦♦ 
ANGLE  OF  THf  A7C.  DOCUMENTATION  OF  THE  EQUATIONS  CAN  BE 
FOUND  IN  RAOC-TR-75-175  PAMS  AIRBONRE  INSTRUMENTATION  *♦*♦♦♦ 


30 
50 
50 

80 
90 
L  100- 


10 

18 

50 

50 


2SI 
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PAMS  AMS  COORDINATC  TRAMSFE*  €OUATIONS  ♦*♦♦*♦♦*♦*•***♦•♦ 
INPDTl 

TP  ■  TRACKER  RANGE  TO  VEHICLE  IN  METERS  •♦♦♦♦♦**•♦♦*♦♦*• 

TE  •  TRACKER  ELEVATION 

TA  -  TRACKER  AZIMUTH  ♦•*•*♦*•♦♦♦**•♦♦♦♦♦♦*♦**•****♦»♦♦ 

VH  ■  VEHICLE  HEADING 

VP  ■  VEHICLE  PITCH  ********************************* 

VR  •  VEHICLE  ROLL 

OUTPUT  I 

V£  •  VEHICLE  ELEVATION  I  TO  SITE  I  ••♦♦♦♦♦♦♦♦♦*♦♦♦*♦•* 

VA  •  VEHICLE  AZIMUTH  ♦♦♦*♦♦*♦*♦♦♦*♦♦*♦** 

RE  ■  RADIUS  OF  SPHERE  HAVING  EARTHS  VOLUME  IN  METERS  ** 
VYl-  FIRST  YAW  ANGLE 
VPl-  FIRST  PITCH  angle  ICOHPUTEOI 

VV2»  SECOND  YAW  ANGLE  ICOMPUTEDI  ♦♦♦♦♦•♦*♦♦**♦♦*♦♦»♦♦♦ 

CORRECTION  FOR  MAGNETIC  VARIATION 
VH>VH»1S. 

IF  (VH.GT.I60.I  VH*VH-360 

6TR-.OI7A5329252 

RE-6361A96.A21 

yH«yH*OTR 

TE«TE*DTR 

TA«TA*OTR 

VP«VP*DtR 

VR>VR*DTR 

R180«ie0.*DTR 

R36O«360.*0TR 

RORE>fP/RF. 

SE«SIN(TEJ 

CE«COS(TEI 

SA>SIN(TA) 

CA>COS(TA) 

H>RE*(SORT(I.*ROREA|RQRE«2.*SE)l>l. ) 

GAHAN«ATAN( TRACE /IREaTR*SE) I 

CGN-COSIGAMANI 

SGN«SIN(GAMANI 

VP1--GAMAN 

IF  (TA.GT.RISO)  GO  TO  ID 

TAI-TA 

GO  TO  20 

IA1-R360-TA 

STAl»SlN<tAl) 

CTAl-COSITAll 

VR>ATAN2(STAIt  (.697S27S*$GN-C  GN*CTAlt  ) 

IF  (TA.LE.RIAO)  GO  TO  30 

VH1-R180-VB 

GO  TO  AO 

VHI-VB-R160 

VY2*VH1aVH 

XO«-TR*CE 

Y0«0. 

ZO*TR*SE 

ROTATION  SEOUENCE  IS  VP l» VY2» VP»VR  ♦•♦♦••♦*♦*♦***♦ 

CPI'COSIVPI) 

S»i-SIN  vpI) 

CY2«C0S( VV2I 
SY2»SIN(VY2I 
CP«COS(VP» 

SP«5IN( VP) 

CR«COS(VR? 

SR-SINIVR) 


■XOA(CP1*CY2*CP-SPIASPIAYO*SY2*CP-ZO*IS 
•K0*(CP1ACY2*SP*SR-CPI*SY2*CRASP1PCP*SR 
ISP1*CY2*SP*SR-SPIASV2*Cr-CP1*CP*SR) 

■XQ*(CPlACY2*SP*CRACPi*SY2*SR»SPl*CP*CR)AY0A(SY2ASP*CR-CY2*SR)-2 
(SPI*CY2A$P*CR*SP1*SY2ASR-CP1*CP*CR) 

•ATAN2jVV>Xy) 

■ATAN2i-ZV»S0RT<XV*KV*YV*YV)) 

}«VA/OtR 

IVAO.LE.O.O)  VAD«VADa360.0 
)»VE/0TR 
TURN 


SPI*CY2PCP»CP1*SP1 

RItYJ*?SY2*SPASRACl 


M  60 
N  70 

M  ao 

M  90 

3  m 

N  120 
M  130 
M  lAO 
M  15C 
M  160 
M  170 
M  100 
M  190 
M  200 
M  210 
M  220 
M  230 
M  2A0 
M  250 
M  260 
M  270 
M  200 
M  290 
M  300 
M  310 
N  320 
M  330 
M  3A0 


M  3SC 
M  390 

M  ao5 
M  AlO 
M  A20 
N  A30 
M  AAO 
M  A50 
M  AbO 
M  A70 
M  ASO 
M  A90 
M  500 
H  510 
M  520 
N  530 
M  5A0 
M  550 
M  560 
M  570 
M  580 
M  590 
M  600 
M  610 
M  620 
M  630 
M  6A0 
M  650 
M  660 
M  670 
M  68f 
M  69C 
H  700 
M  710 
M  720 
M  730 
M  7A0 
M  750 
M  760 
M  770 
H  700 
M  790 
M  300 
M  810 
M  820 
N  830 
M  8A0- 


SU8R0UTINE  MULTPTH  (RANG  > ISR> ICKHPfOBGAIN) 

•**•*«««*«***}** A**A**A*****» **•»*••*•» *•••••*•**** **•***•**«•»* 

THIS  SUBROUTINE  CALCULATES  THE  EFFECT  OF  MULTIPATH  ♦♦♦♦• 
RAOIATION,  THE  TECHNIQUE  USED  IS  DOCUMENTED  IN  THE 


FLIGHT  MEASUREMENT  OF  AFC  ANTENNA  RADIATION  PAT 

**•*••**•***«•{«•***•*•«*•#*••*** b*********** 

COMMON  ITYPEfNUMSR,  lOATAI 512 1* ITTAC. lOAT AGI 201 »ISTHO> IHO» 

1  CALSG(200)fCALNyl2O0)»CALSEi29>2)>IHf0(Q0I.T|NiT9OI.TIN2(90)> 

I  NUHH0>lH0IVPt0ItN79OOc2)«NUMDtO£L6ll2OleH02tIXlf IRAOARfNUHRNG* 
I  RDRNG(3{20)f ALTACiOAYSiDAYRf ISTDATfROSTTf ANTHT(20)»REFI20I» 

A  FREOIZOIflNULT 
COMPLEX  NN.RRfRTfREF 
NN-REMISR) 


TJERNS^REPORT^ 


cunrctx  nnaKRfRi 

NNrREMISR) 

Hl>ANTHf(ISR) 

TALT-ALTAC 

POL-1 

RAN6E-lilA9*RAN6 


N  50 
N  60 
N  70 
N  80 
N  90 
N  100 
N  110 
N  120 
N  130 
N  lAO 

3  ill 

N  170 
N  Iso 


2U 


20 


25 


30 


35 


t>0 


45 


50 


C 

c 

c 

10 

20 

30 


40 

50 


FIND  CORRECT  DIST  I 

CALL  ODIST  (RANGE>HliTAlT»01STl»ICKNI>,H2RRI) 
IF  (ICKNP.EQ.l.OI  GO  TO  40 
ARG2-H1/(5280.0*0IST1» 

FIND  ANGLE  OF  INCtO/REFLEC 

ST-(ATAN(ARG2)> 

RT»CSQRT(NN-C0SISII**2I 

FUD  REFLEC  COEF  AND  ANGLE 

I? 


(NN*SIN(SII4RTI 


RR>(NN*SIN(SI 
GO  TO  30 
RR-(  SIN(SI)-RTI/(SIN(SII4RT) 

0«SORTt 1.0*1 (2.0*DtSTll*(RANGF-0t$Tll**2)/(ftANGE*H2PRIi I 
RNAG>S0RT(REAL(RR1**2*AIKAG(RRI**2) 

PHI>AfAN2( AlnAG(RR)«REAL(RRI) 

IF  ($I.GT..2)  GO  TO  SO 
OIVl.O/D 

CALCULATE  THETA 

H3«Hl-0ISTl**2/2.0 

THET A«(a.0Q01385*H3*H2<>RI*F RE 0(  ISP  11/57.  3/RANGE 
ANGLE3-THETA-PHI 

GAIN>1.0*(OIV*RNAGI**2*2.0*DIV*RHAG*COS(ANGLE3I 

DBGAIN-10*AL0G10(6AIN) 

RETURN 

ICKMP-1 

END 


N  190 
N  200 
N  210 
220 
N  230 
N  240 
ii  250 
N  26G 
N  270 
N  280 

\u 

310 
320 
330 
340 


N 


N  350 


m 

N  380 
N  390 
N  4C0 
N  410 
N  420 
N  430 
N  440 
N  450 
N  460 
N  470 
N  480 
N  490 
500 


N  510- 


10 


15 


?0 


25 


30 


35 


SUBROUTINE  ODIST  ( RANGE. HI. TALT.OISTl. ICKNP. H2PRI ) 

C 

C  CALCULATES  DISTANCE  FROM  ANTENNA  TO  REFLECTING  POINT 

C  FOR  MULTPTH 

C 

OISTl-O.O 

OH«1.414*(SORT(Hli*SORT(TALTn 
IF  (OH. LT. RANGE )  GO  TO  60 
DO  10  I*1.2p0 
0ISTl«O.l*OISTl 
IF  (OISTl.GT.RANGEI  GO  TO  60 
H3«Hl-0ISTl**2/2.0 
H2 PR  I «H3*( RANGE-0  IS  Til /DIST 1 
H2S£C»(  (RANGE-0 ISTl  1**21 /2.*0 
H2f0T*H2PRI*H2SEC 
IF  (TALT,CE.H2T0Tl  GO  TO  20 
10  CONTINUE 

20  OISTl-Ot'jTl-.l 

DO  30  I«1»1C0 
0IST1*0.01*01ST1 
H3»Hl-DISTl**2/2.0 
H2PR1«H3*  (RANGE-0 1 STD/OISTI 
H2SEC*I (RANCE-0ISTll**2l/2.O 
H2TCT«H2PRI*H2SEC 
IF  ( TALT.GE.H2T0TI  GO  TO  40 
30  CONTINUE 

40  OISTl-OISTl-.Ol 

00  50  I'l.iog 
OIST1»0.001*OIST1 
H3-Hl-0ISTl**2/2.0 
H2RRI-H3*(R4NGE-0ISI1 l/OISIl 
H2SEC»< (RANGE-OI STD **21/2.0 
HzfOT»H2PRI*H2SEC 
If  (T4LT.GE.H2TOTI  RETURN 
50  CONTINUE 

60  ICKPP-1.0 

END 


10 
20 
30 
40 
50 
60 
70 
80 
90 
ICO 
110 
120 
13C 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
2  50 
260 
2  7^ 
280 
290 
300 
C  310 
U  320 
0  330 
0  340 
J  350 
a  3bC 
0  370- 


0 

0 

□ 

0 

0 

a 

0 

Q 

U 

u 

□ 

0 

0 

0 

J 

3 

0 

0 

U 

0 

0 

0 

a 

u 

0 

a 

0 

a 

0 

0 
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LIST  OF  ABBREVIATIONS  AND  SYMBOLS 


Item 

a 

ADAP 

AFFTC 

AM 

AMS 

ARPMP 


CPU 

CW 

D 

db 

dbm 

dbv7 

DFVLR 


E  . 


F 

E 

E 


A, 

0 

A 

R 


f 

FM 


MIIZ 

GHz 

s 


Definition 

earth  radius 

Antenna  Data  Analysis  Program 

Air  Force  Flight  Test  Center 

amplitude  modulation 

airborne  monitoring  system 

Antenna  Radiation  Pattern  Measurement 
Program 

central  processing  unit 
continuous  v;ave 
divergence  factor 
decibel 
lino  loss 


Units 

statute  miles 


db 


db  above  a  milliv’att  dbm  ~  10  log  -  -  - 

db  above  a  watt  dbw  =  10  log  (P  )  -  -  - 

V/ 

German  Research  and  Test  Facility  for 
Aerodynamics  and  Space  Flight  -  -  - 

incident  electric  field  intensity  vector  volts/meter 


free  space  electric  field  intensity  vector  volts/meter 
reflected  electric  field  intensity  vector  volts/meter 


resultant  electric  field  intensity  vector  volts/meter 


frequency 

frequency  modulation 
frequency  in  megahertz 

gigahertz  (10*  cycles  per  second) 
gain  due  to  multipath  signal 
gain  of  receiving  antenna 


IJz 

Hz  X  10* 

db 

db 


253 


Item 

s 

g  (0) 
k 

L 

MHz 

n 

NATC 

PAMS 

P 

v; 

P 

r 

R 

R 

RADC 

RAPCON 

RF 

SAMTEC 

UHF 

VHF 

VSWR 


C 


e 


‘^2 


Definition 

gain  of  transmitting  antenna 
earth  gain  factor 

factor  to  correct  for  atmospheric  refrac¬ 
tion 

line  loss 

megahertz  (10®  cycles  per  second) 

index  of  refraction 

Naval  Air  Test  Center 

Precision  Antenna  Measurement  System 

power 

pov.’er  received 
reflection  coefficient 

|R| 

Rome  Air  Development  Center 
radar  approach  control 
radio  frequency 
range 

Space  and  Missile  Test  Center 
ultra  high  frequency 
very  high  frequency 
voltage  standing  wave  ratio 
permittivity 

permittivity  of  free  space  l/(36n  x  lO’) 
relative  permittivity  of  air 
relative  permittivity  of  earth  surface 
phase  change  due  to  path  length 
conductivity 


Units 

db 


db 


v'atts 

dtm 


NM 


farad/meter 

farad/meter 

degrees 

mho-meters/meters ^ 


254 


Item 


Definition 


Units 


^  phase  change  due  to  reflection 

4;  grazing  angle 

w  frequency 


degrees 
degrees 
radiem  per 


